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The cbhject of this thesis has been twofold. The 
first ckiect was to develop FORTRAN versions of 
several existing APL programs which were designed tc 
analyze univariate data. In particular tne programs 
were designed to test for exponentiality and normality 
аспа баса апа, Бу secticnirg or  jacknifing, oktain 
estimates of sampling variances оѓ ѕагүріе moments. 
The seccnd object of the thesis was tc use these 
programs in a simulation study of first-ccıe 
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inter-arrival times are exponentially distributed Lut 
dependent. The dependence is introduced ty using the 
mixed Loving average autoregresive structure 
( EAS3MA(p,q) ) for exponential sequences ilatroduced by 
Lewis and co-workers. Four models cf correlated 
queues are introduced, giving autoccrrelated and 
cross-correlated service and arrival times in varicus 
degrees. The simulation study gives a quantitative 
idea с? the effect cf correlation cn the mean waiting 
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I. INTRODUCTION 


The orject of this thesis has been twofcld. 


The first object was to develop FCEKIEAN versicns of 
several existing APL programs which were designed tc aralyze 
lu vdriate Gata, in particular to test rcr exponentiality 
and normality of the data and, by sectioning ог  jackrifing 


obtain estimates of variances of sample moments. 


The second object was to use these programs in a 
sinulaticn study of first-come first-served queues in which 
the service times and the inter-arrivel times are 
exponentially distributed but dependent. Th< dependerce is 
introduced by using the mixed moving average autoregressive 
structure  ( EARMA (p,g) ) ГОР exponential Sequences 
introduced Бу Gaver and Lewis (1977), Lawrance and Lewis 
(1977), Jacobs and Lewis (1977) and Lawrence and Lewis 
(1978). Four Schemes are intzedueed which give 
autocorrelated and cross-ccrrelated service and arrival 


times in various degrees. 


Since nc analytical properties of these queues can be 
derived, their properties are studied ry simulation. 
Because the EARMA queues are neither regenerative nor 
Markovian there are several froblems in the simulaticu study 
of knowing when one can assume that a steady state has teen 
reached in the simulation and of cbtaining confidence 
intervals for estimated parameters. This prcklem is handled 
Ey generating a number of independent sample paths and 
 Шратіпа statistics obtained from the treplicaticrs at 


Several pcints along the sample paths. Only the waiting 
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time prccess k; НО ШО еса this thesis, not the state 


space prccess N(t). 


Unfortunately the simulaticn study сі the EARMA-type 
queue was hampered by the size and speed СЕ tne 
IBM System 20/67 system at the Naval Postgraduate School. 
Consequently the results were not as extensive as had teen 
Moped fcr. Fairly detailed results are given for the case 
where the service time process is autoregressive; one case 
Шеге tke traffic intensity is 0.99 and the correlation 
parameter 0.98, was intensively investigated. This gives 
scme idea of the length o£ the transient in the queue and of 
tne inflation of the mean waiting time Lrecause cf the 
correlaticn іп the service prccess. Some investigation of 
the case where the service times and inter-arrival times are 


cross-correlated are given. 


In section II we describe the basic histograı and 
EEUIStIlCcSs computation package  'HISIGS/F' and apply it to 
tao sets of telephone errcr data. мН! сессісоп IJ ве 
erting method and the plotting Subrcutines  'NOFMEI', 
'EXPLT' аге described and we test them with various 
generated data. In section IV we describe the assessment of 
МЕТ ра 11 Ку Suproutines 'SZECTN', 'JACK' and we apply ttex to 
telephone data set 1. In section V the Surbrcutine 'IIST' is 
described and we apply it to two sets of telephone error 
data. In section VI we describe and analyze the EARMA(E,9) 
model as well the program wazlch че developed in order to 
simulate it. Computer program listings are provided after 


Section VI. 
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11. ВИЛЕН Це вз 165/81 5125 


сэвтэ ae ош» жс чыъ LLL ALII | SEE ca s s NEM 


DESCRIPTION 


Here the Sanroutıme "HISTGS/HISTFS' 15 presented. 
ENSSICS/HISTFS'!' is used for estimating the protatility 
density function from a given set of data and at the same 
pie computing some kasic statistics. Easically this 
Se@rOoutane is the library FORTRAN-Subroutine HISTF/G, which 
was created at N.P.S. by 2. И. Robinson. It has teen 
шодатлед Еу adding the new entry point 'LIMIIS'. 


The new entry point may be used for wild data ог for 
data having a mixed distributicn. Therefore it is a useful 
Peer tor that user invcived with this kind cf data. Simply, 
the user has to define those values in the range of data 
(sectioning the data into several disjcint sections) that 
he/she Ltelieves are useful cut points. Fer each secticn of 
data the user may have the prcbability density functicn as 
well as the basic statistics ( as HISTF/G ¿ces ). Alsc Dy 
sectioning the data over its range, the user may have a 
histogram cf the random variable when it is conditioned to 


Le between given limits. 


Пс ошшке: о sections that the entry pcint 'LIMITS' can 


accept is as many as 5C. 
A ccmplere description of how 'HISTGS/HISTFS' operatss 


is given in the subroutine. On the other hand, a summery of 


the subrcutine is given by typing on the terminal the 
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Шап DESCRIBE HISTGS under CMS (CamEbridce Monitor Systen, 
Ref. 4). Ia DESCRIBE HISTGS' the following response 
will be given on the terminal: 


SUBRCUTINE HISTGSZHISTES 


' HISTCGS/BISTFS' gives you a histogram cf a set of data 
along with the estimate of a set of basic statistics. Also 
ШЕ д4іуес you the fcllowing options: 

1. Displaying a smoothed empirical prctability density 
function cver the histcgran, 

ЖЕ Sectioning the data into several disjoint sections 
Mo more than 50) and then having a histcgram anc the 
estimate cf tasic statistics for each secticn, 

Би Scaling the histogram, 

4. Displaying just the histogram (witt or withett the 
Beneity functicn) and no statistics. 


The estimated basic statistics are: Mean, Median, 
Irimean, Midrange, Geometric Mean, Harmcnic Mean, Variance, 
Standard БЕ вава сп, Coefficient et Variation, Mean 
Bwsation, Range, Midsgread, Third and Forth Central Mcment, 
Skewness, Kurtosis, Betal, Beta2, Maximum, Minimum, and 


Cuantiles. 


CALLING SEQUENCES 
MEL — HISIGS ( X, N, NBAR ) 
CALL HISTFS (X, N, NBAR ) 
myer LIMITS ( X, N, XLIM, LIMS ) 
CALL . NCSIAS 
BL STATS 
IL 21Х5 (SCALE) 
CALL  NOFIXS 
CALL РЕСМАХ ( PSC ) 
CALL NCERMX 
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ARGUMENTS 
X àrray of data values 
N Number of data values 
NEAR Number of bars in the histogram 
XLIM Array of cut points 
LIMS. Number of cut pcints 
SCALE Vectcr of two values to scale the data 
PSC A real variable (between 0. апа 1.) defining 


the maximum value of the probability axis. 


шенЕ пе entry point NOSTAS 15 called befcre callinc any 
МЭ nc entry points HISTFS,  HISTGS, or LIMITS tken no 
statistics are printed, otherwise the statistics are printed 
Ey default. 


By calling the entry point FIXS the histcgram scale may 
be fixed and remains set unless it is reset Ey another cali 
BUNPIXS or allowed to float (the default) by calling NCFIXS. 
Ihe reason for fixing the scales is so that when ccmraring 
cre than one batch of data, we can get сспракарју scale 


histogrars. 
By calling PROMAX, the maximum value of the prokakility 
E can be set to PSC, which will rémain in effect 


ШІСІ! NOPRMX is called (the default). 


More infcrmaticn is given in tne suproutine. 


USING 'HISTGS/HISTFS' WITH TELEPHONE DATA I AND 
ШЕБЕРНОМЕ DATA 2 


Here 'HISIGS/HISTFS' has been used to plct the set of 
Telephone Data 1 and Telephone Data 2. The algorithmic 
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procedure fcr both sets of data is the same. That is, the 
EN »ernt  'HISTFS' is called first and then the entry 
puc “LIMITS'. Four cut points have been used (1, 2, 141, 


86000) sectioning the data into three disjcint secticns as 


fcllows: 
Sect1ca I: Consists of data points i such that 
i) 
ES i 2220 
(1) 
Eection 2 : Includes the data pcints Да по пе пат 
1? 
PS X < 141. 
Ci) 
Secticn 3 : Includes the data points 42 such that 
i 


142 < X < 86000. 
(3) 


No grint-out for secticn 1 is given because cf the 


ccnstant value of tne data pcints. 


duties 1 through 6 give the histogram and the Easic 
statistics fcr both sets of data ( For the entire set ала 


Mmeeesecticns 2 and 3). 


We observe that both of data sets are sc unruly that the 
Mee OF tke entry pcint ‘LIMITS! cf ‘HISTGS' 1S suggested. 


Observing Figures 1and 4 we may conclude that botk data 
sets appear to Ге the same, having =сіє geometric or 
expcnential distributicn. But as we will see later this is 
men true. Thus the histogram, in that case, which is 
merarned by using just the entry point  'HISTFS' is not 


helpful and leads us tc a wrcng conclusion. 


'LIMITS' is applied now to give us a useful answer : 

Observing figures 2 and 5, obtained by using the entry 
Ent '*LIMITS', we may conclude that the data sets dc not 
have the same distribution ( at least in that interval ). 


Iü addition data set 2 is far from having a geometric cr ал 
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expcnential distributicn, as it can be seen by the histogram 
of the entire set of data. ла” а cencdusion "which 
cannot Ее cktained Бу without use of the entry pcint 
ВЕ ВМТ '. 


Observing also Figures 3 and 6, where tke histcgram of 
the 3rd secticn frcm each data set is presented, we may 
conclude that this section o£ both data sets шау have an 


exponential distributicn, but not the same спе. 


Maem nost important facet of the data which appears is 
tne modes at apout 124. The mode is much acre apparert in 
Telepnone Data 2. In fact closer investigaticn reveals that 
this mode was due to leakage of dial puises; this leakage 


caused bit errcrs at multiples cf 124 bits. 


The above is а quick and informal analysis of the 
Telephone Data Sets obtained by using the entry pcint 
ДКТ of 'HISTGS/HISTFS', and the conclusicn is that koth 
Паста sets are so unruly that any formal characterization 
БІ their distribution is not so easily obtained. ES 
addition it is apparent that a mixed distribution gcverns 


meth sets of data. 


Having quickly analyzed bcth of the telephone data sets, 
u use of entry pcint  'LIMITS' of the SUbEGUtInE 
EIU OTGS/HISTIFS' is cbvious. 
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III. PLOTTING SUBROUTINES 


А. INTROLCUCTION 


Graphical methods of assessing the validity c£ a 
probability model, and of estimating some kasic parameters 
(especially lccation and scale parameters) have been widely 
applied over Many years. A Variety of different 
prescriptions have been advanced for tne plctting positicns, 
that is, fcr the set of values at which the crdered 


cbservaticns ina sample should te plotted. 


Plotting the crdered data is ап infcrmal and quick 


Visual methcd for getting an impressicn of how the 
Ae tripbution of data looks. Depending сп what is kncwn 
atout the distributicn of the data there are many 


preobability plotting methods. 


The Flotting Subroutines presented here,  'NORAMPL' and 
APLI, use the methcd cf plotting the crder statistics 


( E ) versus their expected value ( о ) called 
Ш CORES . This methcd does not require knowledge cf the 
emmu lative Cistributicn Functicn F(x) (except for its 
assumed continuity). Therefore it may ke very attractive 
Mees chose whe kncw ncthing about the data tc te analyzed. In 
addition the indication of the graph may lead the analyst to 


further useful results. 
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Using the Plotting Subroutines 'NORMPL* and 'EXELTI' one 
may analyze the resulting graphs obtained tc get an idea of 


ee distribution of the data. 


7 


CAD ESCRIPTICN OF THE PLOTTING METHOD 


The idea behind this methcd, as abcve stated, is tc flor 


X versus Ё(Х j: 
Ci) Ci) 


It is used where data ч 2: ЕХ arise as 
n 


independent observations о: а continuous Bandom Variable X 


a3 distribution function which is believed to have  scme 


particular form Еб). Then the ordered data 


x Х я аге plctted against Ef X : 
E О со“ E P 3 L Т. 


It the model for F(x) used to get BCX. ] is true then 
1 
the resulting plot is a linear (regression) relationship and 


ЕНІП plot, within limits of sample variations, asa 


Eragmght line. Barnett [2] gives the value cf ЕХ yd as a 
i 
ШПсат Celaticn of the plotting position Ey" Tuscos. 
i 


EX 1 5 mnm + sx 2 
Ci) (ji) 


here m, and s, are the Location and scale parameters. The 
Mean and Variance of X are nct necessary. As we see the 


above relaticn uniquely defines the plotting position as 
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ВХ. то) 
(1) 


X = ------------- 


Ci) S 


Ci) Е 
ши сет 


ШШЕЕЕТОГЕ, aif the model is correct the linearity will 
 Баріу їе evident. In addition if the resulting plot not 
being linear leads us to reject the model, we may cet an 
шПісасісп сі «һе type cf departure present and thus ап 
indicaticn of what action tc take or what alternative mcdel 
to use. It should be noted that the linearity of the graph 
suggests informally that the model is true, but rct tc 
accept the ncdel. Additional formal tests must be done іп 


Eder to accept the model. 


ЦИА SUBRCUTINE NCRMEL 


SUDECUtine  'NCRMPL', is used tc test a samrle of 
independent cEservations of a ccntinuous Random Variarie Х 
having a distribution functicn which is believed to have the 


Ncrmal fcrm. 


'NORKEL'! sorts the data into increasing orcer to 

Сгеа<е а set ci Order Statistics. Then the Crder Statistics 

are pictted versus either Ncrmal Scores 

F ns 50, n - number of data-points) or irverse of the 
-1 

Standard Normal Cumulative Distributicn Function (F 77, 


ШЕ п > 50). 


ЭЙ 





Evaluation of Normal Scores is based on ccaruting 


th 
the Expected Value of the i Under statistic БГ |e 


Ci) 
This evaluation follows: 

Let 
| | nc | 
- be the density function cf the i Crder 

1 
Statistic of a sample 2 i. Unter X fror a population 
n 


ШОШО the Standard Normal Listribution Е. Then ча) 15 
I 


given by : 


п! 421 Ша Да 
Хэ.) 5 -------------- Ғ (x, -F (x )) (O 
Ж E! (сг тг) X со) | x Ci) X. (i) 
n! ша ост 
з -------------- Q Z 
"sp 2-1)! ' 
where: 
: 2 
1 А су / 
2 = Е (х A ечи оде. со сто она e , 
Х (19 1/2 
(2 n) 
Ха) 2 
1 -u 2 
ENS E #2 Б оо = d j 
201) 1/2 
200 (21M) 
Q = 1-E 
Ncw 
oo 
EI X МЕС ах 
(3) 
“од 
n! dd i-1 n-i 
(1) Zee) X 24 P Q dx 
(1-1)! (п-1)! 
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The value of expression (1) has been evaluated and 
is given in [19]. This table gives the values of Normal 
Scores frcm 1 to n/2 (for nS50). The values from (n/2)*1! to 


n are evaluated using the symmetric expressicn : 


EIX = - ЕХ 
( TER ` E 
Also 
E( X Е D= 2К®1 _, 
(x+1) 
"МОБМЕГ uses the  arove Normal». sSccres 1: п < 50. For 


n > 50 !'!'NOR3EL'! computes its сып Normal  Sccres using the 


program function  'INVNRM', which gives the value x cf the 


inverse Standard Ncrmal алс Ор (x = F (у), 
ШИЕ үу 5 1). Мен са лом. ОЕ this is the fact that if the 
data are really normally dreeributoed then 


EEX J Ер 2(1-1) у. Therefore 'NORMPL* can accoıcdate 
1 
any sample size. 


ОИЕ аас опи to the plot given ty 'NORMEL*, the 
Nilk-Shapiro test for ncrmality can be evaluated by the sane 
ЕГО де. A detailed description of that test is giver in 
(22]- However a summary of that test is presented here in 
crder tc help us to describe how the frcgram works. The 


Wilk-Shapiro test is based on the statistic: 


2 
b 
и = Er 
S 
where: 
b = X X ) К 2 
= a - > =т/ Dec 

1=1 52 (0-14 1) (1) ‘ 
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t a (X -X 125 K=(n= 13727 п ова 


n 2 
BE» (X - X) 
121 T 


The coefficients (a ) are defined by: 
ae 


СИ и ate ZC; л qan 
+ an 
where: 
m = Е[Х 112 (That is che TNORMAL SCCRES!) 
j (3) 
u = Coví X X m ИЕ e 
j. L Ci) съ?" , Ј , 7 , 
Z  ' -1-1 


where: 
Ш = (2 ,a 201 
17 pi а 


That 2s the nX¥n Covariance Matrix. 


< 
Il 
=== 
с 
е” 
~ 


There are some approximations asscciated with the 
evaluaticn of these coefficiencs. However fcr n < 5C the 
(а ) are tabulated in the above-mentioned paper ard the 

l 
program uses them by FCRTRAN DATA initializaticn staterents. 
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ЁСГ:  n»50 the (a) are computed by the program using 


1 
the f£ollcwing approximation method: 
a = "mm e, ral, у 
1 ql 


where C is given by the following Least-Square equation as a 
mumetion cf n: 


2 
C 2. 122 + 4.0930 : 
Ec. for n=1 
2 (1222) 
а 1⁄2 АУ EET . 
2 NGA 21) 
Ncte that а = а Е 
д nec 


ЕСШЕ properties for the W Statistic can be given 
here, which have been taken from the above mentioned paper: 


(a). ч 15 scale and origin invariant. 


(b). 


и has a distribution which depends only on the 
sample size n, for samples from a Normal Distribution. 


^ 


< 
(с). Я is statistically independent of S and of 3 for 
samples from a Normal Distribution. 


Em ZI DI 
ня | = ЕБ EIS |] for any г. 
(е). The maximum value of W is 1. 


(1). Тһе zinimum value of W is na / (п- 1). 


I Ihe half and first roments of W are given by: 


2 
е Ru E92 55) 
E[ W ] E БОСА ЫН – : 
C Г(аи2) 2 
РОА 
ee 
C nci 
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where: 


2 t =] 
R =m Ç nm 


(h). For n-3 the density of W is given as: 


E 2-172 
£ (w) » -z- (1-w) ? ч 2 2215 


IA 
x 
IA 
EX 
. 


PE Program Structure 


олы ae zn u a a A En ee En a en 


'NORMPL' is a FORTRAN-callable Subroutine with єаса 
НЭГ returning (optionally) either a plct fcr a given set of 


data or a value for the Y Statistic or both. 


Me ргостая 15 divided into three farts: 

дете part 15 the cortrol prcgrad of 'NOFMEL'. 
In that gart the "Normal Scores: for n = 2(1)25 апа 
n = 26(2)50 are stored along with the Sukrcutine Ey CATA 
Statements; on the other hand, the Normal Sccres for r > 50 
EX conputed by calling the prcgram's functicn 'INVNRM'. Іп 
this part the user specifies “ust a plot cr just the value 
НИША ог beth. The Library-Subroutine 'PXSOERT' is used by 


the progran. 


The second part of the program ccısists cf the 
EuDroutine  'PLOT' which accepts any data set to be tested 
екот дај)1 Су.  'PLOT' itself scales the data accordirg to 
Bs crange and plots the scaled data along 110 equal spaced 
Mesttions of the X-axis. No plot is given if the data has a 


constant value. 


IE cg Id Patt of the program 15 the INTER! 
Subroutine. ШІ Іс Subroutine the ccefficients (a) 
1 


ШЕЕ п < 50) of the W Statistic are stored using FCETRAN 
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EN n:toralyzdtdon Statements. For n»50 the coefficients 


ate computed by the Subroutine itself.  Ihis Subrcutine 


HE calls the 'INVNRM' Function to get the 'Normal Sccres' 
Ж Ж БЕ псевйб ісг the ccurgutaticn c£ the (a )'s when n > 50. 
1 


A ccmpiete description of how 'NOKMEL* operates is 
given in the Subroutine. However a summary is obtained by 
typing on the terminal DESCRIBE NORMPL. When the user types 
the command DESCRIBE NCRMPL under the CMS environment the 


fcllowing response is printed on the terminal: 
BUBRGUTINE NCRMPL 


' NOR MEL' takes a set of data, Sorts it into 
ascending order and uses the created Order Statistics for: 


EN Plotcing x versus either ‘Normal Scores" (if 
(1) 


50, п = Sample size) or Inverse cf the unit  Ncrmal 


IA 


n 


Samulative Distribution (that is: 


1 
X i) versus Normal scores or F и (пъ) 52) 
(1 
to see if there is a linear fit. 
МЕН Computing the value of the W test сат ее 


EN 5һарітгс test for normality). 


HU IS call como: 
CALL КСЕМЕГ IAS ABRES Н», К) 


ARGUMENTS 
X IS the array containing the data 
SCORES Is a work array of dimension N 


Is the number of data values 


ж ш 


User's option 


Bor; К = 1 Ao опу, 


on 





| 
м 


The W-value is given , 


| 
UJ 


Bown or the above are given. 


Mcre infcrmaticn is given in the sukrcutine. 
3. Interpreting Ihe Output 


If data are really Ncrmally distributed tten a 
ООС ас line is expected from the output. This is an 
indicaticn that gives the user a first feeling of the 
Ши глосш:јоп of data. ши the plot is not lingar then 
the user has to reject the normality  assumpticn. In 
addition tne shape of the plot пау lead the user to make an 


alternative decision of what model he may use. For example: 


(д) Эпаре СТІ the form of Figure 7a suggests thar 
data are skewed to one side (specifically pcsitive values 


lu 
In that case an Expcnential zodel could ke an appropriate 


one. The histogram of the data is given in Figure 7r. 


(5). Shape of the form of Figure да can be 
interpreted as data not as dispersed as the Ncrmal 
Шише стари тт1оп. This set of data may have a Unifcrm cr a 


апап ат Distribution, and could, for instance, arise if 
measuring tolerances or components hand-picked to lie within 


tolerances. The asscciated histogram is given in Figure 8b. 


(GE). Shape of the form of Figure 9a indicates much 
шсге dispersed data than Normal, but the data is symmetric, 
EroEabl;y Semen Gg наст оо Е “exponential ос Cauchy 
Distribution. The associated histogram is given in Figure 
Bh. 


It should be ncted here that a straight line does 


nct always arise frcm Normal cata. Figure 10a for example 


38 





ес ап indication that data are Normal (Straight lire). 
But these дата have been generated by computer from a 
Symmetric Triangular Distributicn in the range (0 to 2). On 
ЕНЕ Other hand figure tia plots a data set cf size 50 again 
meemea Triangular (0 tc 2) distribution. Eut in this case 
we observe that the graph is not a straigtt line. These 
results dsmcensctrate the informality of that  Probarility 
Flotting Method and in these cases the user has to ccrtinue 
testing the data by other formal tests. Nevertheless, since 
this method is primarily an initial screening device, ease 
Application is important. Figure lir gives the 


associated histogram. 
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4. Using 'NOSMPL' with Generated Data rom Varicus 





uuetributions 


The Subroutine 'NORMEL' has been used here tc test 
for Normality some data sets which were generated Еу tne 
computer frcmz varicus distrikutions. The plots of these 
data sets illustrate the advantages and disadvantages cf the 
HNcrSng Subroutine 'NCRMPL'. 


CEserving Figures 7a-b, 8a-b, 9a-tb and 12а-Е we may 
meee the follcwing: 


(a). Figure 12a shows a plot of a sample of 5€ data 
MES frcm a Normal Distribution (N(1000,1) ). As it was 
expected a very straight line is fitted. In additicn to 
Mere plot, the Subroutine  'NORMPL' uses the opticr to 
evaluate the ‘WILK-SHAEIRO' test and the W-value is ccıpfuted 
and printed cn the same figure. Comparing this value with 
the percentage points cf the W-test which are given in [22], 
for sample size n=50 we see that we may accept the normality 
assumpticn with a significanse level zx c кти е 125 


gives the associated histogram. 


(eye she lgure 7а plcts a samrle cf 200 data rcints 
MS Unit Expecnential distribution. Here the nonlinear 
graph is otvious (as is expected ) and suggests tc 
reject the ncrmality assumpticn of the data without doing 


any further formal test. 


о осе ба gives a plot ci a sample cf size 
200 from Uniform О Го) generated data. The shape of 
that plot indicates data not as dispersed as the ncreal, 
therefore the user is told to reject normality, getting at 


the same time a suggestion that data may be Uniformly 
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Ectributed ( because of that particular shape of the 


graph ). 


(d). Figure 9a елсесесте strongly departure frcm a 
Normality assumpticn. Indeed this set ог data have teen 
generated from a Cauchy Distribution and its shape really 


indicates a Symmetric distribution. 
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INORMPL!' was used with the ccst cverruns data (61, 
шог седе year 1950 to see if there is a linear fit.  Lccking 
at Figure 13 we see that data шау be Normally distribute 
Since the plct is very close to linear. Of ccurse 5сшє data 
points deviate from the straight line, but these points may 
be considered as outliers. 
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IEEE SUBRCUTINE ЕХРЈТ 


ШІ Description 


o Outin EXPLI Tis andended tc be usec for 
testing a sauple of independent cbservations cof a continuous 
ош Variable X with a distribution function which is 
relieved to Ee Exponential, by clotting each of the Order 


Statistics X versus its Expected Value, ard by estimating 
(1) 
the parameters 4, and 1, Of the observaticnas ` ( ML and E 


~ 


Encroduced later). 


This Subroutine sorts the data into increasing 
emer, ortaining the Order Statistics. The Order Statistics 
are then fletted versus Exponential Scores. These scores are 
evaluated by the Subroutine. The Exponential Scores are tne 
Sepected value cf the Crder Statistics ( Ех 1). They 
are more easily derived and computed than the one Scores. 
That is: 


Let X be a random variable Expcnentially distributed 


with parameter A= 1 


Enon it is krcwn: 


DUE AC , ar x? 0, 


И 
ша 
Í 
(D 

~ 


20 up coU. 





Ë m 
Ihe Density cf the i (одоше зса сл ес is given by: 


Е _ n-1 - 1-1 И n-i 
Е (хз) и Ло ) А. ( по О, 


st 
Thus for the 1 Order Statistic we have: 


(1) n-1 (1) 
É = 1- (1- 
E ТО; niue ) ) ж 
(1) 
Е 
(1) 
= пде : 
Therefore ar is exponentially distributed with 
parameter > = n x 
And for %= 1 we have: 
1 1 
ЕХ = --- = --- 
( a)? n X n 


It can be croved (see Feller [7]) that the n random 


ER ыг mx -Х -—— MER -X from 
(1) (2) (1) (3) (2) (n) (n-1) 


an Exponential random variable X with density functicn: 


variables X 


- 2 
PIN RE e : x > 0 , 
Х 
are independent and the density of the random variable 


Y = X mi is given Dy: 
(К) (k- 1) 


AMARE y 
Eo шил uns , 1 


That is, Y is an Exponential randca variable with 
parameter 20251122 
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Méetetore E[Yj is given by: 


1 1 
E[ X E ELS o >< ------- , For il де 
(К) (К= 1) A(n-k*1) n-k*1 
Ihis fact is used to derive е Заир: = 
1 
That is: 
mex ) сег (X — ) + (x ` Е 2-7 
(1) (1) (1-1) (1-1) (1-2) 
(X EX ) * X ] 
(2) (1) (1) 
= ЕХ 1 01:12:11 =x Шэн 
(1) (2 (1) 
E -X 1+ EEX -X 
(1-1) (1-2) (1) (1-1) 
1 1 1 1 
+ --- + ------- + ------- + . + ———-—--- 
n noo | n=3#] пеещ 
: 1 
= Ж, Е. 7-2 27227! . 
к=] п-к+1 
'EXPIT' derives the values of ЕГХ i) J usirq the 
(1 
осе formula and then it uses its Subroutine 'EPLCT' to 


pict the Order Statistics versus tneir expected value. 


ЕНТО to the plot, which 35. ап informal and 
ema ck test шог exponentiality the estimates of the 
statistics y and (2 provide another informal test since the 
values СЕ 1, and n ha ve a constant value ог апу 


EE cuential distribution: 
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Ey definitior ), апа ps are given by: 


(X y) 
mi 
3, = Бин! , 
О 
(X TN 
т 
Коо з. 
6 
Ihus 
1 5 
ши 3 Е[ (Х-ү) ] 
` D С ax EX Ч | 
= --- = А + - 
27 | үтү 
; 1 b ы 2 5 2 3 
(1) "ccs (X ]- 3gE[X ] * 39 8(3] Су). 
O . 
Now 
3 7 3 
мк Ја f (x) dx 
2 X 
oS - АХ 
- \ x Ne dx 
қ. 4-1  -лх 
[ (4) E A е dx 
88 5 ГІШ 
Ihe above integral is that of the Gamma density Ешре стон 
with parameter QO = 4, and since its value equals 1 че 
have: 
га 3! 9) 
ЕХ = ---- = ---- 
3 5 
A A 


ЭЭ 





2 
Ccaputing E[X ] we nave: 


2 “О г ДЕ 
E[X ] = | X ле dx 
о 
muc З 20 
Е [ (3) ) Х A е dx 
42 - СО) 
_ 2 
22 
E ЕЕ ТЕ 
Alsc E(X] = > Пе : 


Substituting in (1) we have: 


x КЕ 
(<= 
= Е ча» «и <= «к» <= = = = 
б, = El E 
1 6 1 2 1 1 1 
= — — = ( = ~ = = — 3 = = => мо E 3 — | — “= – = <= = = = = = ) 
237 Der ог. 
ни) со 


Iherefore tne value of ү, equals 2 for any „Ехрспеп за! 


E tributicn. 
Alsc A шанс сошриетас аз follows: 
ü 


Е. 
ЕГ E 1-3 


da 


4 3 202 3 4 4 
ЕХ -4үХ «бүх - 4UX + y] VR 13 


|| 


ü 3 2 2 3 4 4 
(2) 12 = ay EL X JEt ONU EL X а а Zu cd. 
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put 


4 < Ц = \X 
ЕХ ] = ) ии де дх 
Q 
"X О=] 5 ах 
ШӘ) J x A e dx 
M > ( (5) 
4 
= 24 / A . 
Thus 
à 1 í 25 Ë 1 6 Р 1 : 1 1 1 
EE ,_--. - 21-14 мэ... + ---} - 3 
2. ~ 
с > AU >. м л M 


Iherefore ЕВЕ value of un equals 6 for any 


EU cpoential Gistributicn. 


ШЕ Басс 1с used by the Subroutine 'EXPLT' tc give 
the user cne more informal test ror  exrcnentiality. Ir 
estinates the values of y, and Y, and ie 
ШО Ос \,* € then the data may not have an Exponential 


AS tribution. 


Меша виша the values or j; and ү,, "ЕХРІТ" uses 
y 
unbiased estimator fcr c, El (X-p) J and E[ (Xy) ] using 


the formulas: 


n 277: 1222 
sum (X--X) / (n-1) ) , 
qo I 
n "EE 
n (X = X) 
E[ х= * zen EM Веи 
I (n-1) (n-2) ‘ 
3 Я D T 
-2 с у еу 
Ч (И на) E m ) 3С (0-1) (20-3) 
МИ |) IN NEL, 22 22--------------- - ----------------- 
Ше (1-2) (п-3) п к= 2) mes) 
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> ЕР БИЙР сэ а ie занар a s s Fm ш> 


F ОБОИ Саше 2 ЕХЕЇТ" is also a FCRTRAN callable 
ENMoutine with each call producing a plot fcr the given set 
Of data and at the same time calculating and printirg the 


estimated values of the parameters y, and y, of the data. 
Basically the program is divided intc two parts: 


The first part is involved with all required 
@emeutaticns calculating the Exponential Scores ard the 
various statistics needed ЁС compute the 11 and 1, 


parameters. 


ШОЕ = =сова part is the Subroutine 'EELOT'. PERTOT! 
Scales tne data points according to their range and plots 
tne scaled deta along 110  equal-spaced pcsitions of the 
X-axis. No plot is given if the data have ccnstant value or 


if there exists a data point less than zerc. 


A complete description of how 'EXELT' operates is 
given Dy the Subroutine and a summary description is given 
сп the terminal by typing the command DESCRIEE EXPLT under 
the CMS environment. When the user types DESCRIBE EXFLT 


the following response is printed on the terminal: 
SUBECUTINE EXPLI 
Шери“  aKes а set of data,  scrts at Into 
increasing crder and plots the created crder statistics 


versus  expcnential scores (expected values of crder 


Statistics) computed by : 
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i 
ЕГХ ] = = ( i/(n-k+1) ), Е м 
k=] 


nct computes the estimates of 7 апа ын 


parameters. 


pedis called Dy : 
ШЕГІ EXELT ( X, SCORES, N) 


ARGUMENTS 
X Is the array containing the data 
SCORES Is a work array of dimension N 


N Ts the number of data values 


Mere informaticn is given in the sukroutine. 


ЕЕЕ Іпбсегргесіпа Тіс Output 


сэржэ = s шю =» we Se meum тз 


Tí data have an Exponential distribution then the 
Hot will tend to a straight line. But it should Le ncted 
that the inverse is not generally true. The linearity cf tne 
graph is an indication only and gives the user a first 
Meeting of the distribution cf data. If the plot is not 
linear, however, this suggests chat they are not 
Meeponential. In addition the shape o£ the plot may lead the 


user to make an alternative model selection. 


Cn the other hand the values ер: the 
statistics 4, and 4, may be used to test informally the 
Messer ibution of data. These values are an indicaticn cnaly 
EU rhe test for expcnentiality. If the estimate of à, has 
the value abcut 2 and the estimate of Y, has the value about 
6 then TE suggests that data may Ге exponentially 
distributed. But a departure of these values suggests that 


data do nct have an exponential distribution. 
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. Using ‘EXPLT' with Data 


The Subroutine "EXPLI' has been used here,  plctting 
data generated ру the computer from varicus distributions, 


Eugcrder to get a visual sense of the behavicr of the  plct. 


Figure 14a gives the plot of 10€ Uniform (0,1) 
variates.  Crviously ncnlinearity governs the plot (as we 
should expect) and thus it suggests  derarture frcm the 
Exponential distributicn. Besides the nonlinearity cf tne 
NES the estimates of j,and Y, (0.1 and -1.15 respectively) 
are far away from the ccrresfending values (2 and 6) of 
exponentially distributed data. Figure 14b gives the 


associated histogran. 


Figure 15a plots variates generated frcm a 
iN ngular Symmetric distribution (0,2). The nonlinearity 
Ши the values of the estimates Y, and = пасе оо 


reject the exponential assuurtion. 


2 
Figure 16 shows the plot of X generated variates 


with 10 degrees of freedom. For tae same reasons as  abcve 


the exponential model is rejected. 


Figure 17a plots Normally  (N(1000,1)) distributed 
generated data and the nonlinearity of the plct is apparent. 
Eesides the nonlinearity, the estimates of Y, and Y, 
ОЗ апа -0.63) suggest not only ЕС reject tne 
ENDonentiality of data but in addition lead us to test the 
data о хогта11:у. Figure 175 gives the associated 


histogram. 
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Figure 18 plots expcnentıally generated variates 
with parameter A= 1. As we see the graph is linear ané the 
user gets the idea that the data possibly are exponential. 
Besides the linearity here, the values cr the estimates of 
Мада у, (2.06 and 5.4) are close enough to their true 


values tc indicate an Exponential model. 
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ШЕ СПЕНКОШТІНІЕ LIST 


ИІ DESCRIFETICN 


Anotber Subroutine presented іп this thesis is the 
EN uUrIne 'LIIST'. Its purpose is to list a given set of 
data in ascending order and, taking advantage of like 
 Пггепсес іп +һе data, to print the ordered data in a 
compressed form. This feature becomes highly useful when 
listing a large number cf data points that contain папу 
repeated values. iig is also аА ес АЕ се но Ба еј] 
ШЕР 1с151265 іп supposedly continucus data, and a 
probability function estimating routine fcr data which is 


known to Ee discrete. 


A complete description cf how 'LIST' cperates is given 
in the Subrottine. However a Summary is given by typing 
ESSCHIBE LIST. When the user types the ccammand 
DESCRIBE LIST under the CMS environment the follcwing 


response will be given on the terminal: 


EUENCUTINE 1151 


Mest’ sorts a set of data into increasing order and 
= а 3-cclumn print-out as follows: 
ie Column: Serial number cf the first occurrence of 
this value in the crdered List 
2nd column: Value of ordered data-value 


БТС її: Frequency cf cccurrence cf the value 
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ШОО column: Percent for the value 
шэг сс1ї:: Graphical representation cf сое Frequency 


mer each value. 


ie as called by: 


CALL BIST. (X, N) 
where: 
X Is the array of data 


N Is the numcer of data-values. 


More xnfcrmation if given in the subroutine. 


Mee UNTERERETING THE OUTPUT 


The rrint-out is a visual representaticn cr the data and 
each data value frequency. For data рсіпіѕ havirg the 
same value 'IIST!' will print this value once with the rumber 
СЕ  cccurrences. AwNS-coblumm output will te rrinted ard its 
iaterpretaticn is: 

ШОО First column gives the serial numter of tne first 
cccurrence of this value in the crdered List 

2. Second column gives the value ci ordered data-vaiue 

3. Third column gives the frequency cf occurrence of 
the value 

ИЕ ҒОГЕП column gives the percent for the value 

ЕЕ: Fifth column is a graphical representation cf the 


frequency for each value. 


a} 





ШОО LS5 54 2, 2, 5, 6 bea given data set. 


Then the crdered data will be: X =X =x = 
(ly (05) 


6. Then 


Mme = 2,3 =3, I =X = TES 
(4) (5) (6) (7) (0206) (10) 


the print-out will have the follcwing form: 


s D NUMBER ORDERED DATA FREQUENCY PERCENT РНОБ. СБАРН 


= ИО G “a a O A e A Á s “s m ЕШ mam m° sam eee eS —— s s ÀO “h me m mas sm O V GU RUD TUNE ЕР TO s À Ó s m - - 


1 1 3 T oro 
4 2 2 22 ++ 
6 3 1 01 : 
7 5 3 23 ... 
10 6 1 3 + 


If there are no data-rcints having a common value then 
2-1" gives this indication and prints crly the ordered 
data. Itis happens when da ta have a continuous 


Ae tribution. 


C. В LISTE" WITH TELZEEONE DATA 1 AND TELEPHONE ГАТА 2 


ШО ОШО Пе LIST! was used With Telephone Data 1 
(Figure 19) and Telephone Data 2 (Figure 20) and a brief 


matysis сі the output follows. 


ШООК гд at Figure 19 and Figure 20 we may cet some 
informaticn cf each data set. It can be seen that toth of 
the data sets contain a large number of multiple occurrences 
cf the data value cne and data value two. As we can see 
reading the fourth column the occurrence cf cnes is 1°% for 


Telephone Data 1 and 24% for Telepnone Data 2. 


78 





Comparing also Figures 19а-с, 20a-d we can see that 
Mitiple cccurrences happen in the range 1 to 24C for 
Telephone Data 1 and in the range 1 tc 132 for Telephcne 
Mata 2. 


Furtheracre, quick visual informaticn concerning the 
rangs where we have multiple  cccurrences can be cbtained 
from the Prcbability Graph. Thus we can see that for 
Telephone Data 2 there is a regicn from 113 tc 132 where the 
multiple cccurrences of values is larger than in neightoring 
Mmeemons.  Thererore 'LIST' gives the user a useful tccl fer 


analysis cf data. 
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À. INT BRCDUCTION 


In the previous Sections Prcecbability Бест С 
Subroutines were presented fcr use in informal estimaticn of 
the form of the distribution of a set cf independent 
cbservaticas. In the present section subrcutines computing 


the variability or some basic statistics are described. 


4e kncw that X (sample mean) is an unbiased estimatcr of 


the populaticn mean  ( = x f (x) dx), where X is a 
Х 

Memputaticn with unkncwn distribution с) ) and urxncwn 

Earum and X , X , ..., X is a sample cf independent 


1 2 n 


cbservaticns o X. Furthermore we can compute the 


Шаглапсе СЕ I as: 


Х 


n 1 
Var[2j Маи | 
т=1 п 


1 n 
--- Уаг| 22 Х | 
2 1= 1 D 


n 
n с“ 
м Х] = === « 
2 n 
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Thus we see it is possible to estimate the Var[ X], even 
: а 2 | 
if the population Variance (©) is unkncwn, by simply 
< 
using үсе sample variance 5 as follcws: Уаг[Х] = S sn, 


[4! 
where S = 2 (X -X) / (n-1) (an unbiased estimator of 5 
: 1 


< 
Therefore 5 /п is an unkiased estimatcr c£ Var[X] also. 


Furthermcre, especially for a Normal population, we can 
ertain a ccnfidence interval for the mean using the 
t-statistic. 


Eut, alttcugh fcr X there exists a direct assessment of 
variability, for estimates of other population parameters 
such as Skewness, Kurtcsis, Coefficient of Variation, ard so 
Merten, this is not so simply cbtained. Thus several methods 
have been introduced to oktain assessments of variatility 
and for two cf them (Sectioning cf data and the Jacknife) 


the Subroutines 'SECTN' and "АСК! will be presented. 


Б SUDHCUTINE  SECIN 


I Tescription 


к. анир = c-—- m сто 


'SECIN' Subroutine is used £or assessing Variatility 


of estimates from data based on the Sectioning Method. 


іне гасіс idea of this rethod is: 


Assume we have n independent ortservations zt y 
Е, X from Зоо рида оп оп ао ссн d$ drstributron 
22101 вх). Let O be a parameter cz E (x) апа 
O(n) = ШЕ i) be a statistic which estimates Ə , 
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ла 


Now we want to estimate the variance of Yn) Which is a 


new Randcm Variable and we work as follows: 


ПИ Divide the sample into г disjoint sections. of 
size m (r Should be such that n=mr. If this is not fossible 
scme of the latter data is discarded). № 

EE Есг each section fcrm the same estimate Ә; (а), 
ЕЕЕ ІІ, 2,... ,г 


Bee Ссври е the average cf 9.0m), that iS: 
- r 22 
Día) = --- 29, (а) , 


dcing so, we have an estimate of O and 


= ME S 
g(9()] » E(--- E? (0) 
г 1=1 1 
1 жә 
epum 
it 
- E[ 9 (m) ] , 
alsc we rave: 
25 1 m ^ 
Var[9(n)] » Var(--- > 9 (а) ] 
r 1=1 1 


D fe 
= --- Var[ Z 9 (m); 
1-1 1 


2 

r 

| E 
ш --- Таг( Ә(ш)) 

т 

1 1 Е Б 2 
E es X (Q (m -9(n)) 

Lore] yi 1 

1 2 
= сыы 

Т 9 (m) 
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Thus we see that Varf E and Var[ 9 (п) т л” ре 
unbiasedly estimated from the sample variance of О (а), tne 
main advantage and the purpose of the secticning of data 
nethod. 


Eut the main disadvantage of this method is that we 


would like r to be as large as а order tc make 


the variability of the variance estimate 5 as small as 
possible. This however пау be worse for the Hias 
Meepetties cf the estimation procedure. Therefore the 


ЖИЕС cf r is a factor that should be considered. 


It should be noted here that if (m) are 
approximately ncraal variates (This can te tested Ey the 
previous described Subroutine ‘'NORMPL') then confidence 
intervals for the unkncwn parameter can be obtained based 


Gm the t~Statistic, in the following way: 


22 (1-4/2), (г-1) 
L 


where: 

ЕЕ 15 the fean of the sectioned data statistics 
crtained from the column named 'mean' cf the second table 
(‘estimated ¡parameters of the sample parameters!) cf tne 
Program output. 

1782 | КИН 
S = VE 15 the standard deviation of the 
9 (n) 
Sectioned data statistics divided by the square root of the 
Humber о: sections, obtained from the last column of the 
same table of the program output. 


t the (1- 2/2) quantile cf the 
щеше 2) (ret) 


t-distributicn with r-1 degrees cf freedon. 
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!SECTN!, using the Secticning method, estimates the 
ег!їо4104 statistics: Mean, Median, Variance, Standard 
Leviaticn, Ccefficient of Variation, Skewness, Kurtosis, 
Minimum, and Maximum. (The formulas which have been used tc 
comfpute these estiunates are described in the comments cf tne 
Subroutine.) The first table is then printed by the 
program, containing the values of these parameters. Then 
(SECTN*, using the computed estimates of all sectiors for 
each parameter, estimates the Mean, Median, Variance, 
Skewness, ktrtosis, and Standard Deviaticn divided ty the 
square roct cf the number of sections. The second table is 


printed ccntaining these values. 


here are three restrictions in using 'SECTN':; 

ШО The nunber of secticns must пе пс greater thar 100. 
If it is, a diagnostic messace is printed and only estimates 
from unsecticned data will be given. 

2. The number of data values must be greater than 3; 
otherwise a diagnostic message is printed withcut any 
calculaticn. 

3. The size of each section must be greater than 3. If 
it is less than, or equal tc 3, then the program gives 
estimates fcr tne entire set of (unsectioned) data. A 


agnostic is printed. 


complete description of hcw 'SECIN' opsrates is 
given in the Subroutine. Furthermore, a Summary is given by 
ping cn the terminal the command DESCRIBE SECTN under the 
CMS envircnment. “he following response will be given on 


the terminal when the user types the above ccmmand: 
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SUDBCUITNEC SECIN 


Пеко intended to estimate а set of kasic 
statistics cf a given set СЕ observaticns using the 
'Sectioning СЕ data nmethod'. Also for each estimated 
statistic, estimates of some basic statistics such as the 


mean, standard Deviation, and so forth, are given. 


It is called py: 
CALL SE CIN ху NS.) 


Where: 
X Ше сие агган ок data 
N Is the number cf data values (must be greater 
than 3) 
K Is the number of desired sections (no greater 
than 100) 


Note: k should be a number which minimizes the number of 
data points that will have to be discarded. 'SECTN' rlaces 
the data intc the equal size sections discarding any data 


left over. 

Fer k$3 or k>1C0 or (п/к) =3 only estimates from 
unsectioned data will be given and no estimates fcr tke 
estimated statistics will be computed. 


Неш 15 ехрескед 15 п53 . 


Mcre informaticn is given in the sukrcutine. 





3. Using 'SECIN' with Isiephone Dat 


== dms —m ош ор 


Set 1 


'ЗЕСТМ! was used on Telephone Data Set 1 to assess 
the variability in the Mean, Median, Variance, Standard 


Ес, Ccefficient of Variation, Skewness and Kurtcsis. 


The 672  data-points cf Telephone Data Set 1 were 
broken dcwn into 16 disjoint sections with 42 data-foints 
per section. Because cf this break-down nc data-pcints were 
discarded. Ccmpare the values cf row "!unsectioned! (see 
Eure 21) with the values of the correspcnding statistics 
раседа by the Subroutine ‘HISTGS/HISTFS' (sea figure 1). 


The values are the same. 


Now if we want to assess tne variatility in ary of 
the above parameters we proceed as follows: 

1. We take the mean (9 ) of the parameter, whose the 
variability we want to assess, frcm the tabie under 
ШТАТЕ EARAMETEBS OF SAMPLE  PABRAMETEBS'. 

2. We таке, from the same table, the value (S) frcnm the 
last column (STD. DEV./NS**.5) for the same parameter. 


IO Sing the t-statistic with k-1 d.f. and the formula: 


У + 55 
(21522721 , (x-1) 


we get a (1-4)% confidence interval for the рагашетег. 


Example 
їп сгдег to find a 952 confidence interval fcr the 
Skewness cf Telephone Data 1 we take: 


О = 4.9343, 211 an 20-18. 
„975,15 


те гоге the 95X confidence interval is: 


4.9343 : 10312%2,131 = 4.9343 + .64595 
===> (4.2884, 5.5802]. 
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It shculd be ncted here that the use of the Variance 
estimate  frcm the sectioned data in order tc get confidence 
intervals of the parameters is based cn the normality 
assumption and the independence or the estimates frcm the 
Sections. The normality will depend on the number of 
data-points in each secticn, which should be kept large. 
This reguirezent, however, conflicts with the need to make 
mages number of sections large to reduce the variability in 
the estimate of the variance of the statistics. The 
skewness estimates of each of the 16 secticns, can te run 
through the normal plotting routine to see whether the use 


of the t-statistic confidence interval is valid. 
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C. SUBROUTINE JACK 


This routine is used for assessing variability of 
estimates of data based on the Jacknife method, ané for 
Meaucing the bias in estimates also. It is particularly 


useful fcr data with small sample size. 


A big picture cf the method has as fcllows: 


Let їг 22 пи X be a sample of n independent, 
n 


identicaily distributed observations from a population with 
unn distribution function Р (х). Alsc let Ә ре ап 
x 
unknown parameter or F (x) tc be estimated. Furthermore 
X 


suppose a method (biased ог unbiased) те ата:гарје “ог 


estimating 2 Then we proceed as follows: 


le Divide the sample size into r disjoint groups, each of 


size m (r should be such that n = mr. If this is not 
possible scme of the last data-points wili Ее discarded.) 

2. Compute the estimator Ә (а) ОБ 9 based on all n = mr 
cbservaticns. = 

EE Compute the estimator Ө (а-а) based on the n-a 


th 
cEservations, having deleted the 1 guo 


^ 
4. Compute the so called 'PSEUDO VALUES' У (г): 
ЕЈ 


9 (n) 


d: 


27. = (r-1) 9 (n-a) 2 Тс) Do SHE а 
nd 
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Then the Jacknife estimator is defined to ke the average of 


Eseudo Values: 


1 


< < = 


Мн 


^ 
9 (n) 


m шаа 


9 (n) 


1 225 108: 
XM SI O(n) s = 9 (2-1) ) 
E 121 i 


ла r-1 КЧ 
E O(n) = puc Ex Э. (а-а) š 


Havirg so defined the Jacknife estimator it can be 


proved that this eStimator is an unbiased estimator cf 27, 


except fcr terms a „and higher order, assuming that tne 


kias has tne form cf Ef 9 (n) ] = © + а/ц + terms of Higher 
E ЕШ n. That is, the Jacknize estimatcr eliminates a 
n bias term. Namely i£ 


~ = -2 
E 9()1]-9*an «оф ) , 


then we have: 


r-1 


Ма) ] s E(r9(n) - E 


O (n-m) jJ 
18: 


и Мун 


а 


© ut pone 
= rE{ O(n) ] = бес 13519, (в-а) 1 


= E r-1 ra Е - 
= r(0+ an * O(ín ) ) EM шалга 100 


га а 
эг гэг. 227) 

n nir-i! 

ra poda 
О o CN 0... 

n nir) 
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г9- rd + ---- - -5- + 9 +... 


ШІ m 


що. 


A special case of Jacknife estimator can be 
Mentioned here, called the ‘Complete Jacknife Estimatcr', 


where r-n and ms1. 


оспе properties of the Jacknife estimator follcu: 


-2 
ПИЛЕ is unbiased up to order n . 


2. The ''2seudo values! can be used to obtain variance 
estimates for tne Jacknife estimator since they can be 
considered (see (18]) as approximately independent and 


Г А ^ 
identically distributed. Thus  (r(r-1)) = (9, (2)-9(1)) 
i= A 
^ 
Should be an approximate estimate of Var[9 (r)], 


and 
722) 


Should be approximatelly distributed as a t-statistic with 
ESI d.f. This procedure is particularly useful if the 
number n cf data-points is small, but it must be used with 
care. 

U Por large n it can be shown that ісі the ccuplete 


Jacknife estizatcr, under very general conditions we have: 


IN ^ 
Vari 9 (0) 1 ---> Mag 9 Са s 
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Р? Program Structure 


Subrcutine  'JACK' is a FORTRAN-callable Subrcutine 
E tekes a set of data, groups it into r disjoint  grcups 


ЕП 
and for sach set of the r-i group computes and prints in a 


table the statistics: Mean, Median, Variance,  Stardard 
КОЛО Оп, Coefficient of Variation, Skewness and Kurtosis. 
The same statistics are also computed for the ungroupeé set 
of data. Then, using these statistics, the  Jacknife 
estimatcr of the above parameters are computed and printed 
in another table alcng with its Variance and Standard 


Beviation. 


The program is divided into the mair program ar tne 
PeeeOutine JACKES. The main Ercgram groups the data and 
successively calls the Subroutine !ЈАСКЕ5" in order to get 
the estimates of the arove parameters for ungrouped data and 
fcr each group as well. A table is then printed containing 
the estimated parameters. The  'Pseudo-Values' for each 
parameter are also computed by the main program, ard are 
used to compute the Jacknife estimator for заса parameter 


and its Variance and Standard Deviation. 


memwavyoid divisicn by zero the number cf observaticns 
must be greater than three, ctherwise no cutput is civen. 
For the same reason the expression (1-1) (n/I) must evaluate 
tc greater than three, otherwise the program will give only 


estimates cf ungrouped data. 


A complete description of how 'JACK' operates is 
given in the  Subroutine. A Summary 15 printed on the 
terminal Бу typing DESCRIBE JACK under the CMS environment, 


which responds: 


ЗӨ 





SUBRCUTINE JACK 


Jack is intended to estimate the statistics : Mean, 
Median, Variance, Standard Devwratione —COeCEfLICrent of 
Variation, Skewness and Kurtosis of a given set of 
independent  cbservaticns, using the Jacknife method. Іп 
addition to the above parameters which are given fcr each 
group, the Jacknife estimator along with its variance and 


standard deviation is given fcr each parameter. 


ШЕ іс called by: 


CALL ШЕП ( X, Хо, STAT, N, 18) 


where: 

1s the array of data of dimensicn N 

X S is a work array of dimension N (returns crdered 
data) 

STAT is a two-dimensional work array, dimensioned by 
(пе 7) 
is the nugzber of data-values 

IG is the number of groups pius one (-2r*1j 

Ncte : r must be such a number as tc minimize the 


number of data-pcints that will have to Бе discarded. JACK 
places the data into the equal size groups discarding any 
data left суег (last observations in  crder of original 


шшви+). 
Nc cutput is expected if N $3 3. Furthermcre if 
IG < 2 or (IC-2) (N/(IG-1)) < 3 only estimates for  ungrcuped 


data will be given. 


Моше поета васп is given in the SuLrcutine. 
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Tc assess the variability in the Mean, Median, 
Variance, Standard Leviation, (ен ис веб Of Varlatron, 
Skewness and Kurtosis cf Telepfone Data 1 the 54 гсиїїпе 
ШООК" JS now used. In order te avoid to discarding any of 
the data-points we use 16 groups with 42 data-points рег 


group. 


Observing figure 22 we see that the values of the 
parameters for ungrouped data are exactly the same as the 
corresponding values of the same parameters of figure 1 
waich was produced Бу using the Subrcutine HISTGS/HISTFS. 
Now in order to assess the variability in any of the akcve 


parameters we have to use the t-statistic with (r-1) 


mee alcng with the values printed under 
'ESIIMATED JACKNIFE PARAMETERS' using: 
O+ st 


acm Irc) 
«here 9 is the Jacknife estimator апа S is its Standard 
Leviation. Fcr example to assess the variability of the 
Skewness with a confidence level “= .05 we have: 


387 7.37321, S = .900118, t zm E 
10751165 


Штеге=гоге а $5% confidence interval of the Skewness is: 
7.37321 + .900118*2.131 = 7.37321 + 1.91815 
===> (5.4551, 9.2914] . 


This shculd be compared with the point estimate 4.9343 and 
Ae confidence interval estimate (4.2984, 5.5802] obtained 
from sukroutine SECTN. The data is so skewed that one will 


inevitably have trcubkle here with the Jacknife procedure. 
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A. THE EARMA(P,Q) ERCCESS 


The starting point for these processes is the defirition 
cf a first-crder autcregressive model fcr a stationary 


sequence cf random variables {X }: 
i 
EE = rx DU ЕО ене 


ШЕ Пе marginal aistritution of the X is fixed to ре 


exponential wich parameter A for all i; 


БЕО < х) = 1-е 2 ОТО ЕСЕ 
ДЕ 


then € is zero with prcbability r and  exrcnential(5) with 
í 


шазшара 11%у t-r. Thus 


IX. м.р. г 
au 
M (1) 
1 
EX TRE н.р. Ier 
к= т 
where (E ) is a sequence of 1.1.4. exponential (2) random 
1 . 
variables. The process defined by (1) is the exponential 
autoregressive process of order 1, the EAR(1) process. The 
: | . J 
Berrelaticn structure is r(j) = г, ј = 0,1,2,... . 
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De Order moving average  expcorential  rccess, 
EMA(1) is defined as 


ФЕ. м.р. 6 
х = j (0 < 6 


Г 5 12: Jo ОЕ 2 ә ) 
1 
ФЕ + E ÓN 
de = 1 
and has ccrrelations (1) = (1-2) and Ра) О. 
The general EMA(q) model takes the form 
ДЕ “.р b 
q 1 q+ 1 
E + Ë м.р b 
ШЕГІ тле] Ч 
EE 7.7377 /  .... (2) 
1 
(Е A, E Св Wep. b 
gu о Tigri 2 
/ + B Tc E, p * E WOOD B 
g i q-1 i-1 1 i-q+1 ша 1 


H. 
O 

H 
© 
їл 
a 
ч 
с 
~ 
~ 

N 

IA 
mah 

l 

И 


2:27 a CS 


1 2 q 
SE e ende 
4 
Di = 1202 ми. 212 ae, 
1 q i su 
(1-2) ... (1-4 | 17! 
Ч 1 


The serial ccrrelations are given by the equation 


а | 
z= ob. (2250 s 4) 
v= J У У%7 
(3) |. _ 
? m) = Ы ) = 
0 (qt? < j <@) 
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ШО сспуегт this EMA (q) process into a mixed 


autoregressive  moving-average exponential prccess of crders 


1 and q, called ЕМА(1,4) we replace E. in (2) with 
1-9 
| ERA : м.р. г 
al - 
А. - 1 
1-9 
ГА жЕ. W р» (1220) 
15421 12521 
Bor further results con these processes and their 


properties see Gaver and Lewis (1978), Lawrance and Lewis 
(1977), Jaccbs and lewis (1977), and Lawrance and Lewis 
(1978). 


E. USE CF AUTOREGRESSIVE (EAR(P)), MOVING AVERAGE (EFA(C)) 
AND MIXED STEUCTURES (EARMA(P,Q)) IN MODELLING QUEUES 


Ecusider for simplicity а queue with a single input 


stream and a single server, and а first-ccae-first-served 


(FIFO) ЕР Се “discipline. Let S , 1 = 0,1,2,... , dencte 
1 
| | : ‚th | 
the service time for the 1 arrival, and let х, 
1 
| | MS 
1 = 1,2,-.. , denote times between arrival of the i and 


| st | 23 
2-1) custcmers. Ав is usual we assume that the first 
customer (with service time 2 arrives at time zerc and 


finds the queue empty. 


If the (S) апа {Х } sequences are i.i.d. exponential 
ЩЕ д 
random variatles with parameters ^ and Q respectively, we 


have the M/My1 queue. 
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Now let 


Е be exponential(2) and independent, i = 0,+1,... 
E 


ОСЕ Е. . 


6, Бе exponential (a) and independent, i 
i 


We want to model queues with correlated (autocorrelated 
and/or cross-correlated) service and inter-arrival times, 
the service and inter-arrival times both having marginally 
exponential distributions. We also want the queing model to 


include the M/M/1 queue as а Special case. 


There are five simple possibilities which we pur fcrward 
here, based on the use of EARMA (p,q) processes, givinc what 
we call an EARMA 34D/MD/I queue. 


Zn what follcws let р and p ke, respectively, the 
5 X 


crder of the autoregressive components of the service time 





sequence {S } and the inter-arrival time sequence (X ), and 
i al 


let 4 ,4 ke, respectively, the crder of the moving average 
5 4 


component of {5 } ара (р. These ¡parameters can take 
x зі 

MENTES 0C,l,... . ЕР. = "брака = 05 the sequence is 

independent; if p = 0, the process is purely moving average 

and if q = 0, then the process is purely autcregresive. 


The five possibilities or cases are as fcliows: 
ЭС {S } be EARMA(E ,g ) Over (š , 2 ; s... 
1 = ад 1 1-1 
ШЕС 1} = Е). 
st 1 
Thus the arrivals are a Poisson process and the service 


times are autocorrelated, i.e. a mixed autoregressive 


moving average expcnential sequence. | 


106 





3) 


4) 


РЕ ЕЕ) 1 = (0,1,2, ... 3 
2 al 


ШОО Бе БАНМА(р „а ) стег (6, <= , ...). 
T X < ot 1-1 


Then the service times are independent and the arrival 
process is non-Poisson because of the dependency between 
the inter-arrival times, i.e. the arrival process is a 


point prccess with EARMA(p,q) structures. 


mc (5 } Ге БАНМА (Е „9 ) over (Е, БЕ were 
i 5 8 ah 1-1 


Mec [X ) te EARMA(E ,g ) over {&, 6 oo Wa eis) 
1 Z X 1 | 


The service and arrival processes are autocorrelated, 
but the two processes are independent. The marginal 


EE triputions of {5 | and {X } are still exponential. 
p 1 


To ccuple the two processes, with resultant  dererdence 
in the arrival and service processes, the simplest 


procedure seems to be the following. 


Let {S ) Le EARMA(E ,q ) in the following sense: 
i SEES 


iSc) 15 EARMA(p ,q ) over (Е, ЗЕ, ос v xn NE 
E E 1 E A 


Then if Х = =, таа Од ие ауел тра (5 is en 
T 1 1 
autcccrrelated sequence and also  cross-ccrrelated with 
ши - съ, i.e. О ТЕКСИ ОУ те ce devendsme 
n E аи 


sequence cf random variables with  expcnential marginal 
МЕНЕ ions although {X } is still Poisson. Note 
i 


that the S sequence is an autccorrelated sequence 
1 


because Or the cross-coupling, but not a pure 
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Let 


s (p q.) sequence. It does have an exponential 


marginal distribution however. 


Moran example let p = 1, а - 1. 
5 S 
EN Е 
0 0 
S = 6 ж мер. ВИА 
(|7 ^s >. 4, 1222 
- Жі wep. (1-4) 
S. 1 1 
S = Ë W.D. А = ГА W Е 
2 2 2 У 5. 2 1 Е 
= ДЕ +A М 2270 = ГА +25 я 1-г 
6, 2 E 5 1 A 2 P 
с, = ba яз р Á А = ГА чор» TE 
1 ST S т Le 
- БОНА wep. (1- ) = rÀ TEE ч.р. (1-1) 
Si d S 1-1 ді 
In these equaticns one can replace & with X and then 
1 1 


the cross-coupling between the sequences (S ) and {X } 
1 1 


becomes apparent. 





Interpretation. We have positive correlation between S, 
i 


and q previous inter-arrival times (when ра - 0). If 
S 
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4!) 


5) 


' 
Ше © s (j = i,i-t,..., ee are short, then S 
j i 


will Ее short if all the a previous inter-arrival times 


are short. This mo ere the server tends 


dels the case sh 
b 9 


F+ mo 

to srsed up if the queue gets long. Cf cource he also 
slows dcwn when the queue gets shcrt and it is not 
immediately clear what the effect on an average waiting 


time will te. 


Similar to 4 but (X ) is EARMA(p ,q ) over 
at XZ X 


=. ч, АР шоог: + Е. 
i СА 1 a 1-1 1 1 


Interpretation is nct clear, but it would have the 
same effect as balking in the input stream. If service 
times get long (and presumably tne gueue gets long) then 


inter-arrival times get long. 


А more general possibility which allcws one to mcdel 
dependency in the input stream, in the service stream 


and dependency which couples them is as fcliovws: 


Let (0) be a unit expcnential independent sequerce, a 
1. 


driving sequence. 


Met (2 |) be BARMA(E ,q | over (Е, В „...). 
1 S07 1 are) 


mete ik } ke EARMA(E „ЧА ) счег (65, e эг 
т X X Т ца 


ШО {5 } Бе EARMACp ,q 
2) С 


OV eG ТОО 1 
Nc 1 1 1=' 


Met (Х | се БАНМАСЕ ,q ) over (X , 9D,, OD, ,...). 
L C, Ca 1 m TER 
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This is a very general scheme which reduces tc all of 


the previous schemes as special cases.  Fcr example if 


Ё = q = = q = 0, then we have case =. I£ 
Ci C, С, Ca 
q 21, E = (9% 2 = 07 p = q = 0: Е m (and p 
с, с, С; S S X C, 
and q are not equal to zero) we have case 4. 
С 
1 


In all cases the basic equation for M/M/1 queue waiting 
time W still hoids: 


я 10.0 ) = 220217 


These equations can be used to generate successive W 5 
n 


in a simulation. We have not tcuched on hàcw the correlated 
Sequences are started. This is somewhat arkitrary as it is 
not known (except for the EMA(1) and MA(1) processes) пси to 
start then so as tc produce a staticnary expcnential 


sequence. The problem is aggravated, as in cases 4, 4 and 


5, when we are dealing with Livariate exponertial seguences. 
Then the  zarginal rrccesses may be stationary but not the 
BEvgriate prccess. 
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ШІ РЕОСБАМ STRUCTURE 


> General 


Ic simulate the EARMA(p,q) model, a FORTRAN Prcgram 
has been written tc take care of all the existing cases of 
that model and also of the M/M/1 queue. Tepending cn the 
input values which the parameters (as described below) will 
take, the program can be used for creating in each run a 
particular case of the EARMA (p,q) queue. Eecause of that 
generality cf the program it is not suggested that it be 
used for thcse EARMA cases where either the arrival tines or 
the Service times cr both are not autccorrelated ог 
cross-correlated. The reason for this sugg=stior is 
efficiency. It is suggested for use with the EARM2(p,y) 
cases where both arrival and service times are 
autocorrelated and, furthermore, cross-ccrrelated with any 
crder. The program is divided into the main programs and 
ео the Subroutines BETAS, AUTOR and EARMA. 


2. Main Program 


ШОШ fain. program performs ail i/C operations and 
calculations required TOL the desired statistics. 
Specifically it is designed: 

а) To read 
1. All randcm number generator seeds required to 
generate the sequences of the Exponential and Uniform 
Vaciates to be used during the executicn. 
2. The values for the variables: 


N : Number of arrivals to be generated. 





Пена сет Of replicaticns. 

РЕ. 2 Arrival rate, Service rate. 

Е, СВО : Taking the values 1 cr O (1 fcr the 

coupled process). 

11. 11, KSl, KS: The order of moving average for 
arrival and service processes. 

ШОТ, СХ, 051, 05: Taking the values 1 or 0 (ТОГ 

autoregressive case). 5 

Arrays with the values of Ês for 

mcving average part. 

RHOX1, RHOX, BHCSÍ, RHCS : Parameter values for the 


tri 
>< 
` 
aj 
>< 
== 
~ 
to 
п 
` 
UJ 
(л 
=» 


autoregressive parts. 


ШЕ 10 calculate and print 


SUMX (1,3); SUMXM (1,3); 
SUNS (I; J); SUMSMULIJ): 
W(I,J); AM (I,J); 

То) ИЙЕ (1,4): 

(1,4): DM(I,J); 

БОСНИ = 1,2,3,4, о 
where: 


(In what follows k gets the values N/4, N/2, 3N,h, N 
Ме = 1,2,3,4 respectively and J-1,2,...,9M) 
k (J) (J) th | 
SUBA (E, J) = Z X j X the r arrival for 
r=1 r E 
th . - • 
сће Ј replication in tae EARMA queue; 
k (J) (J) : th е 
50 ХМ(1,4) = =, (Z š) ; (XH) tae E arrival 
r E 


en 
for tne J перса стони таии шене, 


K (J) (J) th : 
SUMS (I,J) = = 5 у = tne E Service for 
r-1 r 


th 
the J replication in the EARMA queue; 
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К (J) (J) th 
У ИУ о) = E) , (SM) the r Service 
- [a ДЕ 


en 
rormwthe 4 replication in M/M/1 Cucue; 


(J) (J) (J) (J) 
Я(1,4) = Y = max W + S m 0 
1 k НЕ k-1 1 
мэ | th | | eh 
- waiting time of k атна ав) replicaticn 
of the correlated queue; 


(J) 
ИМ (1,Ј) = Qu The Same as NN but 


Есі M/M/1 queue; 


k (J) 

WB(I,J) = 2, Ч /k 
r=1 r 

K (J) 
aME (I,J) = — (W1) IR 

=] n 
PM = iSUMS(I,J) — SUMX(I,J)P^/K - 1/RS + 1/8Х 
ПИ ИВЛ) = [SUMSM(I,J) - SUMXM(I,J)/K - 1/Е5 41/8Х 


К (3) 
SUMSM(I,J) = 2. (50) 5 
r=1 ЈЕ 


Bere: Between the group-variables (KX1, KX, KS1, KS), (CXl, 
051, QS), (CPX, CPS, CEG), (RHOX1, RHCX, RHOS1, ЕНС5) 
Но (211, ЕХ, 851, BS) there are the following relations: 

1. The number of elements Of each variarle (array) from 
B-group shouid be the same as the value cr the correspcnding 
variable froca K-group. 

ИЛ The 'O' value cf a variable fron Q-grcug implies that 
the correspcnding variable from R-group is nct needed. 

meine value '1'! of CPG dominates any value ОБ СЕХ, ог 
ers. 

Any valid combination in the values of the  arove 
Bazjables specially of the Q-group, C-grcur and К-агсир (0 
cr not 0), generates a particular EARMA(p,g) case and is 


accepted ty the program. 
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Ни? G 0 U prroaCch in the main program is tc 
generate at once and for each replication three independent 
exponential sequences with parameters XL (arrival rate), SL 
(service rate) and 1 (unit exponential) which are stcred in 
the arrays EXPXL, EXPSI and EXP1 respectively. Ther the 
Subroutine AUTOR which is involved with the autoregressive 
part of the model is called (if it is reeded for the 
particular EARMA case). А lcop follows which is executed as 
many times as the value of N. From the 1сор the EARMA 
Subroutine is called which is involved with the acving 
average part cf the model. The number of calls of Subrceutine 
EARMA depends on the particular case which is being rur. The 
statistics and all infcrmation that will help us to analyze 
the model аге alsc computed inside the lccp and stored in 
tne arrays described aktove. The program continues execution 
until the desired statistics for all М replications have 
been caiculated and gathered. Then tne cutput part of the 
Eram fcllcws which gives us the values of calculated 


statistics on paper and on punched cards fcr analysis. 
3. Subreutine BETAS 


This is a simple Subrcutine and its ¡purpose is te 
return an array 50 8Х1/Х/51/5 with elenent-values as 


иом: 


T 
SUMBX (1) = = EXE]: TA value of K-grcup 
L- 


De used for the choice of the order іп the moving 
average. It is called at most 4 times during the execution 
and only if the K-group variables have value greater than 
B'. 





4.  Subrcutine AUTCR 


'AUTCR' accepts in each call an excnential secuence 
(generated previously) of variates and transfcrms it into an 
autocorrelated sequence. It is called cnly if any of the 
Q-group variables equals 1. For each Q-grour variable having 
value 1 it is called M times. 


D. Subrcutine EARMA 


Шин TS tne main Subrcutine of the program and it 
ши рееп created so that it can be used fcr all EARMA(E, 9) 
Cases. Because of the generality quite a few parameters are 
transferred and in each call just one value is returned. 
[its Value is the DU arrival or the Бэл Service time as it 
has been acdified because of the roving average. It is 
called only if the particular EARMA(p,g) case requires 
Moving average; that is only if any of the K-group variatles 


is greater than 0. 


Scme CPU-times have been gathered in the course of 
the simulation of some EARMA(p,q) cases, using tte IBM 
360/67, FCRIRAN H compiler as follows: 


N M K-GROUP 051 C-GROUP TIME 
2000 500 0 1 0 ВН. ЕС. 
500€ 500 0 1 0 DA а сокове 
10000 500 0 1 0 ISSO. SES. 
1000C 500 0 1 CPS=1 5S9 Vin.: 0 Sec. 





From these results we see that the value cf N is the main 
metor that affects the CPU-time (linearly). We see also 
that one cf the C-group variables may increase the CEU-time 
ЭГ. Че ignore the factor of M (it is another main 
factor as N is) because че may step thé program at any 
replicaticn, getting the last random number generator seeds 


and tnen continuing ancther time using these seeds. 


D. BATISIS OF SIMULATION FESULTS 


Since ne analytical properties of the EARMA (o, q) ucdels 
can be derived, a simulation has been done to study their 
properties. But as we have noticed, this program is a 
general program accomodating all EARMA(p,g) cases. Since 
our study is limited to only the two cases below, in crder 
to save CPU-tine a modified program has been usec that 


Simulates these particular two cases. 


The mrcdel: 


(X) = €, 
1 1 
(57! is EARMA(1,0) over ааг Үл 06-13 
i 1 i-1 
must is: 
ГА м.р. г 
i-1 
5 = where А = E 
1. Q 
im + E шар. (1-г) 
1-1 1 


Tc simulate this model we chose the values .25, .50, 
C 


К 99 for traffic intensity t and tne values .2-7, 0% 


116 





МО. 95, „96 for тһе correlation r, in order to cover a 
representative range. тыу ыас follows the String Smetter 
Emus fcr the run with t-tt and г=гг. That is the run 
ВО stands for the run with t=.25 and r=.50 For each t 
апа Гг we ran the program and the required statistics were 
gathered, analyzed and plotted using the Subroutines 
ENIGS/FS, NCRMPL and EXPLT. The sample statistics from 
HISIGS/FS were also taLlulated separately frcs the figures so 
that one could cbtain an overall picture of whether or not 
the waiting times had converged in mean and in distritution 


Eu ue liziting distrirution. 


From the analysis of all runs and plots we have the 
Sing results: 


1. Simulating +һе туре 5#25#гг апа S#5C#rr model (with 
ши 25, .50 апа .90) it was possible fcr the W and W to 
reach the steady state for N=2000. On the cther han for 
БЕ type S4254rr and Sé5Osrr (with rr = .95, .98) it was 
necessary tc go up to N=10000, in order for the W and W to 
reach the Steady state. Thus we see that the high 
Memmelaticn affects the choice cf the value cf N, requiring 
N to increase as the correlation is increasing. But not 
cnly the correlaticn affects the choice of N. Therrrarftic 
intensity affects it much more, since for the type S4#95#rr 
and S#99#rr the W and W do nct reach the steady state for 
N=10000 even if r is low. Because cf the C2U-time 
requirements we restricted ourselves tc a detailed study of 


the particular case 5%99%98 increasing successively the 
value of N ug to 320,000, where the W and Ñ appear to Бе 


close to the steady state. Thus we may ccnclude that high 
Matric intersity and/cr high ccrrelation require a large 
Value of N in order to achieve steady state. Figures 23a 
E cugh 23d, 24a-24c, 25 and z6a-26d justify the abcve 


 псіпсіспв. In figures 23a-23d, where we are dealing with 





| th 
the waiting time of the 2000 ШТ arrival for the case 


t=0.50, г=0.25, we can see that D has convergec to a 


value of about 6.3, and furthermore a straigtt line appears 
in the rer under  EXPLT. We also see approximate 
ccnvergence for the correlated queue for the W „обе 15 
presented cn figures 28a-28c, and its vade has converged 
ая та1џе cf about 3.03 (figure 23d). Note that the 
mEuroduction of correlation into the service time has 
increased the average waiting tire from 2.5 tc about 3.03, 


cr by approximately 204. 


ќе cannot say the same fcr figures 24a-24c «kere, 
because cf the high ccrrelaticn (г-0.98) we do not have 


convergence of W and W, although N=10,000, and for figure 


25 where, because cf the high traffic intensity (t=0.€9) we 
dc not have ccnvergence of W and VW. Figures 26a-26d also 
show us the поп-сопуегдепсе of W, W, because of the high 
correlaticn and the баби века вас intensity (t=0.93, 
r=0.98). This is the case for which we eventually tad to 
Carry the simuiaticn out to about N=500,000 to ooserve a 


stationary qusue. 
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2. As гог the steady state of the value cf E(W] anóà E(W] 
асе М was required for high t and/or r, also for the 


ВЕСУ stete of the distibution cf W and W a large value 


Mes tTequireá. informally testing their distrikuticns, 
Myeusing the Flotting Subroutines and the values of their 


parameters skéwness and kurtosis, we may conclude that: 


fee ite distribution of W, given that W>0, has an 
exponential form, if the plot results in a straight line 
Midget EXPLT and the / and {, Farameters have the values of 
about 2 апа 6 respectively (see figures 23a-23d). This 15 
true only if the steady state has been reached while it is 
Me evident what their distribution is if the steady state 
has not been reached. Thus cbserving the figures 2Ua-24Ud, 
where the on has not reached the steady state yet, we 


ВО Еазп ncthing about che distributicn of 3. 


Note that for the M/M/1 queue the steady state 
(би сл, о W, given that 420, 15 eaponential. The 
distribution is not kncwn for the correlated gueue and this 


Was cne of the objectives of this study. 
Ory. The distribution of W has the ncrimal form if we 


get a straight line under NORMPL and if the values of 
skewness anc kurtosis are arcund '0'. The поста БЕН 
аззишр а сп hceläs for the steady state cıly, while the 
ncn-steady state distribution looks like an exponenticnal. 


Figures 2127а-27с give us the feeling that ron tas an 


exponential fcrm, since a straight line appears under EXFLT 


(see figure 27b) and the parameters 7), jj have the values 
2.3 and 7.3 respectively. These are close to the actual 


values 2 and 6 of the exponential distribution. We can see 
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mom the picts that as the Y gets close to the steady 


state, it leaves the exponentional form and ccmes closer and 
Moser tc the normal fcrm. Figures 28a-28c, 234 shcw us 


gua t IE оо has already reached the steady state and 


furthermcre its distrirution has eventuaily taken the ncrmal 
ШОШ since a straight line appears under NORMPL (figure 
28c) and the skewness and kurtosis have the value .7 ага .8 


respectively. Compare its ОБ ripu im with the 
Nu c:rzibution of NOUO of figure 27b where W is not in the 
steady state. Again we may state here that convergerce to 


the steady state of the distribution of Y and W requires 


much larger N than convergence to the steady state cf their 
mean values for the same t and r. And alsc the higter t 
and/or r the larger N should be. As an example we can see 


that even thcugh the steady state (in value) of Ñ for 


54252398 type has been reached fcr N=10000, the steady state 
in distrikuticn has nct been reached yet. See figures 
z4e-24g where we can observe that the distribution of LAS 
starts to leave the exponential form and to домен the 


normal rcrn. 
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Ccntinuing the analysis of the results we present 


s che results we got from all runs, concerninc the 


values of E[ WR] for various values of t, r, N. 


N E Е E(W](correlated) E(W](M/M^/1) 
2000 225 225 22720 В 
2000 229 „50 .4326 23:32 
2000 229 . 90 27-73 23522 

10000 „25 295 1.0576 „3332 
10000 225 -98 1.7077 2392 
2000 „50 225 1.2106 „9985 
2000 „50 „520 155813 22985 
2000 250 „90 52.0785 9985 
10000 2:50 295 93196 -9935 
10000 == 2 21.0542* 295995 
10000 9 5 225 23:5651 18253528 
10000 205 -50 ео? 10.3928 
10000 25 „90 123: 9350 ЇЇ 75-28 
10000 42 295 196.2460* 1223923 
10000 295 . 98 342.9240ж 1323528 
1000C „99 225 62:3285% 5593371 
10000 299 250 71.9370» | 
10000 29 290 193.9113* 522337! 
10000 299 25 22119 ЕЕ | 
10000 m 2:98 426.6750* 54233711 


Note that these values have peen scaled from the values 
Seeained in the simulations by multiplying Бу 55. Thus for 
t=0.50 and r=0.25 tne table gives E{Wj(correlated) as 
1.2106; the value in Figure 23d is 3.026x0.4=1.2106 . The 


values marked * have nct converged. 
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A study on these results gives us the following 
cbservaticns: 


ИЕСІ given traffic intensity t the value of E(W] is not 
constant, as jit is fcr M/M/1 queues, but it depends cn the 
value of the correlaticn parameter г. Thus we see that the 
value of E[# ] increase, as the value cf r becomes larger 
and larger. Furthermore from figure 29a, where E[&] is 
plotted versus г мг сап see that the rate cf increment is 
Moremconstant (therefore 1s net a linear function of г) but 
it increases with r and eventually gces up infinitely аз г 


foes to 1. 


2. For given correlation с, we see that the value of EL #] 


increases with t, the sane as happens in the M/M/1 queue 


(see figure 29b). Eut regardless of the value of r, E{ 4) for 
М/М/1 queue is always less than for the EARMA model. Thus 
we may state that the waiting time of the autocorrelated 
service time model is greater than the corresponding waiting 
time of the M/M/1 model. 
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3. The values of У] for the type 58499898 model 
increases with N following the curve in figure 30 a fact 
that is known to happen in the 4/M/1 queue also. We can see 
from that figure that the rate of increment comes dcwr as N 


increases ard hcpefully should come to zerc as W reaches 


the steady state. Ccmparing then the M/M/1 queue with the 
correlated queue we can see that the M/M/À1 queue ccnverges 
faster than the correlated especially when r is high. Thus 
we can see frcm figures 32a-32c that both the value and 


MestLibuticn of the EG for M/M/1 queue have reached the 


Steady state, but the ссггезрспдапа и cí the correlated 
queue (figures 33a-d) has not 20 least in its 
distribution). (See the skewness and kurtcsis values for #1 
EAR, W2 BAR, W3 BAR and W4 BAR in Figure 33d. These appear 
tc be decreasing tc zero, but have certainly not reached 
zero by N=10,000.) Jooking at figur- 116, ме сап cbtain 


the result that neither 00 fer che correlated queue, = ог 


Ф 000 fcr the uncorrelated queue has reached equilitrium. 


mus is because the skewness and kurtcsis values are 


ExUnrztely nct close tc zero. 
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Ancther result obtained from the simulation 
concerns the comparison of convergence cf W and W. Frou the 
programs we have run we observe that W ccnverges slightly 
faster than W (at least for low t and r). Also the standard 
deviation of Wis much less than of W. In figure 24d we can 
see that М has already reached the steady state but not 
и : Ccmpare also the values of their  stendard 
Bon. (.005 versus .2). Figure 34 gives the results of 


the type S#50#98 model for N=10000. We can see that while 
% is already close tc the steady state, the fact that w has 


not reached it yet fcrces us to continue the simulation by 


increasing the value of N. 
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mM model for the crcess-correlated queue is: 


мај СС 
i? A 
i } 15 БАБМА(1,1) Over CE. ELA ОАЕ. 
2 l i i-1 
Mat is: 
Е. ч.р 2 
3i 
5 = 
n 
BE + A м.р (1-6) 
1 1 
where 
20 Wep. I 
1-1 
А. си 
3 


EA. ЫГ 1705 ор. 
ш | 1 


Ada ln we chose the values .25, .50, .95, .99, fort 
ПИ .50, .90, .55, .98 for r and we run a total cf 20 
шиг (0 ccver all ccakinations of t and гс, naming each run 


as SC#ttirr. 


Analyzing the results obtained from the  sizulation 


we ckserve the fcocllowirg: 


Ш- Тһе value ог ре traffic intensity i:e anc the 


 гв1атісп г іс the main factor for convergence of W and Y. 
mat 15, high Ера лас intenci апас high Correl tion 
require N to be large in order for W and W to react the 
steady state. Thus simulating the models with «5.50 ала 
(135 for-N-10000 it was possible for W and W to reack the 
meeady state. In figures 35a-35d, 36a-36c we can see the 
Mst1ticaticn of that conclusion. The values .107€ and 


87957 of m represent the steady state of 


, апа ў 
10000 10000 
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the SC#50#25 model. On the cther hand, regardless of the 
Ш Пе СІ t the modeis with r2.98 did not reach the steady 
state for N=10000 and regardless of the value of r the 
Beeels with t2.95 did nct reach the steady state for 
N=16,000 . Figure 37 gives us the non-convergence of the 
ncdel SC499425. Because of the high traffic intensity and 
in crder to cet the steady state of that model we have to 


continue simulating by increasing the value cf N. 
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ВРШЕНА ОАЕ trafic intensity and correlaticn affects the 
convergence of the distrikuticn of W and W also. 
Furthermore the steady state of their ПЕТЬ нес, 
requires larger N than is required for the steady state of 
their expected values. In examining that ccnclusicn, look 
TOL example at the results cf the simulated model SC#50#98, 
where although W and W have already reached the steady 


state for N=10000, their distribution has not  converged 
yet. The above results have been cbtained from Figures 


Bea-Joc, 294-394 чїеге Ч are presented. 4e 


10000 Би "10000 
can see from these figures (328a-38c) that although the value 
Ef Boo has reached the steady state (see Figure 39d), 
ШЕ distribution cannot be characterized since the plot 
Ber EXEIT (figure 38b) is not a straight line anc the 
parameters (153.34 апа ),=11.0 are far away frcm the 
desired values 2 and 6 . Ihe same results can be obtained 


from the figures 39a-39c which deal with Y : 
10000 
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Bro test the distribution of N and Я we also used the 


plotting sukroutines and the values of their skewness and 


kurtosis parameters as we did with the queue with dependent 
service. 


Ihe results obtained from the analysis are as 
Шо: си5: 

ШІ. Given that W>0, we may say that Я has an 
exponential form, since the flots under EXPLT Subrcutine 
Et ina straight line and the values of У, апа аге 
around their actual values 2 and 6 respectively. For example 
the SC#i0#¢25 model (which is represented Су the figures 
35a-35d) , where the steady state has been reached, gave us a 
straight line and the values [ES and 4.7 
fcr j| and |, respectively. ќе have to state here that the 
expcnential distributicn comes only if W has reacheé the 
Steady state, and we cannot have the same result if the 5 is 


И na transient state. 


(21) Есі the distributicn of W we can say that the 
normality assumption is an appropriate спе, since the 
graphs cktained under NORMEL resuit in a straight lire and 
the value of skewness and kurtosis is arcund "0". The 
results, for example, obtained from the SC#z5#25 model gave 
us a straight line under NORMEL and the values .17 anc -.05 
fcr the skewness and kurtcsis respectively (see figures 
40a-40c). The normality assumption holds if the АБ has 
reached the steady state. At the begining cf the transient 
MSN? is nct possible to characterize its distribution, 


while when Я is in a transient state we can see that its 
. 2 , . 
distribution takes an exponential or a Х TOC TE This 15 


left little by little as W comes close to the steady state 
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EN Cchen eventually it takes the normal fcrm. The ncdel 
SC#25#95 is a representative example cf that result. In 


that model, although Y has not reached the steady state yet 


(N=10000) we can Observe the values of skewness having as 
Ко, 1.€ 1.5 and the values of kurtosis as 24.7, 17.2, 
ENDS for N-2500, 5000, 750C, 10000 resgectivelly. That 
is tne values start frcm a high level, then ccme down, they 
Fass fren the actual values of ),,), Есі ап expcnential 
model and then continue decreasing and hopefully in the 
steady state they should reach a value arcurd '0' where the 
value сЕ skewness and kurtosis fcr the normal distribution. 
The same assumptions are obtained analyzing the plcts of 
that model also. That is for N=2500 the plots are not 
linear under both EXPLT and NORMPL, for N=5000 a straight 
line appears under EXPIT, for N=7500 and 10000 the plot 
leaves graäualiy the straight line under EXFIT and begins to 


give us a straight line under NORMPL. 
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4. Тре following values for W have been crtained frcm the 
runs with various values of t and r aná for N=10000 


N t г ЕЯ ] Е] (МИМИ1) 
10000 .25 „25 „1169 .3332 
10000 . 25 .50 падащ .3332 
10000 „25 „90 . 3538 ¿3882 
10000 .25 . 95 .4075 .3332 
10000 .25 . 98 .4679 . 3332 
10000 .50 . 25 . 3831 .9985 
10000 „50 „50 . 5648 .9985 
10000 „50 „90 1.4589 .9985 
10000 „50 „95 2.1109 „9985 
1000C „50 „98 3.7196* „9985 
10000 „95 „25 10.4672 18.3528 
10000 „95 „50 13.4784% 16.3528 
10000 .95 .90 29.1312* 18.3528 
10000 „95 „95 75.1180* 18.3528 
10000 „95 „98 145.3484 x 18.3528 
10000 .99 .25 35.5986 * 53.3371 
10000 .99 .50 43. 1616* 53.3371 
10000 .99 .90 98.9556 ж 53.3371 
10000 „99 „95 139.0167* 53.3371 
10000 „99 „98 219.7920* 53.3371 


Analyzing the above results we may state conclusions 


@= LoOllcws: 


(i). Fer a given traffic intensity t the value of W is 
nct constant, but depends on the value of ccrrelaticn г, аз 
happens in the queue with dependent service. We ckserve 
also that Я increases as r increases, but the rate of 


increment is larger as r goes to 1, and very small for low 
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En en _ ян ( 2 


values c£ 5 “цегегоге 4 15 not a linear fuaction cf r. 
ШЕЕ ЦІ gives the plots of R versus r for t 5.25, .50, 
.95, .98 and also for the M/M/1 queue (r=0). Furthermcre we 
can observe that the value of W for the correlated queue, 
ccmparing it with the result for the M/M/1 queue, is not 
always greater as in the case cf dependent service, but 
depends cn tbe value ої г. It is less for rs.50 and greater 
ЕСЕ r2.90 


(ii). Fer a given correlation the value of W increases 
as t increases. Furthermore high values of + cause higt rate 
of increment, which becomes larger as the correlation 
increases. Figure 42 gives the rlots cf Я versus t for 
ЕЕ), .25, .50, 29/099 95 and .98 . Observinc that 
figure we can see also that for a given t the order o£ 
ent of W is for г - .25, .50 (0, .90, „95, .98 
respectively. Thus we can conclude that the waiting time 
fcr the cross-correlated queue is generally greater than the 


waiting time for M/M/1 queue when the correlation is high. 


5. The last result obtained from the simulation of that 
model was the convergence of Y and W. (Ihese statistics 
BE сопуегсе in value and in distribution slower then the 
same statistics in M/M/1 queue). The mcdel SC#25#98 reveals 
that Я converges faster than W (see figure 43).  Eut the 
model 5С893450 gave us the conciusicn that Я  ccnverges 
faster than W (see figure 44). Thus we can state here that 
fcr low values of t W converges faster than W and fcr high 


values of t the opposite is true. 


187 








t=.99 


E[W] 

150 
t=.95 

100 

50 
ESOO 
t=.25 
0 25 22) га L 2 
Figure 41 - ОШ UPECWITH CROSSTSCONUUSLADEUOSERVICE TIMPS AND 


НӨ? 01 INPUT.  PLOT OF E[W] VERSUS r. 


188 





Е [и] 


150 


100 


50 


Figure 42 - 
PCISSON INPUT. 


225 450 


r=98 


r=95 


r=90 


r=.50 


к=.25 


2221 UE 


SUSSECWITH CROSS-COERECATSD SERVICE TIMES AND 


PLOT OF E(W] VERSUS t. 


199 





` 9-SH 'G'i-XH '!858GCOS NOY ІНІ AD 5ЧЭ11ЧЭ1111Ч 006-4 WONT GANIFLIO “ (29) ZK ONY ZS NTIMVLIA SADNATVITAICO OFOVIIAT IHL ANF 
(-оза 949 уч) Szali ONILIVS GALYINWAD “(00908 “0051, “0009 “2052-4 ASYD 804 I ° IM P'O + ТМ) $9417 LDOHIIM BD HLIM ^SIRIL SNIIIVM 
“(h S ORY 7 с) SIRIL JDIAYIS OFLVYINNND “(Hh X ами Z х) 00001 =н ONY OOOG=EN 19 53611 ТҮЛІМРЯІНЗІМІ 431971040102 10 56КО110114165Г11 



































ІНІ. 504 SDILSILUIS 314495 40 NOILVINAVE 710181 12561923 154% 32830745 INIL IDTASIS Q31X113H02-SSONHO HII* 10200 - Ен элобта 

E xU ООП Беа = Ооо жоғ т Ма т Ба ОТЕ ПТ С Па 

ес па Т С 64<0220%0- -- 1501590 = 161025005 5889020"0-2 риса D Cree 42110 

2639252 *9 827191170 291611650 | 592921240:  61Р956670 95252180 56394195 

5211900*0 960650270 8 26120070 180160070 892660070 166160070 “AIO "15 

ра солепов © 6525000”) 2221006000 > 0661000"0. 551000 “0 15220000 4934 $ 

Lu Lo ei c. 70 1986090*0.————1£090000*0- — — 22996€0*0.———— 2065€60*0—————095€660-*0— noo — WIN 
— — SG = DIT Р ғо E О тот YYA- M : «уа с м НО Е ее На 

000007604 00 229207146 0*0 00000°05ғ оп 20000%11% 0%0 0% Ото 0*5 0*0 " 50831 
069,270 9612811 1662171  SRüRR*62 RII2C^€F — IFPIO-"^l 5951971 99581701 2620°0 84200: 029170. 681170 “501839 
1216671 - 682666 - 06856572 51(81"” — №9591°1 - 10200"59- 61681291: 11111176. 0555070 ТРЕ 96202 "35 -— 6972") =—— SS3N*198- 
01122207 901:12072-49226"” 0: 29586170 095615277557 1 б р0І?0 о ре Ро Лас ан 

1861070: 021С070 50252290  9Р80020 — 55819°0 06569070 5651040 22600940 504170 192670 226872 0599"2 (чузя 16 
R0**^2*^  88990":0 0984120 1119040 9156140. б6СЬОТО  85181*0 2198070  €€0*6991 616"? 51255599 FZR-ZELE NV QW 

сум едем сы SOE HM 2M 50425 IM со ат "s 2 Xx z x windvd 

A E ue ee БЕ -- с. 10309 2331v13w922NDn ___ BAN. ET e 
MEER D т, AA ЗА Oi ес И У ООС EEE. 7 оо Оо SCS 
>6169170- 855662170- 1€€£9205*96 901€F1*52*621 „1106 1"101 35168615 7652 “бариж! 

TM 7 Wei1680t0 777 7788880030 ——— 51368864177 7 9111268:8 аав 1 ЗӨТ?» бозчизаќ 
0600 О = P0206 S 0 020 O о 
ËI $25*000-0 »092000°C 392510070 04 6910070. »8» 2200 "0 ОТЕРФГЕ 0 ANIOS 
E£SNOOO*O $2£6000*0 254611070 6S8^810*0 6LRERLO*O ТАШЫ: музи 


» 0 2 0 


— и о а Аа A - nn ——- -- ----- ------- 


ча Y M зүйл ЕН ну ем уур [| м MIWMVd 


БЖ ЭЭ — = д -- — —— ----- 


Оо ОООО ОЕ 0702-2:003017781ЇЁ ООО ТОП ОО ПО I TE 
PUL YG Qi 23919 - 16. 574 Та BOT BEL FE IV e 255829 0111-3528 12272 22505281 :5061020- __ удао ___ 0231 -) оо 


065695 56229179  O3taG*2 099266 8686276 299276  19192'5  ЕЕВ11"6 161270 928279 5622%0- 092070- 55344136 
8966940 Ж91Р2%0 — 8325672 11992170 91:99"0 27686290 2169690 02506740 0952953 6918776 48109999: 126179» "AJO 715 
га O OE — TeREN 0. 1707070. 2505070 НЕО ода Е инт | сен 05707 — (49335 - 


£29 (E20 0606070 — 11965210 65 660*0 22015122022280190 2 02 224515220-12365/0:0----40921891.--9855С68 -4110:6699--0208”2066--П.- ма 


——%-Мм 50 з »м Ем боем га м бам ум ср + рмжю 0 Y AAA > AM Y AA Y дф- 





anand 03177135902 


ee ee a Y- a 





— — — MÀ — — —À 











Q 
• с=54 '16'2=Х3 '05%66825 809 385 40 SNOILVOIId3" 00S=0 AONA 03К1 4140 


(">33 ну ји) 51411 ONTLIVS GILYIANND * (90901 "2064 “0005 *906Z=N ISP) #04 1 
“(nh S ARY г с) саита 1214816 аліттайпо (п х ANY Z х) 0000L=N ANY OO0S-N LV SSHIL "TIVAIHNHV!ZINI 031710402 10 SNOILNII3LSII 


AHI #04 SIIASIIVIS 111445 140 RNOISVIDnGUYl *ININI NOSSIDA INV JONANDIS INIL ПОГАЧА 011УҮ1714Ч802-55082 HIIM 30350 - пБ 10514 


• (га) ZX ONY ZS NIIMLIA SIDNTETIAIGA TIOVHFAV IRL INT 
^YM 'O * IM) 52837 INOHLEIM 9 HLI* °SZKIL SRILITM 


жи | |1000000"075 00007)“ JUS 0070777 "006 6006005 "056 1715 1485 


oo Во Е 81123 898661173---- 584686616------308211514-| --4686411943---4184181651-------55341803- 
192591272 05596920 529664710 8269292: 689545 °1 5396265" ==  SS3NM39S. 

Eg oo us coc EE PE те 1256990") — 5»105502*l1 15892190701 0%:869678 952665249 ae 
Ба = FREVZIID O) 455970072 3850105"С 957 655°0 ОСТПЕПЕ "О JETSI?" (993775 
ES Ls. .5651000*0-  —...*461600*2- —  6661614*11- —. 2998661*91-— —GLZ28 0R* C3 — 3329152 0 IL — À— — (kuNv3q4- 
MEN LU eue ee reos р 0 “- #------ #6 1--- --—--— Wdvfi-—*$ M o—— yn €-M—-— — uv ат Ev 








*26*^ 00000*00G 202000*€6^» 23000”005 00032335 0090%006 29929227926 0009%208 3715 1445 








00000°S6% 00009°90S 33999°ES% 00099°00S 00300 




















“Г695656 73605686 t200871 Э6018“| 1854627 23616"2 7005170 — 22060*0 — *620*0 612000 8131"0 681170 “501803. 
-39536*L- —5€z05*1——4L2£»€*| —-£ObvE*1- L8lv**t-—-0Ll [55 * E- —356€ L* 1——-680 71 * 159 69560 *0——- REED ES IE I пару IT — 59 SNM39IS- 
90968722 66692720 СЕЕЩРТО? 62869702 54725791 15175'91 26569721 92919721 P£auy "E 9 nn 
6606047 6192991 15516") 6294599) 12729199 9946470 52595-0 1399550 809697: 9230"| | 1695701 929071 (NV34)S 
BEIDE EZ 12199°@2 6554:2:22 90996512 16921561 16228-81 ЗУРЕО"СТ 20928791 62011666 58216991 Z6[°99ç€ 1207891 МУЗН 
227 мон єй от үгт тэ c ш МА DX AAA | 











———— | 





——4Jn3525-231V3439uJ2X4D. — —————— —— mn —=—= =... 


wa o р = F+ = — — — сат ——  —_ — —— — —-- —— — и а A ка 






































2 9020022“205  2029003"234 | 2000000:026 | 0000000:2098 | 0930323:3245 | 0000032:328 | 37216 148. 
21€8180*2- 114526579 255166972 950061176 192200976 23591666 Sonny 

7 = >" =r" AI) AED... ОС 505507901 — 7 2098909515 A ZOO SS INES 
u u aM — ——--3€*6€£00*0 -—-— 232316400* 0- ———— $92 €£ 02 * 8————- «05€ 985 * .—— ——- RLESE S9 *9—— — $1 6t€807$——— 3 15 
+ _ 2. . $9L1000*9 $352020*2 L2*0L9€*0 25615860 6955162 *0 ,3629Ғ2%0 (Му34)5 | 
2822022*2- 2056022"0- 3619886751 2290152*€1 151612711 526151972 музи 

» 0 2 a ун Y A uva € M 490 см НҮЧ 1 `w — — ніне _ 

- 00900-66%-000С0%205-22222-757-20002-005-02770%75%-1/701977002-2170 Ur mT —UDUU*UUT—U0C27*305 00037375 — 3718 TIAS" 
Es Piles) 19552351202 DESI орача «из ня реет рае e ar Far A сєшаг 

1129217: 96292541 2235171 36122:1 96532241 29932210 122101 (1154171 2650*0-— 64800 — €8$22*0- — 1920*0-  SSINMINS 
760927504 20090791 27011791 92181591 ©3511 ЕЁ Папа? въ 381601 0166:01 9125796 169617 29115326) 2306062. 7030 715 

оеро Setze I Te ӨТТІ 5155 *0— —39397:*077:525 Ds IN az 


_OSZ226°AE 85459*81 3252:£*3l 66601981. 51370291 91597761 5618121 г:526:11 155^0CE€ 0196991 2579966 .120* £091 — — — МУЗА 


м ыы боем AA AAA AA 1. aaaea Ias taIa tt A 











300103 4315133435 


151 





МАНІ 05 01 51УЙ03 43 МУНЈ 5531 51 У1У0 40 (N) W39WN SHE di 
° SIXV-A IVƏILYJA G(31V232S 3399 


t 3 03114433 ЗНІ УЈ 533095 IV4HON 3H1 ONVH 34410 3H1 NO *SIXV-X 
WINOZIYOH JHL МО 431925 3389 №341 SITASIAVIS 343930 Эні “03200044 
51 531151115 833349 30 135 97 ОМ 431405 38 VIIA VIVO JHL 
JIVSN 
“NSAID ачу d^O8V ЭНА 40 H1099 Е + Я 
NJAIO  AINDO SI 3N1A-ıMs JHL с = 3 
МЗАТО AINO SI 4914 V 1 = 7» : 304 
"€ OL Т 4298434 3NIVA ЯЗ9З31ІМІ NV » 
VIVO JO 339nNN N 
АУЧЧУ X SV ATLOVX3 Q3NOJISN3WIG “534035 ТүнНОМ 9М14015 304 
0350 38 01 (7% Wad) AVAAV IVNOISNJ3WIG 319NIS AWW') V 533095 
° М А8 (03М015М3М10 
0311019 38 OL (иж 1139) VIVO 30 Ag219V 19МО15М3Н14 419М15 У X 
5 | МЗИПОНУ 


"( ч “544095 'Х )  ламзоч 11Ү2 
ЧЭМЭдОЦМе (17772 
“Al 1 ТУИН 
303 CM) 211511У15 1531 ОЧ149Н5-311М 30 SNIVWA ЗНЪ 3119402“ 


JN 
c 
°( 3) NOTLINAS NOLLAGITYLSIO WWHON QGYVGNVLS 40  3Suy3ANI 
I 


15 
ЧО 534035 IVWAON YAHLIS SA SANIId-VLVO 30 135 у 1014“ 


> OL 431441М1 SI IdWYƏN 3М110088065 
35048114 


TdWYON JNITLNOYENS 


= = ажы» Bem ао ао мо мо мо «шы» яны» «қы» «мә. а Q G Ss мо са а амь» O ны» cu w = Oe МӘР m Q um ee F  “ “ ee ee ee «нм» лшы» «лы» «ма» «шы «лын» i oe 


WOU OU OU OQ UU OO UYU OV OC OU OU OU OO VO OO VOU U OVOOOVOOVUOOUOW 


192 








LL Мау :31V0 
51М31 °M “7 “4 “4084 40 5МО1133410 УЗОМП 
SVMINVO “Г 50194039  :uaWWvW90ud 


"911594 ОЈ 
ОЗООУ ОМУ ü3'ÀIdWOD N338 SVH LI INV IISN SI LYOSXd 3NI1NOY8NS 


“03310034 SINI1NUYIINS 


“0313 ЗаХА SI 113 ЗУЗМЕЛ У 1VWHON 3uvV Viva JI 


"Ма 43 301УЛ 31 311309) Ol 30909 NI „НАМ, 
93411003975 $11 51199) WVu904dd JHL NJHL 4-3 40 Z=% JI 


1 3 d3J/53U02S IVWUON-SOILISILIVIS 930490 30 
S3IVd іні 1078 ОЕ 93119) NIH1 SI 1019: 3м11108915 NVIDOSA іні 


(1) 
Mac ер EI * (T+N)/Ī = ( X) d ЫЗ АЗА viva 40 
NOTLINAS ошл уо JSALLVIOWN)D IWOTYIdWa JHL SI A IIA 3844M 


(1) 

( ( X) 4 r J 30 3NIVA ЧНЈ 54Ј|УПЛУЛА НОТНМ AWYNANI a 

NOILONn3 WVu90Ud 341 ЭМ111ҮүЭ АЗ 13310493 34у 5301УЛ я Jm JHL 
"3М11003805 JHL NI ІМІНЗІУІС У1У3 
A8 N3AI9 AYV_ INV OST°d МО11103 OYIHA 51 *10A 5МУ1211514у15 803 

531891 УУТУ 1 ЗИ018. МӘ834 М23ЯУ1 N338 ЗЛУН 538025 1VWUON JHL 
"9350 SI aNV Od3i1ndWOO SI 
: 3 ЗНД “М"0 40350 38У (51ХУ-А НОУ) 534025 ЛУМЗОМ 3H1 


Prey CC) OOOO Ку СУС ОС ЗОО С ОК СС СОК ОССЗ ЗС КО СО 


295 





N+ K ) 


SCORES, 


X + 


( 


NORMPL 


SUBROUT INE 


a г Q^ e3CO cun - sO ~ > > 
Ni 2 OAS ST m On, » O = сет 
Ге € N OD è è oè хо eun ey > - . (N 
eno ++ ° ° ° =—-O Qe ep e AN = © ж ~ 
(N eC ° ~ - >>. СО e ОУ "б ОГ- Г” » е(7-7 ern 
Q er e > өӨГГ-ГЭэ © м = NST tm е едеми О ~ ~ eo. 
NIN e eN e mD OND PM •сј = ort o ° 4 9 «OO о = ч 
eM) ONO  ІАООО 9 от = J~ ~ > em) © mo , Q om ww 
СХ оомо eT O01 ч € ro. o Г = = +-0920190 + emt => e (N 
e N о о е ~ о ~ OMe е РА СОМ ЈСО СО ӨГ” ef = (~ • ~ ш 
С! Det rin ON em N e e NNN 9» eH) © же Ге е (Y 
Sr O» те mO D rean Te 4%» ее ЭО М Oo О 
QU) *() eum ~ & ~ • • • eer ОТ e = eg n e ç =G aQ 9е-4 U 
ФОО АЦА AIT Or e ec Oo 9 e y *e«DOON mem mo . y 
T HOM eoo-ro ey ww = е 9 МЗ Ст “ООО 9» өт ОС ° е > 
Фоо ОЛ С" 02-7 го) СОГОСЛ е(ПпСсј е eON "OD е e СЈО е е е (Де ж 
Po 0 e. [e © споју 90~ ee) CGN (Пе ӨТ > > en + ЭССО еч 
Q De eO mF è ж Cat Us ~ 9 ^ e)e N eM eM NI © 
СМ NN =D е е ° ~ ~ о, ~ ° • СО ексчбе ОечОсО e eT OF e е С 
өсү #64 СО М) егі ете еГг-7 OQ = Ht (QC © ОС ес ос — 
чИ СОС ЦА е ОМ ОСС СОС > о „Со •(~ 9 9 of e е пе j 
em oN өү = O O mu #0 • •00 N Си. N re ° e > ecco Ц. 
oU t • (11 06) UY) 99 ору # «) 96 Ф3 Nme Nam # (1) *Q7 Рај •сј 2. 
UY ec ee aN * eCO we ЧО е PHO с оф еСсЈјС СО едбуечСсСјечеч е C 
ver Ur see ем è  % +. ria СЛ Са Су е осо) ..0 .Q e ~ 
Phe) ert LF et «е Or NI ~ + #0 • • m [~ • • ee С) 
~ “Т elim оо СО е WNIONW = ° AIM e- ~ > e e эь e CN rn 
= M ӨС ONO at è è è om an NOSA Ч =< СО ее Q e«t edem ~ 
СО > e м BD ea Mn e=—— =— С) С СО ° exl ee А Og еСј е ө 
O N еч) = е есО emt Oe e C400) eC e ФИЛ! tem 9С0 9e = - 
маја WAI Om # + ФСО ~ % > С с o Ф °С еа е СО ^ >= e 
m NAND T mt = о ОТС < e (N e ° ° “= о > - ~ ~ 140 
ј СА 9СО0 өе ey: оци ~ AT NUN en со е > amili) ~ «г МЛ СС = 
шо =. УК пао Ко egy e... . uyy Orm особу СС Ort eHDORMONN 0 e 
2N e æ LN OMY OM әгтчггіе e ePTAIOW  eOCOYONIUN eg е е е о > 
O o ~ <<) вото e C CO = ° ° = aN eC) e eO1 (y O^ 09 . ет кы 
N pati gc ~ ort = > © = ~ ++ ° (ПО е е ee > ~ d Ш 
Muy ОС ОА ++ #Ф = e ecl. e& eif) (7) 9 өө! Өөр өСізіф ~ (ср ел о" ~ СС 
СО СОТ омсоссОон С = СЈО 7 oTr өү eTA ом fe > > МСОсоСсЗ . O 
ec) ° ° "СО e«Q.OQQ C == өг" mi > С) * ~ •9 ° Уф © em •со ES 
= o) Sater iaa СОЦ) о cm € е ӨГ- о f se Pp. eC) = • • о > Y 
zu) - = ° ° ~ re e (N о ~ Ped e | EM e Q ==. і С e 
- С «О = - > e() = =I) = е е ө С) ес М оС э вечс Oo st - 
Yew NAOC eO e~ = "ОСО Ya е wien ој СО • 999 ^9 вто m 
ie ОС зака МОМ = СЈО O Тесто ë Me e "Г-н е (Me Pp Q ~ 
се О er O CO I ON CO. OCO = < e= ес е еѓ- • ~ з e 9e ee em ef e |— 
DO me 90 ео ° eO + OO ° 9 ~ ec e e =p "Ni. DN e) oed WU 
OUN 4 94 ent = = ~ =@ eee N ОКТ eU * eco e 2 
= UU > eg > e ОМ) Ген Г =ч NNV ANNE tem е з e e Q 
=| ОР e е оС е e CO mm o mem СО О СО 9 9 9,4 9С0 9 9-4 9 e e - Ne 
= “и ҒИ Ос ССС = ооо ө е • • œ~ e (N e eO eft em = CD 
Qag OF си СО (~ = • - -... + > -0 ОСОСЧ СУ са 
Bec c e em NOO * ООС = = = = e ANOA (Со ҸО СУ о иОС ОСЧ ed “11 0 O OO ~ 
Z G= CO am 9 6 9e 0 P= UA CGN DONO ePMOOWSPOIMSMPMOM Ss o ert Фе °= o oè 
> ЖОШ == = e e7 NNNNA NT ӨСС о е е e ЗОО ве о ° 
Z PAW) ~ ~ е е @ @ © 9 ө ei о 0С) ч = оо -— “ = ^ ^ = @ ^ viu! 
мез үзе DOOQWMNSTMmODW = = e ec) € = ~ ~ е ој СЈ СО е СО СУ - C —Ə 
Ml смсс оО ою чоок ro ^ ^ eU) ОМ ед е ee со УУ УУ О 000 09 О 2 = 
С и) о ео оч СО TAT TINTE СА 00 Сп у 4 СО о е е е o è o о о j— 
КО o == ооо © о e Ol=. O ГП ПО ПО е о фо 90) 9 ori opi гі ~ ~ 2-1. 
k ILIN м mi pmr 9 9 9 9 9 68 9 9 ес == e © < 
VZ e^ »^ Ф® ~ == N > || 
чо - к 20 
Um 4 «1 GO 
EC oO a шо 


Ho HHH 36 36 36 36 36 36 JE IE 3E 36 Xe 46 JE dE 46 зыта ен 


194 





хо 
цэ ~ 
= 
+ ~ 
g9 ns 
м — 
m LL 
T E ~ 5 
N < МА 
"m L N Та! 
2 2 ш 24 
сэм + СЕ ~ in 
+ —ч CJ « 
~ ~ + CO - Nom 
` CN uJ + ф ° 
! a - ~ + = 
cN Z= = = = < 
pu S > - а + 
+ - » e =j 3€ 
c - - 2. ста. 
| >< > < — un UN 3 
zu Q Де (NI OO om 
~O -- = = O + ANIC 
20 m m= — 38 ОО 
~ • || ~ (NI! ~ + ~ со (~ 
ш + = -ч «1 Y ес! 
WwW x © | = со. (Mm 
LE ы u = = > = < ар ++ 
с < Oe — x е О 1752 ° 
K u. < 1 — + ° == O mm 
о 9 реч О СО == > о ~ ч Um = оа. 
п ч N + ш e © edt + 
u >< = O lI =j а 4 || ~ СО 
- | • Ш =“ X ш» cc e ~ См СЕ 
OY о ZX (N wre ,- езі со ~ e. O O 3“ 
а ~ Or N O eu) e P + +Z -- + a (С) = ог 
Му LL. ~ Ј N — СУ uu. - =) со ооо 
=) =) = =Z + SI] J == 2-1 = зе ° е 0 
I < Te ° • Шк QO +) AZ mm N 
= I ~ ш ш u 0) [er ш.) = ~ + = — O> 
> = Оп“. ес о N ONZ > + a OO © 
Om X C 6С) • mU m оч ZU + = СО. 
(9 ” ч 2 OZ СО ~ ОС 92 0) О m = > м $ 
~ eZ ° 4 27 О с Оше (О Fd °С) 
< ~ ag N ° LL. jj = + LL. í| CO ы. нұ? ~ = © | се 
ле це — 0 2-4 UNO = -Џ; СО x ч QA 
— OO o 45 e ° <[ | че «Г я. | (Ke «Тө | Те! • + 
О~ ООО Le + Im v Te ---сс со 
= ә 1142902 2+ 1 Ze wm = + | e 22 еса 
RP (yin ni | <C < — — (N w ө,ф n О eC? ст 
° O см - 97 0 9 "| — (N | == • di | Oe mà 9 еч 
WON м IZ ° ae | — — = i— C N Zum Qe = сео ` (ej 
-U) ez w“ == | өм Zj w~ OW Jj —%— N ы oe по 
ex<xZun— Cc «000 UN w NVD- ° Z wu — = • 
а =# с) OZO нишо -14UWOG m= LULL CO 59 ~ || _ ~ 
Z it~ Ec C) «ccu coco Y I caw ч. со 
= ee) Oo = ga _- 2 + = = 
— aoa f= J ш =) ш см. с) 2 D жиокк 222 =з „+ с. — ) 
„I JA O - = < м С FO < < > EO 
ESLCIuJUu.zuuu U 9 22 2) Шеик 1212 соса Z us Lu c Wd A 
C) C0 00 —— C Dm O TC) 20 очи occu се се — (уш 
X 
Чч unn =ч © no C Tr WO 
N m . N Num “Т + mo 
с 
ы) с. 


195 





PECOT 


SUBROUT INE 


SAMEX 


1 
pepe vc. 


Мо СО] ~ 
че” (1.4 ын 
Л O ee 
Ке е Pu" 
mY) » Өс = 
(ue = = 
mm CY > = 
ОКО OR 
= (2016 -. 
ч- O+ == 
mim ы  - - 
"ON = ЦМ 
NAAZ == = 
— -< Ze e e 
>< U, el <I 
PD» 
x= d 
Ц. > ~ 
e (Сү өө = = = 


-- ELY 5$ СОЖ = 


d ow 
> > > 


UJet » э СЖ 
СС ес 6 | ө ө - 
Сл сеје w => o 
QOO е-е 

Т е ото е 
eQ -- — 
— < CY j.e NN 
=Z“ ux 
— О > LL) 
> > < 

= 

r LL. 

at + De 
O2. 

cu P eu. 
СО «t «I 
созсо 
шч < 
zou E 
м Ооа < 
ЕЈ Јо О 


= .TRUE. 


SAMEX 


“N 
N 
сэ 
c 
O 
1 
= 
ГЭ 
дэ 
ще 
et 
< 
за 
2 
> e 
Z 
< е 
=j С? 
са m u 
(N «1 ° 
= || m 
==“ Die 
Por u^ 
ne) NN > 
> “0 == 
~ || < ——r 
ж ~" өэ үэ4 (1) 
РЕ “ОС 
17221755: 
< ~ 
к= СС мн  — ДЈ 
eC» e. => 
< 2 
гч Сс Ц. в 
"ч С NO = 
Z 
C C ©) 
O O O 
e) см 
= (Y 
O 


CROSS 


GRAPH(65; J) 


За 


196 


5 
5 
GO TO 45 
.2- 110 Х( 1) /ХКАКОЕ ) 


= Ф 


) 
10 


0 = 
O e 
~ O 
9 ———C e 
> w t 2 DO 
a Ye ~ 
лж ич Ow 
ч ш ш к 
e EEE ELEND — 
z~ O85 « ~ 
=” ZZ m 
Жми. + 
Ў < 3 í < 
— <í > CJ > = Н.С 
== <> 
иная 
> zn 
II Ц.АЗ« « 
~ СС ђе we O 
wen je ЕС 
<T <1 <T tU U ЫЈ 
Ц. гу си па | ва ПЕЛЕ У 
нм. ЭС >< >= >< >= 7) 


) 


TO 


1)) / YRANGE ) + 3 


( 
) 





цу 
N 
O 
> 
- -— 
C wy 
М СО N 
== A 
ys t 
e — = — 
е 
! © C 
=> + «=ч O 
~ СО ~ 
e ц 
= с 9 oL 
ЗЭ) =— _ (9 am 
~ <= == t 
I ww e om < ~ ~ 
-i ww æ <U CY: с 
EE 1 — Dom >< ч • e e 
Eu Ээ =! ЭХ LLJ =i м 
ІС" <1 UJ <и || (2 ~ = 
| СЕ. ~ HE ° eA - 
Се ~ >x < >x O— но 
(2% > f Z) LL та Н 
< ~ ~ омоси г >< 
= = е т4СУ C > = 
= ~ 69-42 ғ = > 
-ш ще + ОЧ > ~: 
OQ = O. x I ~ ~ “ЭС 
- c Рој ом Y ds pre uu 
эхо | со ОС ыы” 
ш < —5 Nt eo со >< <Ç 
. NO + ч uc > = 
~ Ec) тс чаи us 
A ~) < ep 9? 9L wU) — 
Wu ~ <I ll) en > 
цу “LO О әм e < (SC — — — 
0” моо O mmm PC UI pm) am 
~ = < (ÇN АТ ~ о || =) Z Om gem apes == 
mit ON = “~ XKAN e= lU Or Ot m 
O +O сого ~  200-с = со о 
СЈИ низ | ~ (~ ~ i Io Xx e e Z =) ны оо О 
-” фа ~ <Т 22. «ос | ~ > > e 
ен wV UOU ОЈ UIST qu. —OQQO-O 4 
Ne | ыы О У ПС ем] не ~ ~ 
x= -Z-— XxX | | се = = 


uJ P- eu. «orc ocu 09 «1 ~ 


Б ÀiLIC)UJLLI cc UD Qc 


ПИ =F ох СОКО U. ULU. Ə—F- Ə 


au 
pS O 


ZA "Y Z Prem --- 
OO tLouuodu.zxaxogocdcccuic tu 
OXO 22.2 ши. 0602Ж 505 с 


NO Um 
с” + 


nr со 
— а 


197 


~ 
>< 
с. > 
~ 

~ p ы 

м Nu 
т РЬ- 

vO — 

m" 


! y 


NT VALUE жк ж 
g 


> ~ 
N гэ 
ad em 
мез С) е 
«ч QO 
o. Mile: 
N eLL 
C == 
се O = 
uj ! 
a ul 
ог. > 1 
CJ < = 
== 
U. jt oe 
(89 Ju e 
~ u % 
e < e 
ped Ll 1 -—— оо 
< с) 
сме» жеше 
IN ODU 
"Zo >> «іц 
= о“ (3 J 
- o хоч 
> — © IOG 


ed eU) mst СО 
{= | 
“Сос <чщш++ > 
MANONG > өө 


e Qm X 
ХО Qr ege 
У + ох ж 

— = цд го 


от соо N 

хх ы, e 

me X ISSN 
(О мм ~ 

a uas 27 еб uas uas es 
е 


== eee 
II ILAIA 
== E222 
ccc Coo 
ao 5090936969 
HEU. EIE 


32 + 
Om ODOT 
оо OOOO 
CIN TODD 


END 


Х, М, К, Wy A 


( 


WILK 


SUBROUT INE 


INVNRM 
D1(159), WD2(2CS), WD3(209), HD4(8), 


C 
N 


~ ~ 
СО) «= 

ц 
eod 
СУ м) 


| ~ r 


ex 
> > ~ 


Y) онь онь 
eel 
“~ ~ СО 


= <<< >= 


А 
М 
САТА 





ме ls x Ооо а 


- ~ ы» ~ оо ~ ~ бо е е е 
Sa Q O rq I NIAN rm 1 09 e 
ODT ADP СО (SY €N О СЧ е4 СУ 
ОСО СОМ Т Сені ссе 
(eI mures oocooodo 

• . е ое @ е е е o 0 е o 
> P> Фе те ~ ~ ~ ~ ~ ~ 
€ QU eA O00 P eO О С С 
CO 4 uU «Q (uf e CO) «Q9 umm 
CIM МООС um «run 
ОЧМ тт счеч О О О О О 

9? e 9? е 6 е % өө Ө е ө о 9 
er » 9? оо ~ ~ ~ ~ ~ = = 
LALA PO OM ен ТОТ 
r= Го СО Ге CN (OY T — O = DIO 
СО DFU et UV. OOO Or 
Q OOOI OCO CN —— OCO OO 

$$ 9999999992999 
Pr pp ph Ф» % ~ ~ е е ~ m 
WO NOU Q^ ure СО Од —+ 0499 ej T 
Qr-o' gu'guvuomor--rg 
ог от сосчососв ее 
Ne QA (NO NN == #4 СО СОСО 

6 Ө 09 e 9 9 9 6 е 9 9 өө 
9 о. 9 ^"^» ф ~ ~ ~ ~ ~ ^ 
РА ОМ оо = O CA O eun 
c q (NGO O CO СО О OO CUA Os 
+ сч РУ СО O O (X Q C сс 
МО со оне eACO em ST € N — m mM А 

€ 9 оо оо e оо е е е ө 
~ ~ ~ ~ ~ фе е е ~ фе ~ P” ~ 
rn OC со г СОН ООС 
HOW eat UIP UI СО ПАТ АЈ е 
van RMN Or TN 
(NJ L'AXACNICNJ UA =4 O (Š ) CNJ — 4 —4 = 

Қ в @0@ е 6 в 8 00 o ө 9 ө 9 
9? 9 9, 9 е е ~ ~ = ~ >» ~ ~ 
DO UY O GN UA CO U) Ch ++ De 
SY Qe 0l 00 040 Мөн д O O 
МОСС О М сә Ш eA MO 
„( Су Ој СО O = O ++ CN CN ча 

e è % ө өө е 6 ао oo ooo 
Ф Фф. 9 % %9 9 6 % ~ EP = = ~ 
P- 0e G^ P- 9 Q^ 00 СО е OU 
r + O uy CN (N O O = f СО cO +f 
SOP COS et et CU) = e 
e ql UN NY O (NI AO MOS INN СМ 

e 9 Ө ә 8 ọ 6 86 9 È o o е 
re er PF FFF ~ ~ = % -~ ~ 
САФ CY rn O О ДУО ЛЕ ићи СО ХР 
мост O- fto rnc mop OM 
GO СО 00 CU C Un T f LN LA LA 
O Q e0 0 T rO rN CN CN 

9 ө 8 9 % д © @ о ө ө ө 9 

~ е е ~ ~ 9» 9? 9» = ~ ~ P» = 
РАСУ СОЈА с0 (~ су СО 0 О о 
MLACIO © Ш с сп О ОСЧ Сх 

QOCQ со Q0 m0 e m O 
(Мен Су СО СОК го СО СО 


~ - ы оқ 
MON O = QN 
соге см с о\ 
SÍ OA OY 
O Q e СЧ 
"а , во» 
~ ~ = ~ ~ 
СО 50 ~ со 
er «Tuo 
WO со гп ос 
O OQ e С 
e ө 99 9 
re ф Ч о 
СО ОД «rete 
CN O OLUN ON 
СО СО) Mea СЧ 
Q e m 0J«T 
% 9 9 е % 
~ ~ = е е 
|D CK OCO ГТ 
SY iy СО 
Ссс сог со 
el el HN O 
"е е o 9 
re fF a a 
(Oen Cu 6100 
QO «OQ u^ C^ 
СМС н СУ СМІ 
ete СЈ СО 
9 өӨ 6 о ө 
~ @ ^ e ~ 
DONUT 
СОСИМЈСЈОЈ 
ш со uU со > 
m =i (N „ф СО 
t o 9 е 9 
= = ЮҠ" ~ ~ 
QJ 00200 
CM UO e P ОД 
CO e OO eun 
= (N cy CÓ O 
е . 9 е 4 
~ ~ 9 ~ ~ 
СО сс 40 02 
~ со QO un t 
е ЦА rO (Y t= 
NOOO 


Cr FPF е ы о 9% Ф о FF PF ~ ~ P” ” e" 
Di СО СОМ А ДУ LA MIA OR O NM 
ИО СООО СЈО чл СО (N CN + OO 
ANOS TINY SY GM 4| uvuegquna«sonmuveom 
O [Or NOA NO" NO TH Омо чо ч 
@ @ % в е е ө 9 9 6 со өз о % 9 9 9 ө 
~ pp Өө Ф4 ф ~ ~ 9» Ө о. е ~ ~ » фе. ~ % ~ ~ 
СО с eq] un TY rr nuy n O u un moo «T e 
СМО ЈУ 0 «QUO Ч0 гс wf СО (N = Cry ССГ» {СО 
MO Y СО XP ар СО СО ЧУО О о ~ о ил 
O СО NO "MO NO өвч Dra NO rn Oe 
e o 9 68 è 9 Ө 9 Ө oo % 9 * о ое 9 өө ө 
* "^ "^o^ Ф 9 ~ ~ 9. = ~ ~ ~ ~ > 
СОО СОГ ОСС Опе С СО со 
99 со Г 69 СО СО DA QO cr — CO О -+ О 
Y OD О ЧО еј VOM СЧС со ош +) 
Оччо. О_О. ООСО ч 
9 9 9 9? 9 9 9? 9? 9 9 9? ө 9? Ө 9? % 0 9 9 
- = ~ $$ ее ~ ~ ~ е е о. 9 ~ ое ~ ~ ~ 
СОМ ~ Ce ССО О ЧО (nG ОЈ па У) 
ЦП ~ ОИ 02 O A ALA OWN 
O CJ e О CO CO — Q" =+? (3 O C" C) D СО 
O (N CO CO — СО лет DO MN Ur DINO rt 
e e 9 9 9 в ө 0 o Ө 6 9 98 Ө 9 Ө 9 9 % 
Фф» ~ ~ 9» T ~ ~ ~ о. 9 ~ ~ ~ ~ о. е = ~ ~ 
C UU O^ enm СЧС МУ О О со С^ сп со е Сч 
VAR СО СО СО МО СООСУ Сосо ОШ чи ССО С 
DEIS TUD A DA MO PI ч О К = O 
Q — (N СО СО Сы С СОСМ СО ч С\) СО яс + ООСЈ 
6 e 9? 9 өөө è 9 өз ө 9 в 9? @ 09 ое @ 
>”... ор оо ~ ~ ~ ~ ое ~ ~ ^ 
CN QO 9 UP 00 eUY Om mU е = NO <? ШАУ ОГУ 
С) ОУС СЎЛ ЦУС SPIN Of DF" DM WN Or uh 
ох бо анг ог оо бич о СО СО 
= p CNJ (2 — CN СО =4СМ O AN O СО е (КО O CN 
е ое @ 9 0 9 6 9 9 9? 0 е Ө 9 ге 9 
9 9 о = е ~ е ы % ~ ~ ~ ~ pa. ~ ~ ~ е 
cN OQ ЧАЛМА O O Td eco) ex «qe 00 
PMNS OP n O r Q (9 + СЛ СО (~ GOD 
eO 6106 (NI CO tn rO OO MORAN e DO м 
ті + Се СС ес С) по $ OMHMOHOON 
e 9 9.0 9099 9.0 % е 0 со е 9 9 0 
= $» ~ ~ Oo ^» ~ » ~ ы 9 е ое ~ os ~ 
Ln OO `Q Os G" O C: uA QO O QO I + O += Оса O O 
nee mio S— ОД е ' tn (N =O €N O O +f 
CO AC wes +f OTO + O O ay Се О 
e 0400 00 e СО еа Онч О О СЈО mm O m0 


e o 68 6 6 Ө 68 6 9 Ө 09 Ө 6 9 68 о 9 9 % 


Pb %» ое ~ о @” ~ 
СО QUA UY (NO O ON 
OO Koh drm uu 
"ЧЧОООООООС 
Се СОС) е (СО СО СОС) 
e е 6  Ө 9 9 @ 
9 9» е » е & е е 
MAO OM Doni 
DANTOTO 
(CN CJ e — С) СОС? С 
O m O e Q eQ e 
8 68 9? е е е % 
9? 9 9» ~ ~ ес % о 
осе хом о 
r ON (N .O о О ~ ~ 
MONA St tet 
Q e СО ен СО ее СО ее 
9 9 % 9 ° е е ө 
Фе = 9» ~ ~ = ~ 
мо ЦП СО СО СО СО е 0 
Po (a DD Pu el CO 
гус сс СМ СМ СМ СМ! 
Cu O rr O et CU) тэ 
• . @ ө ө ө ө ө 
фФ»оф үт 9? ф оф Ө 9 
Шчо Си О О 
Ге СЛ ее СО ЦА СО“ ем 
лос 
OO O m O ч 
6 e Ө 6 Ө 8 9 9 
~ ы» е е е е е ~ 
гг час СА uA 
г Хлг 
QO OQ U^ SPIN GI 
СО өн СО өвч СО өө4 СО чі 
6 e 68 е е е е @ 
Pf е ~ ~ ~ ~ ~ ~ 
NIN em ч О 
ODO OM MAM D 
О 0 ока о: о 
ON O eO eO e 
6 8 Ө Өө 6 9 9 9 
e фр ~ ~ ^ e ~ ~ 
pt LO et CO СО СО С) ил 
NO NDDON 
co rS r— O .o Q tr 
ONONO mC -+ 


- ж Ф % %7Қ, 9» F?o о„ eo > P A e 


са со Тоа 


QM СГМУ спс сч О сч сМ С 


ГУ Ш 


+ +n f. NADAN QO О КМУ ил ОЈ СА СПО С Осо с ~ ИО 
пеела СО че СО 40 а0 б0 50 оси о осо бог го О 
rN r GO O X МО СО пса Он О © С) ос О ен Су ЧО) ЗЕеч СО СЈ СО СЈО СЧЈ 


• o • • • 
ref ~ ~ ~ 
см roc! 
о 000 42 
СОЈ СИ ч 
сс СО о 


m @ 9P » ~ ~ ~ ~ ~ ~ ~ ~ ~ p е = е 
Реф ефи О СО Qn Ol СМ СО Од СЈО КСС O дсасч 
ооо те Qo О Осо «uu esf ФО еее ЧО чмо ч а ТОСТ 
Qf eU Рио А СО Я О N CIUS OI e CJ em CO) eme sd C O 00 eU UO FIO OO DO OO 
г" че уен О ЈУ У ОГ О ОДО e «NT DO ANO ANO NO A TONO чо < XO £010 04 C009 


€ «9 9.09 ө e 0 ө е е оС) е е 


Ө а 9$ 9 999 9 909 9 0 9 6 9 Ө 9 9 е 
те @% @ 9? ~ ~ p ~ ~ ~ >» oe» ое ~ е ее 
+ г O O O ЧА СУ MM Mm OO OOOO 
О ++ СО О К- 0 ЧО СО СО е-4 00 00 Сп ~ О со О ~ + 
cog mod] Q0 0) -O 000 er- со Ne 
NU СО A DO OO mr OA NO ч О СО 
6 o 9 oe 9 9 9 99 0 е 9 @ 9% e • е е 9 
~ ~ е е оо ~ ~ ~ е ~ е pp ~ Ф% ~ е ~ ~ 


ОУ спин О СОМ Оч ИМ О 


36 36 36 36 35 35 35 3 € 3€ 48 3t 3 30 30 30 36 40 30 3 +T 4 Tc 3 4640 4$ 2€ 36 3 3€ 6 АНЕ > 


198 


... . э . % e 
~ ~ ~ фе не 
CO f Oh СПИСА 
(N = O O (N (N ч 
е“ Ол ОД О со а К- СО 
“NN СО САЈ Ср СО СЧЈ 
ооо о e 9? 0 © 
> 9 9» е ~ е m 


DIO FM С 


e eO • ө ое е е е а е 


кек 3€ 33€ 3€ 





т.е е ~ ~ ~ = ~ - от. 
OND AIR OP Онч 0 
сс СО СОМОЛ С СС (У О 
O T O O == CN AT TOD 
MO O m O чс ч О чо 
è 9 9? д 686 6 Ө 6 98 9 9 
~ = ~ ~ љ~ е 9» о. е ое 
NO PODIA MN егі 
OP HSH 00 € 00 COLO e СО 
СФО чеч СМСС Ц СО My 
CQ eO eO eO MO M0 
оо. в 6 0 6Ө 9 0 8 9 
ea 0 ес от е фе ~ ~ ~ 
mM Om O~ M = çN QQ Nu 
ООСО Осме С опхо 
чече ИЧ ОИ ОМ 
ог [697 С) ее СО е-4 СО СЈС) 
ee в 0 68 9 0. 90 4 9 
> е ~ ^» 9» p» >> ~ е е 
NINO М О М) 
T Om NDO ANATOA Ww WF 
QUAIGE FTV FD ORIG 
C) =: C. tO Mt) tN 2 
б е 9 è 8 ев е е е ө 9$ 
~ ор. ~ = ~ F&F FF KP ~ ~ 
DA OLA mM СО СО OO ++ сч 
C eU LO «Q OS «Qn СО 
MFMA ОО О ИУ ОХ 
Qe O0 -O-OOoOQwNO 
е 6 6 өӨ 9 9 9 9 9 9 @ 
Ф $9 9 ~ е = Qe. е  % ~ 
NIM (DUN OUND eat Со ен СО 
єММТ СҮСГ-Т ЧУЧ ouo 
FLAT DINO Q (N OO > O 
СО е СО) ич СО) сед СОСЧ СО СО ее 
9 9 e @ 9 9 9 ө 9 9 е 
~ ~ 9- ое е е е е е ~ ~ 
r GCN QA Q (n С Дуе Са 
‘ADNAN TAD e P + 
+. обоје], По O Qe 
О СО ОМОСОО о 
9 6 8 Ө 68 9 9? o 09 е г 
ет ф ~ ~ 90 » ~ ~ е ~ е 
+ см СЈО чм О (тусч 
O Y DL 4:00 e otf 
NOD OA NUDO O my 
O “O NOCUOM=O ч 


.0314,.0244,.0114,.0104,.С035 / 


~ = оф ее ее ое оо. C 


MO + О У СО ЦАР 


9 


I in G G: O O O G O == хо 
O O О сч ооо О Q-—oOCa ou 
ONONONOO 0 я: О 


"е е е 9 9? е е е е 
~ - ~ ~ ^» ~ ~ = = е ~ 
IANO OND ч\с О 
D OG (NCA O Q a O осо 
= qr ORD UNT 
осо СЈО (m Q СО e O0 e 
ое e 68 e. е 9 ọ е е o 
9 9 O6 о. 9» 0 ое о. % 


Ошооо юш сч JT 


14 0459, 


e 


DOPE OP O^ CO e «QU «wo 
СО О ООСС О О со со ыға 
ON O MO O ri 0) Orm < O 
9? 9 9 9 e е е е е е е СО e 


HH HH Tc ЗЕ ЗЕ ЗЕ T€ ЗС зе 


+ 


— 5 
er ~ 
~ С 
UN 5 — 
ee ~ 
= = 
eco + _ 
— өне = 
СО ~ ~ 
о O 
NO ~ - 
-— E < 
= — = 
== - 2 
+ < 
-— e ше 52 
~ | о 
ж 27 + 
Гл amu a. == 
еміп = © 
ee eee ФО ч (N 
“(ту (О е 
че O Zj = 
<= — O LL. сү 
pet са С? 
Om = „2? 
"ОХСООО 
4 e cc СО 
са~+ > o ~ 
== 2 ZN 
с = De 
~“ ~ Цу РХ о 
“+ << 
Om] mr Zu 
W WOO NO Fe 
O WIN и ОО<СсЕ 
= eN 1 Сос 
ы гс Ш С 
- 2 «С mei 
< СЈО) „zn II 
p< і < ccm amem. 
=) Ша! Hog ee 
I LL. Q e me 
LU »-«XQC соса 


e] 
өте 


199 


Al 
e 
хо 
e 
О 
=. 
* 
4 
Еі 
о 
СОТ д 
eco. 
О •0% 
OQO x 
H о 
! H 
= | 
2N «ч 
MS) CD = 
те 


А(Т) = WD(CIA*I) 


) 


50 


ЕР 


М 


NU 


22 CON 


) 


Х(1) 
591, 


Х(М-141) - 
TEST VALUE 
52, P. 


V. 


В + A(1)%*( 


— SUM )*x?2 
В = 

К А 

o R 


) 
) 


= > v) 
— м 
x 
9 PU) 
Z 1) 
= 
+ • 
= 
NIE Ne 
4” сл ~ eww 
= LCS OD 
e I IA — 
zt > 
HW OAL COLD «1 


= 
се 
E) 
ke © 
UJ 
cc ud 





t=f 
НЕ ег РЕ ще 5 УПОЗ OILSILVAS 33430 E 
JH) 40 SN IVA о ЧНЈи JVAL 19VW3 ЭНІ М 93578 0350 51 М%%%%%2 
“1-1 “((14Г-М) /ТЭМП5 ҮЗПА04 АН 538035 *NOdX3 3O NOILVO1VA3 303 


°ЅІХУ-А IVILLYVIA JHL 
NO 09031955 ОМУ 031У711УЛЗ 33У S3402S *NOdX3 JHL GNVH YSHLO 3H1 NO 
"SIXV-X 1V1N07180H dHl NO Q3]!V2S duWV N3dL SOLILISIIVAIS 93040 JHL 
“МЗЛ19 51 5311511У15 33080 ЧО 135 У 2МУ 031805 ЭЯ TUM VIVO JHL 
39vsn 
(€ МУА1 83197349 349 1508) VIVO 40 YSEWAN N 
"SI AVuuV X SV ATIL2VX3 GANOISNSWIG “538025 *NOdX3 9М140165 
Ч04 0450 38 01 (5«1934) АУЗЗУ 1VNOISN3WIG 319NIS AWWnü V 839095 
“М Ag Q3NOISN3WIG 
0311019 33 Ol VIVO 30 ($x*71V33) АУЗУУ IWNOTLNSWIG JTONIS У X 


5 | МЧИПОЧУ 


(М “5380025 'X ) 119Х3 TVD 
S39N3hnO3S 9NI TIVO 


"CVW4AV9 *IVAWVO9 SYU3JLINVSVA FHL SIAID ONV SILVWI1SI *2 
5382395 1У11М3М04Х3 S^ V1VO 30 13S V S1O TN 17IdX3 ANILNOYGNS “I 


35048114 


114Х3 3М1100388315 


VODOOODOOOOOOOOOOOOOOO DO OO OODOOOOOOOOOSOCDOVOOOOOOOCSOO 


200 





— < a TEE ae aw cr ow ee ЧӘЧӘ V еко СО доо е ер и? «тр еи «ие» < ано ако ако «ығ» “ВЭЭ G ако вээ ано ам s s s s cA ано ано ee сии ано е ме s s Gs ано ss s s G PU ако s AA A ne чи ам вэ ӨВ ам s ee өшше 


IL За :31У0 
SIM3! *M *V *d 4084 30 SNOJ193YI0 язамп 
SVMINVO °F 50198039 :U3WWVu 90 Md 


° ç NVH1 SS31 SI N ШЕ 
"ЗПЛУЛ INVISNOD JAW VIVO “а 
"N'***'2'T-] AW3A3 ЭУ 03437 МУНЈ 55317 SI x I 
> ЧЕ 03193IdX3 SI 1П44ПО ON 
SNOLLIGNOD 30933 


*9115 Зи 01 
G300V QNV 93711340) №339 SVH 11 INV 02SN SI ,190SXd, ЗМЪ4ПОНЯЙ5 


03Y1n03Y SINILNOYANS 


7041234Х3 51 114 МУЗМЕЛ У GilifnSINLlSIO ATIVILIN3NOdX3 34V VIVO JI 


"544095 “МОЧХ4-591151ЈУ45 33030 
JU SYIVdA 341 1014 01 0311993 SI „10144, 3М110038015 Wvu909Ud JHL 


Es 
С 


-( (2-N¥(T-N) ) / (X- НО (£+(2-N)xN) ) 


¿ VWNV 9 


(1) тес 
O сорал AS X)WDS +N ) 
3 


5 УТПИНОЗ 341 ЭМ15П ОЗ1УПЛУЛЗ 349 2VWWVO *IVAWVO9 SIJLINVIVA JHL 
* 9391109131510 *NOdX3 AT!v3u ЧНУ VIVO 


UVWWV9 


0000000000 00000000 00000 00000000000 0000 VN Y 000DODoOQO 


201 





SCORES; M ) 


Xy 


( 


PAPET 


SUM2, SUM3, SUM^4;, БТЕ, ХМЗ: XM4 


SUMI, 


SOBRES 


ТҺЕ ЕХРОМ. 


+ X(T) 
SCCERES (1-1) + 1./(XN-1+1) 


— || = 


»xXe U 


~ ec 
«uo O 


We DO 


Пис 17,19; 


zo 


) 


PARAMETEKS 


SCORES; N 


САММА1 АМО САММА 2 


1 


( 


EPLOT 


Н 
при 
= 


NM << 
=> 


MOT 

+ ++ + 
= 
«tU. uu 
LLJ pes 4 pand 
= CO 
> 


+ + + 


$ 
Z NOS 
“>= 3 7 2) 
ww LY) 
1 >< 
IN H H 


аа | 


— | ni 
с 2 < 
| ~ >< = 
„т 
XNA < 
A ti y 
- 7 -+ ° > 
—4 >< wet ent CN e 
| £ I % £ < 
ZA ZN E 
— >< (N >XxX í x Z= 
|- рч со 
===> >> 
x XNXX 
че “ан Ji 
+ 3 ~ ~ 
~ с СА ОС ~ 
ze IMO 
NIF>IZO 
2 3 жжжсч 
at ж > 


» 
LUC02Z ни 


IMFO же ~ 


Ш =i! 
Ne > > D = 
NO NANA lr a 


(6,400) 
(6,401) 
(6,402) 


Z ANS 


SUBROUT INE 
REAL *8 
CCMPUTE 


«1 «i ac ucc Lu cc LL o 
CQ O Zc = = =< = 


= N im со «и 
= ~ N (t) (9 


22 ФШ 
под o 


C > 5: < 
Ç > >< >< 


2300 
KNN 


CALL 
COMPUTE 


2 OVO O 


2197 





? 


= а 
хоо! 


САММА1 
B 
Xx 


SCORES, М ) 


? 


E PEDI 


SUBROUT INE 


- 
e 
eo 
= 
+= 
= 
жә 
= 
= 
= 
e 
~ 
ge 
ж 
~ 
= 
Фэ 


>» 


= u- 
Г ne 
~ e 
v» > = 
H г--- 
m ео & 
o > > 
ке e o 
e SO „ч 


— eo = 
СУ ем м - » 
пе 27 27 s 
меч) с o Yet 
Natit ae e = 
~ СО СО ~ ~ 
>< LL. o = 
<- = ~ о om 
xz===—- = = = 
LOOD о > > > 
> ОСе - - = 
— =» NA — 
Z — = = е = = 
= С) > + + e 
Пе =». 
WHOO+NYX 
С РССССе то = 
nvo нэг 
0 .-- 
Der = = ee ~ (АД 
Л СС» =w = = 


-0.«1ї «146 
pre NN 
га 
~ 00 
x ч? 

= 
~ мн = LL. 
== + De 
O ч 22 
m (SE 
отада 
ZOU 
Шые < < 
ОКО F 
молат ч 
aon С) 
+ 3t 3* 


) GRAPEH(CI, J)-2 CROSS 


2 


„Ед. 


МОС(Ј,10) 


203 


CROSS 


— 
= 


СКАРН(Т,1) 


) 


ш 
= 
= 
< 
e 
> 
м 
om LL 
20-52 
— <= 
NA 
ucc 
Се >< 
О 
Ш М 
Un 
| 
эв 
= zum | 
2 ~ 
чиг 2 
o» t 
LLJ ~ 
ne, Y Y 
=! О 
WY (9 
ш aw Sard 
СЕ 
= ++ + 
eC о 
907) аа С) Qe 
ede Dev 
u зень Q 
xZNN ++0O 
ч” 
уши; ми 
12:92 eo 
бс Z ecole 
2 Са < ~ ~ 7) 
~) СС СС - 
XU» || a 
ZZ 
J H H UP No me et TU 
С 
ши< < ин<г 
OO FO a € 
ZI Lr COP 
<< ш ш Z 
Сс осо ES) 
>< = > >= Q (2 
ха 
ч 


—$ 
ITE (6,201) N 


— 
=- 


17! 
WR 


< 
-- 
= 
ы 
52 
> 
55 
CJ 
! 
57 
еттей 
ж 
< 
> 
| 
vn» 
es 
x 
«t1 
> 
—Q 
цэ 
эс 
ГҮ • 
o с 
- Ш 
rii ° 
+ 
1| беј 
tp 
| ~ 
Н 
Са” ша 
(Y) aq =ч 
C 
e 





N 
rd 
ri 
e 
= 
II v 
"э = 
= >С 
- «t 
ж- CN >= 
“«м еч 
~ m 
ЦУ m ын 
ж = = mj aa 
ШЕ "oy 
EO CIO 
= < а 
~ СЕ ^ = 
Оњо — 0 
CD 3t >. 
< `- > 
= —4 ] 
=) ~ j 
+ |) T m 
Ee Gc 
(AAN «д 
- с 
эхо 
u < 
Z > u 
Z 
~ LA 
om (/} хо 
(um ~ 
«о ©. 
FOO ос 
(64 TU 
GN = > © 
O Orco 
È Ц СО ан 
ELT I 
Жен = Z 
U << O x U — 
ин у=» — (” -Z — | 
= 
O 
e 
и сп 
4 Пп 


АХТ5(Ј-1) * 10*XDELTA 


) 
à 
X 


que a 
er ~ || 


>< (NI 


) 


| H 


IN 
— = 


оа 
х 


< O 
xc 


44 


uy 
СМ 
C 
= 
° Та! 
(9 см 
C 
ae p- 
• ы 
SY ~ со 
СО x 
ш м 
@ — |] 
йі 00 — 
= 
...<1 po ~ 
ш с 
CD << • e 
tm LL me 
zwei Wo > 
SI “эрэ” c e = 
<<Оғ ш | 
wr че D? “иц > 
>x — 2 < = 
оси. на - 
22-12: = — 
e O < и > 
СЕ be Z ~ ~ 
Q< СО 60007 
ceo = G >< 
VA. |CO < 
eq) e Xu >< 
Cd += ~ La 
.ф 9 Of (Y) 
eu О) “о 
ee (SC. e> = 


eu ALU ADA 
21) 2.1 • UL 


eur 9 


mb 


ZOO FO = 


„rm Zum. 
zz «ozai 


ІШ. 


> 


Ox 
OO 
Јо о 


€ о. 


-—— Qu. s —OO.Q9-O0 
>< l|» >< Z — — чи чи 


Nr) 


tá pe = w wa ых ччы [1] С) (11 


oo ш 
"(с < 


UL LL be bb fF 
Se re ст 


ПС ОШ аса 
— Z Z Z rmm UE OSI 


no 
cN + 


204 


2X; 


м > œ 
~ env) 
== 
© ји 5—4 
у га 
=m 
е 
|| -- 
~ ~ 
92 (а 
СЕ em 
m= Q e 
< Q) 
Q. LL et 
e LL) 
с - 
и. ~ 
СЕ а 
us $ 
a 
СЕ || 
O e 
= 
ц. х o 
O > w 
— = | е 
NX æ 
LL! = 06 > 
Ico 
CN Z ~ М ||] •• 
гм N 
— Z "< «uj 
= ~ 
~ = TY 
>< “om $j 
re €) mx) 
Э 111! 4» i 
"ОС ечи > 
MAOONA. ee 
vr OU mr 


40 ӘӘ c0 oe 


Zee rere 
<< «<< << 
Зээ 
сс Coca 
шоог єє. 
ши цици 





СР 


SUBROUT INE 


PURPOSE 


15 ТМГЕМОЕО ТО 
ORDER A SET OF DATA IN ASCENDING ORDER. 


LIST 


S UB ROUT INE 


l. 


= <> О 


- >< 


rr 
== 


еке еі) 
о сова 
а.а. 


CALL ING SEQUENSES 


Xı N). 


( 


LIST 


CALL 


ARGUMENTS 


NAL ARRAY OF DATA (REAL ха). 
INTS. 


COO 


x 


205 


USAGE 


CRCER AND THEN 


INCREAS ING 


EN. 


O 
V 


LLJ 
= 
E 
40 
> 
ә 
< = 2 
= =Z < 
< — >=> 
e © 
Q. uL 
Q O 
ш ч 
C wu == 
ul QQ «IC 
Q Z Qm 
OL eul — 
ONY Та 
шә Oe 
La) az 
DIT UII 
«to YN 
© > шш. 
Li UE Qa 
armo а 
202 LL: LD 
Тэ» Фа 
x CQ vet 
о. ER 
= ш га 
лао uu 
mJ U m 
rau СТ 
OC ша, 
со ци са 
u) Су, 
ШТШ uJc5 
тегт I 
= се e< 
и 
шы >ш шмш) 
> > «>> 
љо (П) р 11 | эмч ен 
MO) 
| ты 
укы уу о је 
aa JS J) ° 


St et 
$ ee Sse 
2602 = 
22 222 ~Ј 
маға O 
DI DE DIDI О) 
Jl 22123 
асоб асве і 
(2 QH-O]Qui- 
о «ч a 
ооо 
сосе Fu 
cem T rr tu. 
LD TOO = 
ни оо 


= 
шии шш 
Wege lt 
--н ос 


SUBROUT INES REQUIRED 


сеосоеос ооо оо ооо ооо со ооо ооо зо чо ооо ооо ооо ооо 





S31111 119139233 


(1)X = 401 
бс DO ( NIWOYd *31* XVWOUd ) 31 
N/*I = NIWONd 82 
ALLVIGVdOYd 3HLÀ HdVu9 ОМУ 
'S3ISN3n2339U4 *ua3gWüN IVIJIS INIA аму GNI 
02 01 09 
а2 01 09 (Nas NE i 
( 80Sd'XvVWOud ) IXVWV = XVWOUd 
N / (1МП021719014 - 8084 01 
JNNILNOD zZ 
1 + ANNO) = LNNOD 
ОТ 01 09 ( dOL °19° (1)X ) 31 
МЕ = 1 22 30 
(117Х «431 
С + INNOI 02 
"0 = XVHOYA 
1 = 11 
г уд Наунд 341 31у25 Dl A3N3QDO39394 WüWIXVW ONIJ 
(NT X) 13Y0SXd 11У2 
1 = WANYIS 
(МХ NOISNSWIG 
/ зжъ / ЧУ 45 “170 
Нү15 Тж 1v91901 
WONYSS “1 мпо02  Y3IDILNI 


эр л» е е ЧЫР s G “s s O «вт ЧМЩР ç. > mu ¿e O ам Q w = 


(м *х ) ISKI INT LNOYENS 


SIM31 °M “у °d deste NOL oad) GO. ам 
SVMINVO °F 50193049 = AIWWVYIDLI 


0000 


000 


лк QUO 


206 





6,200) 


thea 

> 
= се. О 
met EO 


TO 15 


GO 


ы 


1 NSTAR) 


> 
> 
СЕ 
ч 
үнэ 
2 
w 
> 
ae 
а) с) 
ES 
сс с. 
Q. 
[ 
- Lm 
- 2 
ЕУ 
вар С) 
=- 0 O 
oO 
цу ~ «2 
° e. ~ QX 
-- > 72 
+ + > Lr 
== Ш 
= ге с. 
a x ш 
> ш ©) 
O wu ~ 
Эд w~ ж ~ 
Сб ж 
ын | — 
м со Ра. 
Е =) > + 
— Ç ~ —) = Z Q) = 
=O Z mc 
=Z * O Out) 
20 ш N+ ° 
Oy • C 


< Z Z O OLLIƏ 
QO a < + < + ПО 2 
- н к = c 
uU. осо ü. 

12-2- 9 


! = - 


СС weer LL Ú J жь 
otf = == 
Om m ma 


X(1) 


! 


Р 


LA 


© 
> 


WU OU elu etl 
Ze Zum Z СО 


aee reat 
LA 
= 


Іс-і.М ) 


Х(1), 


( 


(6,600) 


ч 
~ 
ч 
- 
== 
х а. of 
од - 
Quo UL 
<> Ie 
- — Ш 
ч < 
к>” > = 
Im ei 
Са СО 
-НУ > 
Omt Ow 
ша 
XI“ uU 
шеа me 
Сз СО ОО +» = ~ 
te c» 
Qa. uu vx 
ве ° ох 
н „г У С > 
< ц. 
9 00 Nu 
ae He Су 
= >< > Z < 
Pa prat quo рээ 


FORMAT 
ж 


О 
e 
N 


207 


== 


>= == 
ЕК СС СЕ 
COO 
u. u. u. 


Q eO 
сос 
Fro 


END 





015. 
OR € 
OR T 
BAST 
ED F 
(REAL*4) CF DATA. 


IS 1МТЕМОЕС ТО 
I S 
I 
I 
Š 
к ) 


T 
S 
S 
A 
Ха Ny 


SECTN 
OF 
ОР 
ОР 
IM 
IC 
( 


SECTN 
E A 
E A 
E А 
Е E 
B T 
SECTN 
A SINGLE DIMENSION ARRAY 


SUBROUT INE 
О 
T 
T 
T 
E 


SUBROUT INE 
PURPOSE 
2 N 
C Е 
С Е 
С Е 
T T 
CALLING SEQUENCFES 
CALL 
ARGUMENTS 
X 


CQOOOQOGOQGOCOCOCOQOQQOOQOQOGCOCCGOQCGCCGCOCOCGCCGOQOCOCQQOOCOQCOCQ(QGQGQOCGCCOCCOQOCCQ(Q COCCO 


208 


NUMBER OF DATA. 


N 
K 


MUST BE NO GREATER THAN 100 ). 


( 


NUMBER OF SECTIONS. 


USAGE 


МЕЈО 
Nu Wi 


510 ) 


( 


AVERAGE OF EACH SECTICN 


MEAN 


ECTICN 0.9. 


DD 
S 


MEDIAN 


VARIANCE THE UNBIASED ESTIMATE HAS BEEN USED. 


STD. 


Р: ЭТЭ?) 


( 


HAS BEEN USED. 


TEH UNBIASED ESTIMATE 


DEV. 





~ | 

c «t uJ | 

p ш - In © i 

uJ - < ro a | 

v) = LO <. < | 

а. O Or | 

uJ = LL v с | 

х O c X De = | 

— O шев Wik < | 

m D ше << | 

= ш zou Cr | 

uJ a 2U> uu с v I 

< «OO rm uJ == | 

~ < OQ u = | 

= ° L Y) ш к шш =з ш | 

2. - =- Q Tau e с: = | 

uL = нэ Бш c = | 

>. ш Се = Z3 D2 O ° | 

22 = O v) LL fe — O) т | 

229 C ш. — MC = | 

= = wur Lo = ° | 

= e ч En аш ш « | 

a = uj == н (3 LL uu | 

се < m и) = • Uu со • | 

= Се = 1» MI a } 

— Z = Q- LL) “ша Wr Vv) | 
= ш хо >. > са E < | 
«ч <) ш С? 23 Осе <t E ° t 
+ O e аб -хо O U. i 
= O «1 «22 "Uu. p О і 
— Y T LLI NN m —ч ЁС } 
- - а 4 ша «но uu о. | 

Ex II a < Lu сас шам • с) } 
Ere по a ee X a p= Z N Ue } 

> е ш іш eu Б utu еа с) < © uU | 
88 © LL bh Od LO dawn > | x | 
С >го жш м1 Ош Z -— We М - - Y | 
шко xc = о DOO m~r 2. ТҰ 217122 | 

е Z ur NE © cr YU 9 < © | 

на vo О = Е 2 Ми LL. = О ке | 
< te Zu =O о аше = — = | 
= Q < qu << >> wO v) . O | 
V) «wn» m | zu тич = c | — CM | 
ES c.t ФШ м OTOKO ш v СЕ | 

| EN I < = Fe D ul OC = Be 9 mm | 
Dr > < nO UN Jn < Q = F а | 

ОЈ ша Ц won ша = Oa ш ui а — ps | 
С мин NY ~ ~ > Uu Let C a > 0O O E > | 
= шаз ul ~ X OL Wi DIM sec >< ш uv) Се ш | 
Се IT ~ ~ OT zo m /ul* ul A ж m Ee o SET 
м ХО > < X Le Ск но мо - а ш Z < | 
a 5 оч Tul uU) > = O 53 2 | 
< | И lI M “O Ты сою m C | 
> че му Zoo ш ш | 
v» — Ju = > а. к (бё A ve ve і 

BE 2 mi WN Joc MO MIE о = | 
O (7 “9 >=> > >— < кеш Izzo а ul E Er | 
uJ uU) D D I T Эш ez VOU r а 200 ui | 

. = С) mc. Y) TI ll = мд Ors = | 

ш 2 > к. < Z= WOT шш © EF Am = | 
ш ш x< Z x utr ul Mar «ос ока < ! 
Шо ~ = ~“ < а > C О 2 O Da x | 
о мо x = X ко uz «o шол) Z о ME O ш | 
EE = == со O fa | 

u YN «<<< ји ја 22 © a < | 

2 2 a= == v? > Q. a T 

| 

| 

| 

| 

| 


ооОоооооооооооооооооомеоооосооомобеоемеомеомеомеомоме 


209 





X3 ND, K ) 
(100) 


( 


SECTA 


SUBROUTINE 


ш 


bm == 
NO e 
Смс > 
YN)» 2 Ге 
аж = 
em СО 
== "Os 
ZUL m 
~ — < 
xO Or co 
«Im 
mz 
€ Ux 
= J «t1 
SO a 
ай” 
< Z 
мес) O 
ше < 
по офа 
-—CJLLICO «t 
ee oo 


101 


= <(— O< ОС H e 
Ше о e e pp ee 
= =] 
бб» о» "o X i 
о о еве е = (N 
>> оо хо 
э а е - 
Q ожзы 
сү © 
~ са 
• 
LL; ° ° 
- kk ud 
• © zl 
e e 
Qe 22 
Zee NYS 
LT YOK 
H 
"I 
uL eu. x ш. 
вх — >< Z m 
38 4t 3* 3€ 3€ 38 


e Р 
Z 
e (N әл. 
— > 
= e 
O О wn 
а = Z 
я e 
| a 
= > 
с> N 
«T u | 
LLi 
~ (N 
СЕ 
CJ am 
[ш 
-- “м 
2 N 
ы O e 
< м- >. 
к - + ~ 
< Y) — e 
e мын Q + 
= <= (N 
e < - О 
LL! | же 
= 2 + ww 
= - х 
ке4 са ~ N 
= L O + 
O c = > 
шш — 4 === 
= LU мм е. (_) 
= (0 = ~" 2 
x бж w 
U. ша “~ X > 
Z < ú DA aO 
UN Qe AN СО 
С) к-р- » A= N 
но E eu. utt) 
> nor O Or 
и) > ШО ~ es 0 
— = 6 x х "ZO 
= Coe | © 
Tez ~ ~ » | 
mx due e X 
WY ur Fe 
C ЕО wan = 2 Z 
ш СКЕЧ > > > Y Y м СИНИ == (Y) 
б< c 266 e— - 
m2. oe ~~ ww at Ye Xm If e 
Yee We WOOD | ZU ANO ul 
We то. 1| 7249 =) oY) am __) 
ON = СМСС wo DXN е» 
ш ~ ~ ~ 772 ! =. > 
ш РМ РЕМ НО E ex 
= Q -—— `< 
= — — == (NI CN || jl j СУ < О. < > 
б..) Ғ => = >< e 16) Lour 
Фа Uc OODONWNN -4<1 
оа ыы 6626422000 «t = uu. 
QO (m = 222 2~ С C =m 
== 
=] 
ОРО 


210 


STATISTICS FOR EACH ESTIMATED STATISTIC 


OO 


TAT(T,J) 


>< 





= ~ 
~ > же >< om 
со ~ = со - - eut 
2 VY) Un - e Oe 
es d ЕС < © ~ СО 
|| cN сч = uy = LLJ su 
м = z = -Ж е L ner 
e Q mL ~ Wow < * e« 
» CJ Q XA = • 2. - 2. 
мез > < с» е ег 
— - | | ~ ~ чи) V) men 
e D - жə Lanz ~: uJ 
B ~ eN N ZTL к ~ zt 
cN че <I >> ШО » «тыш гт 
C p зө? cere > wile SMU I u 
x -- e Core а =ош ш әш 
= ~ « ~ > ~ WW <f >< => Oo ага 
T- e c = > < < > 7. c x 
Cx. WY p amn ше = om ие ee LJ Wi CM) tu 
C Sc ~ Y = ue <ы сз: Boo 
VY) ~ e a — Li) elit æ U DU ИЕ с «ч «X 
= a e 2221.) «я V CUM. шаш 
+ | ~ + || + Th LX .Oe ac ao 
~ СО O < = = e > ee zr 
~ N ~ © = o 22 Qu »x e. LL] T =>xX << AI 
OO ~ ~ - че om = Сте.2 PruxNzZOFZOO 
шие ж 0 72 XC x ~ Wt sa «ШІ y «C 
- A > > ~ „u «DI NĒ» ee (20J—-]O0 
£L Е ~ «че а. Ф O." + Qo Ze my) = œ o OU a j Z LU 
с ғ | > 22 "ne N Ze om Wore *"NOZ LD» < 
e OD" 7n lI ~ “је - ~ ~ < > < CO. OA < X 2>2—— J UvuJ < CD 
x мо ~ ~ N (N Gd YN мше) асыш носи 
- WY e || O ev OO Ii X OE Ede OVU- mje 
” >< ne ~ 2 = 2272 Ol = e Ce ONO UFO 
Со Ух 727” - ~ A) O a uJ © WN М» thet eL) 
см > - >x + KOO & Б eO =m > Z O € O r е uu 
= ax Om m > = х-- = FR oe «е --ж-ц очи ТООСОО, 
HOZ SY wam ЕС инь “зэ? унд =D) зэ ZU twuguzxEeuxzouu 2 
Q - E? == < > = wil] e —Oe O. Z—ü e е С Za» 
«е |” < + = NI = N re a ep "O = > С 2 00 2 
Et m Me =x Zak æ IW Ne == —2 ное еконо О 
се | vx < сом Or X= WOK wf шет ос = 2 У) 
(2) -- ---- а u >C) em Ske © Wee KMOe СТ МАДА uiu ши 
UU eJ гэ ma e exo <a N шпалы a о и 
= Ч >< — e = E ræ NUV) N es ea NNK IVN 
u Ne =j ej mp ~ | Оо жа ес Z > эт. {м 
ч” NN N e ж v» =ч О >ш 3 (M9 #4 00 ЗЕ ШЈ 2) о Ш 
— + о ~ е ХР Ос Конче РК ОК Ее ре 
<t- |l Ш О ~ ~ ~ nn Dur 21048 em ie LL oY) Ud 
za = O O ~ uie "e Хо С.ш 
um m Г-45 O SS o (Y) ed “. “с. xXXN ma "xs «ч 
A om CQ m | — эр ~ е “(ү б ч М М С СС OA 
+)» = ОО О т ол О СО OC е IN е“ е ~ ~ ~ ОМОТ »2 < 
шо. ә) ООС | зечОСО О оО оо ОБС WPO ee Y DALIA АА > 0 за Ае а 
Y = АСКО Н Уә мс N NN tO MU a d O O ч» СО SNe o Ne - 
Y Aena € e - e»x(N е > оС е ге w we тий ен w x 
N че о ~ 0 N piwo oO оса. ~ DO = JO 
— 5 Б -- =. ! E o s w eet Se 
“I ad od < O Z Z x = X zm | | ЕН 
че WW ree a и ao o LL СО << «<<< < << 
“хо онно әм онон он nme e DZ > Z" = zz 
jet fd (nt p= je pag jae pant 1) T MAR LXI ~ | се СС СА СК СЕ Ww a СК. 
ca а С) cx: ud a ud oc «t — Oe Wa De LO eC Amer) 
2320 zur x < = CO c Z Z < Z CO G —w— Zoll u. LL. LL цш чи 
++ +++ кк + $ ++ 
съ f с © u^ O © Ф a NO 
съ + = ст = O OOO о од 
e NIT O о оо 


211! 





y 
= u. да: 
40 
= E 
| СЕ ea 
РШ i 
Nao 
= ќе 
E ~ 
= - 
N 
5% m! 
B = 
jt LL! Ж С 
> -, Хз 
= АТ we Jt 
eO. > ~ %2 
~ • ~ оо = 
чш (N -a са ~ ~ w KN 
+ LL и = жі ме < > ~ 
'= а E 2 өгөе [P же 
ща гэг H- ~N mj 
em) m e = << — | = 
y ~ ~ N нэ me = У< • 
де >< =m = 2. unu» Qu» < 5 гай! 
о.) + Чо = гос | 
e mm ЈЕ Ғ- e Un ~ ~ ~ і < — (G (N 
-х «с ve ир. UJ чи ЗЕ ЗЕ 
ev) ru: - e Ec cm O OS 
vo mt J 400) і — ++ Са = | ` + — a 
Or „m - =: - "бл + Z | те 
UJ «t «t -tu. = = ~ — oS) Ds A = 
спот e - KRIS OW we ee 
к а ++ x ~ ~ ir фана amn OJ ~ 
ROD It t= Ке СЈО ХЕ чи ча Z < чо | СМ | 
СЈО о - ZO - <( 3 Jt + x= “~j Zi IA 
EP ec < ~ | 46 36 З ete Ye 226 6 Б 
NJ me со + Z = NULL U. <> li НО че YN) Y) 
Z= < < > © i DK „ч — — — — < = nn 
Set ZO = = r „ХХ + — FOCO, = 2) = ~~~ 
- "Ох - =~ — СУ! AN e jj — > 
N= y = > = <=... ~ + + Zo ~ + — r + 
eur =O v» ~ < + D —— — Н НОС МН 
x OX Ww = ~ Ze xx = + cuc + v) O CO «tW — > >< 
«АЦ. Ц ч М. ww LY) -2. 4 5 Z = ~ |- ~ о < 
~ ДУ ~ x “=> | | H "12222001! 92025. Pi 
х оос х ; Сме “СО 2 соо NN е = —— — 3 o» 
- 34 c | emn NOOO ии DY momen we YP en 
~ —.—— = = ° = | = СО О“ е е е НІ ІШ n О Х = СО OPM 
== O о. = ооо ~ > w e e eyaz ae 
~ =) | 4907007 NO Бо АЛ | Оо 
= == =) NR ZUNZ XZ a H || LL Z Z Z < —X—— o км? 
5 == се ш, I REE NMP чоо ә о FR | ЕК | | 22 
ос e cC C) аз Y AZ || < << < = = > < < < мг << (О 
O Ou = (1 2 СС) --".2222 U || је 22 ҥ„ 02 
Це LL. LL. LL 9? OXOMXO NNNNMNMN ОС NNN жж паз 
жж + + 
“Т О = N 
о OO = N 
ко, со оо 
Q 


212 





T OF DATA. 


Т TES 
E 
E JACKNIFED ESTIMATES. 


au 


R 
F 
T 


IS INTENDED TO 
Е 
L 
T 
fs 
A 


JACK 
D 
r 
F 
I 
O 


JACK 
E 
Е 
T 


SUB ROUT INE 
T 
T 
T 
E 


SUBROUT INE 
PURPOSE 

C 

C 

C 

C 

G 


ОКО СзО Ооо КО Ооо 


CALL ING SEQUENCES 


) 


XS, STAT; М, IG 


X 


( 


J ACK 


CALL 


ARGUMEN TS 


CF DATA. 


(REAL *4) 


A SINGLE DIMENSTON ARRAY 
A SINGLE DIMENSION ARRAY 


X 


XS 
STAT 


215 


RETURNS CRDERED DATA. 


(REAL *4). 


TWO DIMENSION DUMMY ARRAY (REAL*4). 


СІМЕКСІСКЕРС ВҮ (16,7) 


NUMBER OF DATA. 


N 


IG 


ES STATES 


NO 


=R+1). 


INE ( 


NUMBER OF GROUPS PLUS 


USAGE 


5 $0.1 


P. 


( 


AVERAGE OF EACH SECTICN 


MEAN 


ECTICN O.W. 


000 
Н $ 





LL AVW : 31v0 

SIM3! М ту *d 4014 3) SNDIJD3YIO YIGNN 
SVAINVO “Г 50184039 = YINWVYIJOYA 
91 154и 0| 


З300у ОМУ 031144092 Ч3348 5УН 11 2МУ 03450 51 1905 Xd». ЗМЕЈПОЧЧП 5 
51М3Н331П0334 3NI 1nougns 


"МЗЛТО АЛМО 39 ТИМ 
УЈУО ОЗЧПОЧОМП JO S31VWIIS3 E=>( (1-9I1)/N)x (2-91) YO 2=>91 31 


ет>м 4І 0З1934Х3 51 1й424100 ON 
5МОТА ТОМО Ј 30343 


"ВАЛО 13431 VIVO ANV УМ1ОЧУ25ТО 9010149 4715 ЈУПОЗ ЗНЈ ОЈМ! VIVO 
ӘНІ S32V1d 32Vf %0308У2510 38 21 3AVH ТИМ JVH1L SLINIOd VIVO 30 
YIJWAN JHL JAZINININ JL SV YAGWAN v HORS 39 ISNW T-9I=Y : ЗЛОМ 
“МАЛ 19 0513У SI OILSIIVIS H2V3 303 NOILVIA3O 
“0МУ15 ОМУ ЗЭМУТУУЛ ЗНА ОМУ QOHI3N 3J3INDDIOVí A8 V1VGO 40 53211511У15 
(515501 Ч)М “59534 MINS” A3O 539415 '3ONVIUVA 'NvIGA3W *NV3W) 21598 
Эні SALVWILSS MOVE “SINIVA-00N3Sd 3AO08V SV оз4а9й02 ан 9NISR 
(SVIAWYOS 3WVS djHl ONISO) VivO 30 
135 3311М3 3Н1 303 GILNdWOD ЭЧИ SORISILVIS ЗЛОЧУ SV 3WVS JHL 
“ЈМЗКОН ЛУЧ1МЗО Н18034 23H1 SI YW 333HM 
en xRJINVIYVA / YW = SISSOLYNY 
03104409 А115У18МЯ 
1N3WOW ПУЧЈМЗО ОЧТНЈ JHL SI си  343НМ 
exx” AJO "ОУЈ5 / tW - SS3NM33S 


! 


INV3W| /*A3Q "0715 ЗОМУІЗУЛА JO “ 4309 
( €45*d ) *üG3sn u338 SvH 31V4I11S3 O3SVvIS8NO H3J1 *A3Q “015 
“0351 м339 SVH SLVWILSS G3SVI8gND 3Hdl 32ONVISVA 


OOOO DUO UOU UU OOD O00 0 ДО СО ООСС ЗОО ДОО О ОЈ оС) оС) е) Сре) е) е) е) е) е) Сре) Сре) 


214 





XS, STAT, ND, IG ) 


X1 


( 


JACK 


) 


SUBROUT INE 


- т т зе ^ 
ecc е NI 
- ee Фт е Y) 21 
~ = E = 
= LLJ> < ~ | а 
- те е е е е см CJ 
а > 
ж - ~ ee» s | 
Ф os oe oo o N 
eQu uu) = a 
60 е е е е = 
~ ~ 
y е е ф» a 
= се е еее = p 
e ZO WV == 
= т т е е е + E 
= с = 
* ^^» a» ~ =) + 
= ee se е «wm ать Zen N 
т = өт ә Фф с т Qu ”? = 
Б Oo = x ~ 
Lal ^» ое %9» + савж, Pa Uu» 
en = <le о о Ф е + ” = > 
~ Ze m oL Ze- ч Z ze BE <> 
> т оф т о р е = == ~ - 
e - < > il с 
pos ел ое ое o С С) e С) = 
o = oe Ф = ee - a x — me (y 
< өт т р е а ее ы ~ О ЕЗ eC Uu = 
- Та о = nn Z Rx = 
“9 “o ^ ^ 9^ 9% eed > - (5 AU N iet < ж Y) 
т = өт өт э = (= be mm >< 
- ше х-ожос с = I ea NN = 
= еше---- и O = о - = EN ыг 
e > Z < < — ~ eO e 
Z — тото те aC) Y =) — LLJ O fon) NNO t 
~ Сје е ее е е е ( = N СЭ. (23 maxX >< > e 
~ ЕРО ~ u цэ JU | eC) 
Счет ое е е ао ее Q са t O v зам теі „= 
се > N Z ~ ж | 
--- == me Y) 2 о —_ em =) | эгч = 
—2 ~ ш же ш — N Фе Du Nee ~ 
Ow N es Е- Fr NN m „с 22 7) Om LU =» 
=< < < - © Cd m de ONZA ~ 
“сү . °. > о > Z =Z + Jr нн-—ун — 
E Fr w SES мој CU) ти < X — — >< CX || 
= ә I = ми > E ME I II 70.04 E 
с ° Та uuo WM O erc eNO ll ш о eu lue 
EX = å =x WU «x ш МЕ | МЕ ~ с ~ = LN 
mca. СО | Ne = ~ CANN =e NZU чирот, = 
tj ZOmHOxX ш ш = l! — = — _J || HW in Al A A D 
L ьч ZO FP Ке — 2 -ы < w 
EO ~- ~ Z 2 = Y] Ne. w Ji CN == C ЧАО w 
LL] — «T | || а 4 6. : 2 I! је JO ODON MAX m 
ск it = ee ee NO NO „a 
"O < LL XL O о < tC NOVLU NLOQO O ЧЕ 
Eo мег O O O AZE J ДЈ О O QU 
+ 3¢ 36 ~ + + 
~ - СП 
N = сс 


OOO ооо 


215 


WRITE (6,20C) 





=1, IR ) 


Мұр” 1 


Js 
rf 


ССМРОТЕ JACKKNIFE ESTIMATES. 


1X4 


жаны 
223 ` 
= 
m > C 
=. ~ = шы 
! JO e « 
= см = =. 
= > < э 
e C - ==) -- 
- O се се 22) 
22 > “ < ы 
Фо { > >> ~ 
= (N - = = = 
= = • Y”) мч шш 
Ре чи? ZU Ср = 42 
< «tui ж-е m (Ту 
N > ын — — 02 - 
ЈЕ 22) eu а Lu e" 
* o am Жы > = XM 
— = Ae 4 < 
= ~ -- г LUZ keL 
= e — » х eu < + 
— ~ + д >= Qe ак 
к? e e >< DW ZA JÉN 
е e © шо en «шмо aq 46 
-- - - 22 =ч» ИМ ам СЕХ 
<L ~e CX ч ~ - U ee ee). NN e 
bb | - Из Хөг» ШРАМ - ш 
< СС а. >< «< шш. о 7 omar 
Se pm mt mi = ~ стон ахеч о wa 
=y) j uo m n Ф. аш и VO 7 нас 
~ + Y сс >=. — c а. Ne С] Фе m Ex) 
| — — — (77 ~ || e QUO enx OO e 
Cnet м, - ОО NS CO er x (~ uox 
== | | ~ ~ оо c LU >< (CAAC х эмо 
a d ~ СС CC CC 7» = > > Z = Орде. IUe Z 
Es | w = ~ мо O wen Im LN Ae A DOI + 
+ — w (N) x — и N) Ze ro өт “(Сү LU] e 
— I Ку ~ == xe ms X (iO zum O 
29 дет 34. = > ZXO Ear «шын ДО ДЕШ 
> C5 > «X се шш mimi NI -жадр-це 3 >00 
С -p e oi j < E eN LLUQ = = eraot 
— mann al = = „Ze erui ошаш ошо 
- om at pt OL me (a di Q < OW e eXe СО) е cc 
— ч «ын ~ | ~ == eC Lr ev ZZ € XO € 436 36 000 
<{ кекене 2200 < очи? eA очаг 
- а = оаа не m 2 О ено DJXwv DO 34 UD 
v») poet и | |-- С. СУ + қ аға» — e$ “соц ов ~ Цече = Fe 
em e NN) OC ~ = s= сода Z ñ 4 (С e U 
II =ч) с + "ч O ас = a > e > ex == = 
Ре +MY КП ~ О Zi ~ ~ Рош > ен Жар << >< 
e | 4 0.07 —о 0“ 99.9 * e > әз. x eC. Oed e 
m, - QU ammm О OW чоор AN «лж NNR s< 
= s ra = Н. 1 077272 ноосоо O O a ха о О А 
| < ОФ I ee ej ~) "је СО Ne eff Nea “со ``. - 
— CN e (NN ~ > ~ XN e o ч” ` м, ww "uM Gra 
9 ~ ионы” НЫ) о ONE ma~ (.  ~~О x 
= 11 оон ı ана Ü (D N 
I ini «I«I«r Oz Z w = Ze Е-Е- a Pp = 
fk HOE rem lite ON el < < < III < 
да. а. С ANY) Dk De со < 2-Б О э р >>> > 
---- — <= | СЕ СС ек СС ОССЕ 20 
с. О шсш а ЧО DO Socia 
(ЕЗ C С с кен чоо ец. uU цн uU 
+ + + LAS 
ХО T pa ut цу со © O r ОО mr с 
О г ~ ms N со © оо ОООО 
N + $ SQ со СО со 


216 





* 
803 FORMAT 


IS, 16, 


STAT, 


X; 


( 


JACKES 


SLBROUT INE 


SUMI, SLM2, SUM3, SUNG 


16,7) 
МА, 


= || 
Hee I CN MT 


= 

— >< 

— 

w pt 

<> 

Ф 2 

eu) 

= 

т) 11 
Oi OOOO 
СМ ку 
• Печ o o è 
C sO OO 


eA) 
zei u 
4 — w 


= 15 5.5. 


«09026 Б 2025929 
><) оду 


= 
тч 


217 


QUE 
= СЕ 
-Q * 
NNO 


) 


) 
(ХМ-215(ХМ-3) 


) 
)*6 XN-31 
XN* 


X N— 2 
/ ( 


* ( 
) 


) 


(БТАТ(1541) 


= >< ° 
Qu X jt (y) 
СО СО ers ~ CN] i 
| < ~ КСЈ 
U + 3 
INA 3 3€ 
= | N mt un онь 
az im 
oT Xmas кє 
Ке ew хо 
€ — C JO ч че 
~em | (NI þe E 
жи га ха 
~ > +» 
OY) 3t (0 
2 Z "o 
= Jt = L N 
= + mn (Mar 
Panmezz 
«T u^ -жж 
pe #6705 < 
NNT ни 
DEN 
-Nt k mm 


r+ mn te о ~ 


- = єг 


NNN 
ЈГ 
CP Sa s 
= j- 
< < <1 
rer 
NNN 


Ф Өз 


< = ж | о 
v» ма” 
-— ІІ ~ СЕ 
--> 
ог << — 
с :2 БЕША 
æ жх спаси 
ЗЕ x 





TO SIMULATE 


GRAM. 


OOQ 


~ 
~ > > ы 
oor м © 
Dom ~ и 
ооч Y e 
DOn e w ~ 
"тд 2 ч” 
-Q00 с O 
>x = =. r = 
= >< > - 5 
ш лс? о 
Ф pa) O © 
om Y Pray C ES 
5.4.) ome om ud e Li У 
ФООНЖО • © 
сососоцаАч ФО 
O O O < uN 
"O Xu С = = 


e > оСУуСУЦА ее 
Ç) O O mO < MX” 2 O 
оО Јо 2 9 О 
€ CO C UC) xu 
«44 414. ZO е 
IN=ZWOO w 
292 Mm e 1) 
~ UN enw m > 
>< = ami, (onem) outs 
Gam - ”"ONOO е 
о--посео- 
eC CC uma 
=O COCO е е гО 
СО Ср чо уч от $ + Ç 
IH OW) we Oe we LO 
чао реј Се со ХО е 
> 
о лл 3 26 о 02 5) ~ 
Qu. n uo x 


Е Ех 


SUMPS1(1), 


SUMES(1); 


) 


1(1); 
8*1.0 


/ 


it 


e Ç 
=< ~ a 
N ÇV) ы 
LEN 
æ (2 > 
>< 
Ф > Фа. 
N e "O 
W Xx x 
— (Y Q e 
WY 
> > ra. 
= > >< C 
x Y Са 
е“ $ 
~ == gums 
ом > mx 
m (9 С) (20 
ӨЗ С? 1 Е 
=ч (Хх. Cc 
e po 
- (7 »” Y e 
ох х 
ш О О 
— ” СГ e. 
шы 24 |-/в4 
~ х >< 
x x er Y e 
=>) ==) ж = жж. — am 
== ЕО wy) — = 
x с KO m om = = = 
e d. T U) UN = >< >< = 
e x A UM. LLN m ммм 
x = ~ x "ер мм ~ "d 
2 =— ~ mex e eg) uem e € eX mm 
æ N N NY со М | 1 eue 20-31 1 > 
x OO wow | i рече СО | | pt url end 
> T С =m ӨС я -- С) > 
ох» м БОСС Ор = = O х өү. 
HUC DN өс” mm Pj = = жә” 
mart Om XL eo mam ^^. (С) меа ф 
DW ANCO Dr N =l 
ree CL Се (бучие че 2/7) HU) O suo mo E pd e >< 
==>. + = ANQ VND xXx “же 
NX же e Cg dz do az 
DUu EAS ex = = I = 
pued pd ZA. LOL ar a Y) чичо СП) mwa о HY) 
= -| 
= ~ ний яав “Зэ? ~ 
э ° 
= iy сав = è a “ — = = @ = == — = =a 
CO QO = O PANY SHY RH SONY) ая 
O O о СО „,-JO02ı : OOIIOOA Оо 
.. со СО (М e«NnNF — MN е УСУ ONE 
~ > ~ — ~ elije е оШј ~ ejjj ~ ell] 
MN м е OH го уто Теа [721 го хо СО АП О ОС ШО ФО 
U ч UJ < ш мм Uj ë œ (1) 
O са CIE Qm. CM. Or. ОУ 
< ада а н See «ч Ty 
uj ш ш се шие Ш A LL. 2 Qc <í LL. u J < 
с сс с Ж кс созса О НОС ЈЕ СО c = C 
2 цэ цу цу 
с са 


218 


ОХ, ІНІ, 151, ІХІ, ІҒ5, ІР51, 


10Х 1, 


? 


ох 
эч. 


~ ~ 
a= 
“уж 
za 


u 


О 


KRITE (6,201) J, 


45 J 





К51, кх, кх1, QS, GSi, QX Qx1 ) 


KS; 


( 


“Т 
~ 

ОСМЕ м 

Nit + 

< ZZ Z= 

= 

| H H i Mg 


= CN (Yt >C 
x 


N KMN 


СИ си С 0 
У 


e e» o 
рев дов ОС = 


a V XK 
ожоо 


— 

+ 

> 
hen 

MaS cH. 
сс а <¿ <] < 

00000 00 


© 
= 


= 


11 


m un. 
~ С) О 

~ > 
SS = = ) 

N >< 
e ell LL 
DI > trad w 
LU > > 
pay Pt peed qud 


e (С) С) 
N>TIKH 
с) О 6С СЕ 
ще Да 
сж ее 

o 9 
Ф» ww quel 
J 
(3 е = 

d 
” Фе) ж 
40. О. 
Qu Gh >< >< 
>< >< LLILLI 
шы. 

~ > 
^ 9v еі 
Un» > > 
LL. LL. LL. LL. 


= 

e о e 
Ge OQ 
LL GO LL LL, 
1! 9 o 

e o 
5 a 
WKY >x O 
CO OOO CO tO 


w w at 1“ (С) 


- 


uuu gu c 
bur һе ме іне ме ( 7) 


210 


~ - ~ ~ 
= мм = 
2у N x ж 
U. шш ш ш 
ca e w u 


AUTOR 
AUTOR 
AUTOR 
AUTOR 


ES E 
шиг шэг 21252722 
ч «ч «ч < 
С СООСУ 
>» әң, SEA èë 
O m O œ 


Се CJ C C3 
11.) еі! == LL) — LL] — 
€ ° ° @ 

-л ә кж e 
со vu x х 
а. а. а. eQ. e 
OD Сене јен 

Uu x х 

Uu uou ОН: 
Om 900 04 9 4 
Aaa 3 
E AAA 
«I eT ex e € e 
02 N ож ож 

OO UD 

TE ХЕ 


e "9 —.) ы ә ~ 
Q ONG «ОХ 
UJ em UL ә Чч ЦС? 

© e 9 e 

om e ~ e 

— de e Je J 
NAY) AU 
со Зза. Оа оа. 

x >x х > 

w LL j= LLi = UL: < LL I 


° è • 
4 OO 9 •0 • 
OO оо o 

ии l 
H i ии il 

2. = > 
> > > > = == < 
mt CY. rm C — LET < 
QC UA Qc uu «t z «ta 
чыч. ж.ж 





0. 


lI 


WMBAR 


ut 
po 
О 
p= 
O 
O 
rot 
° 
2 m (Y 
Ep ш 
аните ° 
mle Т 
= ч.4. Q 
ПЛО о за 
I ! < < Oç 
" 
һе O H H H— 
ое. 
-44-4С)99 ы 


моо шш ч. 
мед X LX e > po 


CJ 
e 


е Y) > x 
mi a CJ aam riam Y ars 
YN x 
Q = шш. Q wr a 
mm d == ік 
w. ~ w w 
w<< == ~ w a 
[ZN =ч срж m 
шар ~ xD = 
с о < r~ СС 2 «1 = 
< ule (7 < Wwe X 
ule e U, Lue e U. 
O O 
- 1 - 1 
— 9-4. - oe) J 
EJ = ыы >3 
<<< >> М!) + <> LU. 
ОРО = 22» 2 < 
er UY) ee ж. 
= = e ге 
a Ul — Хх 
e wo x a 
= > 2222-20-22 
ooa O оо o 
9 Y) 
9 ° = 9 9 o > • 
| ee Cc B eme СУ 
DD ш очож ш 
e) Uu. 0 ° >x eLL 9 
и u 
== = e e m ~ e 
©) en Y) > == х 
eld» Ne МУ У 
N 9 ao x x A 
ша > ы C x 
= 676 ett o — em ell 9 
= UO Ow QUO С 
=== ғ > Ге = 
NA еі IX а оде 4 
A tet 00) Дай 9w4 9x, 9 
x ny е x > 0 
LJ со Ш са E 
O e e« e O == e< € 
+ e UD + - х» 
e o %0) ° ве ес.) e 
=QG =D +С Се с ФО 
-" We Ve U We Ш 
~ 9 9 One e ° = tz. 
= тэг U: el (энэн Li! 
NN WM ВС X X едем 
MA LAND eQ. =Q. = >< > CO 
OMY O Ом Gr 
tl e xu H >x жи 
wor = 1| — x< w ~ (О 
Ж ж = Жою 
LJ ш шш] ш 
zu ц орган Huc 


Пе tee KD 


220 


> = 
= = 
° о. ~ а 
= x "4 x 
Y ш x uJ 
ај а) 
> e 
° =i co ° > СО 
=| эв? = ов 
e > = moi 9» оты 
> сл = — ES >< sd = 
Y) >< 
> | фә LL. ~ | e» tL 
4 — =; = 
v = | >< = И 
~ u) ж > ua » 
Z wv = < > = 
~ Y) c uL > > I wu 
— om zZ шкі еш» 22 
N ~ Оо х ~ > 
om Y Qe. > 
— m CJ ~ - СУ 
-— — мо чь 
9 >< 
<> ° O <45 9 < 
= ых < >= јео < 
с чї) > ° сечи 2 ° 
СО , се С dO. qe dg 
LJ «t tu Lin < LU 
>» я ° ede le 9 
e uw) == — 2 == 
— <í = x — (7 -<1-.., o x 
el > 7 (4.43 м - > X< Hal 
< LU оч < W -1224 
Dem ° = > O en” e = > 
m.) < = 9 m= > < = ° 
NL Qro XZ Vm 
0259< Z cox «t z 
za ° mm <і — > e «T 
~ = << је ож _је ° 
Nu = < < 22188) „~ e 
e = > O = > > 
> 4 O Dei ot O Dea 
1.4 • e $ = ° e 
= U) = с 947734 е KK ext e 
Que но OC сошеш Fo а 
е чи 9 (ox Eu 9 = o ens = Ш 
о өш » >= e—uJe 
4.1707 4 Gc e LX e uc 
97920592 902 Ли ежжож ох жж 
“0. GOX a. CO YO Gdore С 
>< II WY il >< і >< || 
wm ~ || "UL „Не Би wr 
= = = = 
uJ - 143 = ul = LL) 


i 2u 2 сипати Cu 
ыыы Ke De Ху» 


3t 
LA 
me 


+ 





е 
сс 
О 
• 
Q 
Z 
0 
СЭ 
ы 
e 
= 
e 
— er 
m СО 
“ир? o 
-4 e 
> O 
s O. см 
Qo. I! < 
u. 
i = ° 
ІК Си 
N F ~ ана 
(N чаъ е СО 
= <(— = 
~~Ош | HO 
i j р = 
`w чыыр >< e 
SO -- X OQ 
KMD | > |- 02 
ага. ае. C 
ZEXX O 0: 4 eè 
LL UU LL шш -4 мэ = 
ZZ — C e° < ~ 
XN++ WY OPUS $e 
— ш ы г 
+ FZZ + (0 е+ са 
>> e e < 9% 


эгээ сс Ozxzaoqo 
OC CC UJU L2 E- ZI «x uui 


LUAIN eC) Ен ОЈ ө » 
«Tuo» Ж д сч 

ці- SHS || СУ 
НН и! | 

zw wo (ч О 
>>>> Fr PF «ч -- 
Се рече m ке «СО 
ccuocunu. «u«u comu © 
INN Se se Sse (ко 


O 
f~» 


2721! 


<= 


м 
э- 


SL + XL 
“St + XL 


171 


/ Y 
VM 


anms, фид 
>Œ 
= <í 
а 


=Z мај 

>> c 
>> mex IZ 
сс Sc > 
xw an í > mc 
Ine <1 _= се Ш 
m=O 1 «220009 
ни << = vv 
< ^^ || ~ 
"IH H 
> ~ am j] 
ww) »22 
XY) есозжилсдчк ед 


emm ЕЕС, 


*e(5”3 el 


=) 


| СО ZZ С0 2 2 


22--—02202z2Zzx-—2z 


== 


VQNNZZNNZTIDO-< 
= 


85 


22! 


O 
O 


см 
N 


IFS1; 


~ 
со 
UE 
рл 
= 
= 
> 
— 
еті 
Y) 
ке4 
y 
= 
ш 
= 
> 
>< 
=) 
он 
~ 
e 
>< 
==) 
е! 
~ Ww 
QD» e 
DU 
и (О 
r 
~ ~ 
= СО 
NO 
UE 
rÀ 4 w 
~ сатты. 
- 
> 
ami © 
су си 
Nu ш 
~ 9 
Oo 
-- 72 
LL! = 
> 
4 
с и 
= = 


GC TO 10 


OUTPUT PART 


Qo 


Реј и кә) ғы) 9 0 Ф ~ 

= aA 
H Ww H H M 

H H H H IH 
а н i= 

е к кен ме 


> = » = 

em PRON (Ge Фъ о. %- фь 

ju pre) == gemi {== с а ор р 

Ф 9» o P^ ферита оф фен 

СУФ СУР но он о о 

— — uns (NJ (Y) T om 
we? wee “ee wee 


X XXU) 
2.2.2 >22 за 
5) =) e Ec 
лилии 


Lodi "a dh di dA d w 


Kem, UE, quu e каршы, SO, e, IT en a, at, 
m (NJ <T' (N + —4 (N (Q ST e 
ті өне т ОЈ (NI YY cO) r en «ft 
СМО СОЈ NAN 
~ ~ ~ ~ ~ ~ ~ е ~ ~ 


0000000000 


“эээ? “ээ” wat xww ` x w m ` Y 


шшшшшшшшшш 
H i e В В н Hl 


=ч эч ум — — — — — — — 
СЕЧЕ СУ СЕ СЕ СЕ СО СЕ СЕ СЕ 
= Z =Z = = = = <= = = 


цу 
өші 





PTE, qu, aue, Gr. que, ёлы лың „ашы, SUE. SE, нэ Gillin, Gh, EM Gl, GM, Gl, SE SER, SEE ғы олы, — жыш, ышы ss а. [a m уулны, 


КОЮ = е >72 — Ft + eB en ~ ~ ~ ~ - = ~ ~ 
то бо оци се оо боб бо дрен ри 46-0173732"2723727272 чо ро e ei em oH em ћир Јиг Дио Дар Дер Ди 72 70 em e 
үзэн qud md md m Se ated et et et esd 4 ~ ~ ~ ~ ~ ы» о ғ ~ ~ ~ ~ ~ ~ жө ~ 
Hou n on WoW aaa o os Lo Lo ||!) 
"uuu H lH f H H H HW JH J lH 
м ч mt fat pao И ПТ и и ти и и s f #H Jl И И и a a 
"кеі кі іе =i i= = =! = Se 
— (=; — i= = ee Qe бон] (m=) = i= (оној фина бан |) s=) 
> @ ~ ~ те. е & гот? е > > 
Ф% ;; Өз ~ (Om fp qa ~ a е оо бо дои ми сти PE) йш; БШШШ ¿Sm дина, а) y 7 
a a ss a фи Фе фи фе. Фе Фа. фе» ee et Sete Фе рь Ge Ge Фи m ul аф 
jum um aq ba mt ФП оф ud ju фит Џеј Џо dem em | em, ar uum, oem у фе. 99 0. Өз до Фани оғ ұлы, т. от м «ғғ, ° ф» 
р те АР он = = = ~ во ~ С (NI эе эн ——— —— CJ f ('N + (N АД iem N А 
САЈ СО T AN P ~ Вч ССП фа СУ СО F (N pt er er "о о Фе @ s о о ~ ~ ~ о ~ P ms ча 


-----ХЖ 7 XE—————————— СМЕО ISE ANNO ИСО Ё 
G GO < 422172 СО СО СО CO Z ` = >< >< Ç) (* — —— 0 — w ( Q) >< >< Ç) (/) —  — x v 5 
оО === свазссоо а 
=== EESTI) SORTIE зора со со со сО 55595222272 
DNNN млолжХхХХЕХХЇЕХ опор 


Ач чы O ga? с“ х“— х“ ыра” Sa! See Qaqa am a? a? Y” o? a? E a Ya? ga? "ga? mt s swt ac мб rums ss unes “Qm Ass N s wama А as ss 


gre GE) co, ¿ey ETT uua, qim, o, ou, qi. y GTO, GO «gr, GO GD GO GD GS GO, GO, q IT, GOTO, GT, GD us, GD, GO cr, GTO, GB, GT, GO cum, o, qn, ¿sum о Хо A се В дет. 
CN СО“ SQ OO CUT AUST SNOOP MOM TO QO ODO OOO OOD QO O CO O CO O OCO O O OO CO C O CO CO OCO PC OOO OO 
+ +f SL LAO LO —— cN CN ry са rQ + +T РО Ши О ОО ОО оо омиили ороо оо о оюн ~ 
МСМ ТТЫ OSEP SES SE SS SES NUI СМ М СС Мс СМ СМ СМ СМ СМЧ СМ СМ СМ ӨМ СМ сиси с СМ СМ ANIN 
oer ~ a ~ ~ ~ е ~ ~ v» ое Өз P» ~ 9: ~ ~ p» o ~ p.m sl то 9%. ~ т; ~ ~ qa ~ ~ ~ ~ ~ asma ое ~ ~ ~ & е е ~ 


0 ХО ХО ЧО 0 0 0 000 00 00 Ооо о Г КМ PD SEP e P P P RR RP 


о "ugs "ugs tags s sss sf y as vs sa “sss Vus a a? a? "a ugs sam qas Yama ss a a as a x a a ss ges s sama ws st ef Wq s a s eh 


шшшшшшишшшшшшшшшшшшшшшшшшшшишшшшшишшшшшшшшшшшшшш 
-- Н Б Б Б Б Б Б Б Б Ке 


кл |ә |ы Quent jr бан seed (ен fu rd Gun Fd PE feng ead Peed emt ed ey frend ed fed peg Pt fend d pant jmd md que d dud ча) р Сеа nq aq aq aq aq q anaq aq, aq d ue Dd 
иг cy CY er: c c си с ск ес сс ог ск СК с СС ек, ск ск СС ск СК гу Се сс се сс сү се Се с се со се сс ек Се Се се Се се ск Се Се Се 
ЕЕЕ Т ЕЕЕ 


222 





«лт, ұл, Ee eg сатир Ee eg Ey A, өв Хальт OR, OE uum, qm, qm qum, SD qm qum, SU ¿m қақ, олт qum, SEN) qm. cum 
a Хо дава, quM, GEN JUN. QU UH qum, uUo S II CTE Ето өпті PA ото) дано цин, Ста иа иы (лсе en, GP, I Sy GG 
~ —< — — `` 


тї АРААРАР АЕО ОЧ ЕГАР АРААРАР E ОО 


0 9 9 9 9 % 9 090 6 9 6 9 6 Ө Ө 6 9 0 9 0 6 6 99 99 9 9 9 
Сы СМ СМ СЫ СМ СМ СМ СЫ ЄМ СЫ СМ СМ СМ СМ COIN СМ СМ СМ СМ СУ СМ СМ СМ СМ СМ СМ СМ СЫ 
—— — — — чч кч em — — —— —— — — e m — — e e A 
LL. LL. U LL. LL. LL. LL. LL. LL. LL. LL. U. LL t. LL. LL. LL. LL. CL. L. LL. LLL LL LL. LL. LL LL. LL. LL. 
оэооооооооооооооооооооооооооо 
ә-езңі езі к ———= ---- emm A n m -- - - «--« 


~ ~ ~ o^» ^e» ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ^» ~ ~ = ~ > ~ 
>< >< >< >< >< >< >< >< >< >< >< >< > >< >< >< >< >< >< >< >< >< >< >< >< >< >< >< >< 
LAY UN UA LO LA LALA LA LO UY UY UY UY LY UY LC LAN LA LA UR UN LAY UR LA LA LN UN LO US 


"му? “эдэ? Stage? Sarg? Saag? ЧР" Чани” Tage “cag? tag Зэв? tent tag ME Чао ЧАР аР at Se? eee Tag eae Зээ? Чи” “зэр? Зар? “эм?” Зэ” tag 


5 ОА SS RS SS SR SSR АА FR, RS SR — мн 
МА ДА б Rk SR SR SR SSS < 
o SR Ree SS SR SR SS SR SS 
``` SS SS SS SSS SS SSS SS 


~ ~ 9 9 о. е е ~ ^ na ~ ~ Ф. ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ о. ~ ~ ~ Ф a 
ч = = >»  р-р е е е ое ое ш ое qp Wo ое o ш о ое Р ет е е ц е C9 


es cv) <г ex 
EI С ОС СС ЕС 
<í <í < < сее 
Hain st СО СО СО СО МАРСУ еч СЧС СО СО со осм! 
e «Too +T СЧТ 2.5. 2.5.5. 2 2.5. 
KK KM XXI = = > >` > > 
bm fm few fm fee fe — — A ји eb ee |- |- |- СО 


332 55 — — — —— ey 
> — > > — <r( <( <( <( <( <( < < > > D> > <t <( <( <t <( < <€( «x 
Оло а NONNZZZZIZIZ 


= = SO pu, port od in — үрээд зі 


> e -” ы - ы о ғ ғ ғ % о. "е 
KXK е е ерес > e = e> >< 
+ + X< >< >< T —+ + +T О ТТ 
О O + -T-T Q О QO O O О ЗОО О О T T.T Хо 0 
400% 40400 
MÍN х х ~ 
° * e 5... ~ ~ ~ ОД ч ч ч ч ч ч, ч ~ ~ ~ ~ ~ ~ ~ Т. 
OO Q9 CO ew ww, , A A чо чено АСК А ФАСО Эм М 
це. ``. ``. ``. “чы “ы “ы ер “чы М ЧА, “ы 
ОФФО ММММ ах Мам ММММ М М МАМ х “ы, “ы, “м 
деде ~ ~ ~ че 7, ~ ~ ДО АД A о ЧАД ЧАО А во NNNNN 


о чац ЧИНУ ee Te? “пи “аар? “эээ? “тын? ЧУ "e? “ego? “aa? "GO aus O “e "O YE “ты? “тын? “тын? “тин? “тыу? "a? E “зар? “эээ” “нь” “чим” "a = 


--- RRR RRR ER RR ER RR ERR RRR EE EEE 
< <ç <( <( <( <í < <€( <Ç( ааааачаа<а<ааа<аа<чаааа<а<аа<а<а 
Mere ss SSS SS SSS TSS TS бо 4. 


EN CEDE CE OC сс сү Сү су гу Сү су су СУ oc oc ac oc СЕ СО Се СЕ ac ac oc oc oc oc oc oc 
ПИТИ ДОДООСОСОСООДЈОДООДОСОДООСООООСОООО 
C U LL LL LL LL Lu LL LLL LL LL 


СО) СЕТ СӘТ TFN MS DDN SOY SA NINOS e СА У Q CO 


DN OP е САСА 00 00 00 00) «pr ST. ST TF UND e em ICON AIM MMO TF T T ST uu 
NNNNA NAINN NNN AOA СЈСЧЈ СА СЈ СЛОЈ ТТТ ТТТ ХЕ ХР 


223 


БИДЕ А, ОТТІ 


e, 


e XK ПЦ ие = Ц. 
u^ ~ > WN 
"е ое е || um = 
пе x ll 522 
= Om е И Oe em 
T $ ed +->< ә 
NT eN AQ ALCALA 


See NKOOUM mu. 


чеч CUNT Le e ete 
Ne lle ANA NN cee 
ww е ~ 
w s (/У > ыг иг тээ зээл 
м< М >< < 
e жила 
|- |- - -- 
<< LAA 
т> 222206. 
СС б OA ACA 
Cc) со сл 2 
U uU LL ULL LL LLL CAO LLJ 
HHHH HH 
OO "A NU O ч 
Qo OOOO 
eJ МС СА 





) 


К, В 


( SUM, 


BFTAS 


SUBROUT INE 


> 
сї) 
E 
СО e 
— = 
Та e 
end — чиг 
me, (Y + 
© О «ас 
га ү" 
~ "m + OC 
aa О ~ С. 
(2 шом 
~ ] — + 
_ e | 
= d 
eng — ww Z 
= ~ > (00 
Dam IND 
Une 'мс ©) 
w с ап 
шош ut II 


n=. DE+Z 
Цј~ С) ча о > 


2 0 =O 

«оси с) Du 

OWN CL — O (22144455 
447 
че 


Е, EXP, Ls В, 150 , 11) 


( 


АЦТОК 


SUBROUT INE 


вать 
— 
=) 
~ + 
(N == 
p = 
CL a. 
— x > 
Q 89) WU 
O + + 
© ж = 
С) 
~ + | 
e 
22) > — 
= імеЧ 
— зае? 
> LL, u. 
e II T 
(A) a 
Qe ~ = 
С) ~ == = е 
e G. ~ ~ 
-— 2X ZU. L. 
OW 
>< > 
2 
ас» a 


- & 
ee o а 4 
a pa ОГ) 

Olm e — 0 
Qe | me 
Q4- 10 
Qa ое mo 


Ц. н С) Су MOM 
Du X w 
= о. = 
zat It I! 
С.с or 
mI == еі === 
zu Il 5 
C) ғ мейлі. 
«122... 
U O im e A UL O 
m COQ U OO 


224 


ш 

=> 
== 

ч СС 
- 2 
— 
Q tu < 
„сш 


= 
= 


19, к, 15, L, B, EXP, F, SUM, VAL, Il, U) 


( 


EARMA 


SUBROUT INE 


= 
СО 
чи? 
= 
= 
WY 
ge 
© 
O 
c 
сә 
“ Z 
ч”. Y 
Цц = 
= 
~ W 
= Y 
O Га! 
са = 
O 
O ~ С) 
~ - 
а - С) 
>x w OD 
uy = 
= 
- SY 
== ° 
СО UU e 
= = 
N e e 
с? 
ш 
= ~ • 
„JO 2 
O YY 
= | ~ || = 
LL! 
[AI А 
Ule «Lu. Ou. 
CC OO m» Mx mM 


) СС ТО 10 


SUM( IR) 


ӘСЕ? 





FCIS) * IQ 


-1Q)*EXP(1S+11) + 


PRI 


IRMI 
L*EXP(1S1-1) 
VAL + L*(1 


L*(1-T1Q)*EXP(1S+[1) + F(IS)*I1I0Q 


VAL 


X 


Сади 
pao. 
EU x 
Си — 
Y в > 9 
+ UNAS 
= шј чиј 
өт 2 || > e 
+ се 
Di gira И ас 
pi ч = 
II ы с Z 


|| ~ СМ < ~ ОС 11 СК 
a > 2 2 
a= be at R= 
“сис „wawz 
өсч metres) et OO 


<> сч цэ 
= N = 


225 





000 


SERV; ЧУВАК, ARIVM, SERVM, WMBAR 


ARIV, 


MB(5C0) 


? 


) 
W 


1 
500 
Ж 


C CY = 
$ UU 

DE 
wd LL LL 
Ir 
UJ c rm 
сс (3 


>K 


NNN 
се) м 
Ше 22 ш 
CY CE = 
сл 
> ~ LL > 
се) = > 
СЕ СҰ uy 
~ — 
ro > 
xXx ul - 
Co ск ш 
Y p<: 
> > > 
22 eu > 
Cie 2 
= 
OOOO m 
СО ОСОО © 
Mirate СУЖОЖ 
~ ~ ~ ~ »” Cro 
LA O uy LO "O + ~ ~ 
wee? Saggy VR аа? се . • 
22 Z= — = 
с ш ш 
> СОҺ-СОСО оо 
2 < << || -+ 
Сома сеща ща té x ces 
са Z í< £ > Z Z >< 


= 


~ ~ 


же 
шш 


~ > 
< 
шшш 


me ке 


— om) 
UN >x 


+ ++ 


pnl 
= =ч 
Zu r 
x 
A Uy LU 


|ке аш, Өл. 
ме == 
y == 
ч О) >С 
eL, LL 


226 


(5,304) 


READ 


WRITE (6,304) 


ж 


үеэ 

ж 

ш 

же і 

— e= "< (/) O 

Lili) — 

2 ~ O 

ши + + en 

.E++ 2 m + 
uc >> ч 

Кени о 


іне 


"n I ji J! 
e === 
C O J > > > >t Е 
A JODIO CY eet í Z kL) 
NOCZZXWXUSZ< 
сч.) > Z n < < ¿n Cn r = 


C 
с 





=> 


>> 
| СЕ 
CW m <) D > cx U, 
«T Unt- «t - 50111 

они а шие 
<= л = 

= (~ 0 


сред) с 


I| | 


MY) SKY) Tun 
EZ I Zu Dun 
=) О ~ ај.) ) 2. 2. ~ 2: сү 
UU Z Z=Z u < =Q O = 


3° 


WRITE (75,302) 


ЭХ, *RHOS 


/7, 


ТЕЗ 


eeu 
е e. 
“Ә 
eL. 
- х 
pre 


ТЕХ, 


12017) 


1 


ОСЕ ОУ 


5х, 


14: 


? 


72 


— > 


мл uU 


| = es 
mim 
> 2e 
e CO ccr СО 
CC UL. «ILL 
< = 


> е e 


о ж. қата ух» 
— =. 

ООШОТ ЛООГО 
тг ~ ~ ~ ~ a е a Ф% 
6 ӛг-ітті регі тағ) үлы үшеф ула) ра 
циа 
едее 008 H H H H H il 
LL. UL. 


ЦА ЦА i waq q bored pd toed font gem 


> ~ 

жж ~ Фа Өз ^» a фр ~ ~ 
LEV LD ee eee oe atte cay oo cay s 
Om dd Domed freed Pd fd me py waf 
“Зи Coen чад мои чин agi) Saggi Чаши? See 
~ У 0 >> 
+ ЕЛЕДОЖХООЋШ 


өт, 


со ° DD === 
æ e ил о 2075 
С) m ° ит 


a Sw» “A tom LL LO 


- х 


~ Пе ММ > хин 


О il 
— (Ni e 


N > 
Pa шг 


po 


> 


Fr 
ша 


== 


221 


> 


о N 


~ LL (ГУ ғә was eee ee See” “евр? о See Senet 
OWN ги. 
~ eun 


*QONZ ЗСО =O CN 

ғо сМ Ш >> Се o eam am m am роь пато бе дат 
(Tq am A NZZ Теч СО зС САЈ ССМ СМ М) СМ СМ! СМ 
Оп». СМ UD) e че СОИ ЧС СО 


ы. 


e 7% 


Дй 


ще (NCN CN CN C IIS Af + 
ча) ~" ~ © Фа Фа ~ %» ~ ое Фа 
9 O O QO 9-0 ФО 
ug Чани” “ым” Sat” “e “чай? "a? "a? 
eere 
<( < < <( LLJLLILLJLLILLJ JLL ШЈ 
Тэр ББ 
CK. CK CK. CK — r q — w w ме y— 
С: СОО СЛСЈЕС С СС ОС СС ЕССЕ Се 
Lu. LLL. LL LLL < = < = Z = = = 





FJ cu GE m HR жен, A, e. A, SP, ты, 
AE) GE E) лла со, жау оғы сқ ы 
~ Ва ч ~ ~ ~ ч 


СОУС O 


ОО е 9 9 ое 9 е ве е e 
CNCON CU C ысы сым сы ы 
== ener quu] pum еі реді gu] Gr quem] gm 
шешшш щщ LLE ЕШ 
ОООООООООО 
A 4 - сэн сэн ут 


re ~ ~ ~ ФО - ~ ~ 
>< >< >< >< >< >< >< >< >< >< 
LLY LEVELS EO IN LOLA LA LA 


bd d dA dA Adi A A A 


У ~ 
~ ~ 
~ ~ ~ O OS 


= дока, с а с у ss ¿s ms Qm 


ри оо о ә ==) 


~ qu ~ ^o ~ ^» e ~ ~ ~ > e 
que] quen] quond qued] ce] qu) queni cod eS фе 


jud eet juae mmt dum en emi jud рано pend 4—4 


~ ое ~ ~ ~ ~ ~ ~ ~ е ~ q 
С A Tey ES, ay Ry, Ey EB GTR У 
jud et aq: a aq Parmi pom dai jmd jud jud jud 
um wap о а и tassi "ugs? "ums чине чине сар Чаши” 
0 хә о Y cd 
>= > ZZ C ж2 5 >: 
= — зоо 

MY) = 
NY 


NSS ISSN SUN 


Ф ^» ^ ^» ~ ^"^ - 
eo =» = = гр е өт е е 


N с 


22 22 
< < 
NIO NNNON 
NN сэ 
Xu XN == 
FRO rra 
Zum 56 
ораа DIA 
NNZZ UNT 


~ ~ > > ~ > & 
х «ххх «е«ххх 
RT SR XS 


mem DIT OO TODD 

e CO US 

9 9 o". — — — “мч зм “м “м См. 

CO CQ CO — “ч, — — — — — — — 
пио и LL LL < O М М — М 
CN QG Q O C2 QO O OQ OQ Ç? G Q QO O — = ese СК 
LALA LAN O О 1 0D at et ee ~ 
FEN AIAN СМ СМ СМ СМ СМ СМ СМ 
тео ~ ос бр по 0» P ~ фе» pp "e "el q e” wn e. ` w tef we” 
OO pe Pe fh m 
"am waqa Se? tee Чаши” “Чаши” ey Se eee ee See eee 

misi Б Е 
LES LC LI LI LU LL «X «C «T «T «X «r «c «1 «c «wt «T «t wx 


Fm fe fe þa þa pete xxx z z === = > o. 

p pant mt ped Pt Pd ot Ot СГ СС ОГ СК СС СС СС СК СС СК СС СС СС СО) СЗ 
асс се се Си си Се се се се се СуСЈ Са) СОС) СОС) СС) Су) С)уС) е 2 
Ес = > < — LL. LL. LL. LL. LL. U.. LL. LL. LL. LL. LL. LL. LL. ul 


СОС)? СО СЧ САЈ СЛСЈ ЧО СЈ САЈ СЈ САЈ ЧО 
VD CJUO ST US e СА СО ГО 
ММ ТТР 


228 


) 


Е, ЕХ, ВК, ІС, К, 5 


( 


AUTOR 


SUBROUT INE 


ж R 
ЕХ (К+ 1) 


U( 10900) 
S 
F(1) + 


rmja uu 


т ч— LL 


Ос с ~ 
zw jj 


= Шел J 
= > =4 J 
C mM clu. 
WOW Om 


F(1) + EX(I+K) 


FUE) 


) 





TIST OFTREPERENCES 


nett; V. "Convenient Probability Elctting Positions 
ШОШ пе Normal Distribution.” Appl. Statist. v. 25, No 
ШЕ. 47-50, 1976. 


Barnett, V. "Prckability Plotting .Mettkcds and Order 


seat istics." др ал Ноа Па опери на pjo М, р. 95-107, 
ISS. 

Beyer, W.H. CRC Handbook of Tables fcr Probability and 
Peeatistics, The Chemical Rubber Co., Ohic, 1966. 
Mecontrci Program-67,Cambridge Mena to System 
(CP-67/CMS) Users Guide", International Business 


Machines Manual GH20-0859-2, May, 1972. 


David, H.A., Order Statistics, John Wiley and Sons, New 
Yerk, 1970. 


e 1949-1965 Time Frame. м. 
Thesis, US Nava Fcstgraduate School, Mcntersy, 
1 


scu roer ria, 


Feller, W 
ВЕ са 1оп5, у. II, 2nd ed., p- 19-21, John Wiley 
ani Sons, New York, 1971. 


— e Se Domnii SE mi Шы» стреа = == Amm ub m mmm m = ЕАУ АЕ. 


Fishman, G.S. Concepts and Methods о 





=< aa «ты» «жы» vs ae es < — сж s cx 


IRSE. 


Gaver, D.P. and Lewis, P.A.W., (Sone LeCOMINngG rarer on 


229 





TU. 


11. 


125 


14. 


11512 


118) 0 


ЙГ 


incr, 


B2. 


further EARMA process results), 1978. 


Gaver, D. P. and Thompson, G.L. Erogramming and 
Probability Hdodels in Orsrations Research, Brooks-Ccle, 
Eonterey, Califcrria, p. 452-527, 1973. 


Hutchinson, D.R., "APL Histogram, Density Estimation 


S 
Оне е5 Thesis, D.S. 


Naval Pcstgraduate School, Monterey, California, 1976. 


and Prctability Elotti 


lo 
la 





Barter, A "Expected Values СЕ Normal Order 
Exnstics." Bilonetrika, v. 52, p. 151-1:7, 1961. 


Jaccks, P.A. and Lewis, pA "А Mixed 
Autcregressive-Moving Average Exponential Sequence and 
Point Process (GARMAT M) SN Advances аа Ла редова 
салату, у. 9, p. 87-104, 1977. 


Шгаго, T.L., "An Investigatio 





n € 
the Expcnential Mcving Average Pcint Erocess". M.S. 
Thesis U.S. Naval Fostgraduat 

ша а тсггја, March 1976. 


Lawrence, A.J. and Lewis, P.A.W., "An Exponential 
Moving Average Segence and Point Fıccess (E*A1)". 


Шола 1 Applied Probability, v. 14, p. 98-113, 1977. 


Lawrance, A.J. and Lewis, P.A.W., (fcrthcoming parer 


on furtker EARMA process results), 1978. 


Lewis, E.A.W. and Cox D.R. "A Statistical Analysis of 
ШОС Опе Circuit Error Data." ТЕБЕ Iransacticns on 


ЕР эр» ч» m > > ар а A A A “Qu mb m вэ» 


Communications Systems, Com-14, p. 382-389, 1966. 


тес, G.R. "The Jacknife - A Review." Biometrika, Y. 
mr. 1-16, 1974. 


Pearson, ®.$. and Hartley, H.W. Bi 
Sra tisticianms, v. I, Cambridge 


Cambrzidce, 196060. 


230 





20. 


21. 


N) 


М) 


Robinson, D.W., “Non-Farametric Quant 


Naval Pcstgraduate Schocl, Monterey, California, June 
1975. 


Rosenblatt, 27 "Curve Estimates.” The Annals of 


M 
Hathematical Statistics, v. 27, 1956. 





Wilk, M.B. and Shapiro, S.S. "An Analysis of Variance 
Test fcr Normality (Complete Samples)." Biometrika, v. 
EE and 4, p. 591-610, 1965. 


231 





ENTIPTATDEDISPRIBUTION LIST 


Defense Documentation Center 


Cameron Station 


peexandria, Virginia 22314 


EXbrary, Code 0142 


Naval Fcstgraduate School 


Monterey, Califcrnia 93940 


Department Chairman, 


Code 52 


Department of Computer Science 


Naval Pcstgraduate School 


Monterey, California 93940 


Department Chairman, 


Code 55 


Department of Operation Research 


Naval Pcstgraduate School 


Monterey, California 93940 


Professor P. A. 


W. Lewis Code, 55Lw 


Department of Operation Research 


Naval Pcstgraduate School 


Monterey, California 93940 


Asscc. Professor G.G. 


Brcwn, code 52Bw 


Department of Computer Science 


Naval Ecstgraduate Schocl 


Monterey, Califcrria 93940 


Major Georgios I. 


Danikas, Hellenic Army 


Thrakis kai thermcpylon 4 


Agia Paraskevi, 


Athens, 


GREECE 


232 


No. Copies 








D14774 
ELT 





Danicas 

CMS histogram, densi- 
ty estimation and 
probability plotting 
routines, with an 
application to the 
analysis of the output 
of a simulation of a 
correlated queue, 











